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Summary. — The integral equations of the new Tamm-Dancoff method 
are derived with the aid of the Yang-Feldman formalism. The inhomo- 
geneous terms of these equations contain the boundary conditions at the 
time t = +00. 


1. — Einleitung. 


Zur Behandlung von Wechselwirkungsproblemen in der Feldphysik wurde 
vor kurzem von Dyson (') ein Gleichungssystem einzeitiger Wellenfunktionen 
(neue Tamm-Dancoffmethode) angegeben, das 4uSerlich den Gleichungen der 
(alten) Tamm-Dancoffmethode ähnelt, von diesen sich jedoch wesentlich in 
der Wahl des Bezugsvakuums unterscheidet. Ein anderes System (7%), das 
dem Dyson’schen System äquivalent ist, wurde aus dem lorentzinvarianten 
Differentialgleichungssystem (2) der mehrzeitigen Wellenfunktionen durch 
Gleichsetzen der Zeiten gewonnen. Beide Systeme haben, als Folge ihrer 
Beziehung auf das Vakuum der Gesamtenergie (?), gegenüber der (alten) Tamm- 
Dancoffmethode den groBen Vorzug, renormierbar zu sein. 


(1) F. J. Dyson: Phys. Rev., 91, 1543 (1953), im Folgenden mit A zitiert. 

(2) W. ZIMMERMANN: Suppl. Nuovo Cimento, 11, 106 (1954), im Folgenden mit 
B zitiert. 

(3) Für den anharmonischen Oszillator untersucht von W. HEISENBERG: Nachr. 
Akad. Wiss., Gottingen, Math. Phys. Klasse, 111 (1953) und K. SyMANZIK: Dissertation 
Gottingen (1954). 

(4) K. NisHsima: Progr. theor. Phys., 10, 549 (1953). 

(5) E. Freese: Zeits. f. Naturf., 8a, 776 (1953). 

(8) P. T. MarrHews und A. SALAM: Proc. Roy. Soc., A 221, 128 (1954). 
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Unklar bleibt jedoch in beiden Systemen, wie über verschwindende Energie- 
nenner zu integrieren ist, was sich als nôtig erweist, wenn beispielsweise in 
einer bestimmten Näherung des 2- Nukleonenproblems alle Wellenfunktionen 
auBer der |2|-Wellenfunktion eliminiert werden sollen (?). Dieser Nachteil ist 
darauf zuriickzufiihren, daB die Randbedingungen zur Zeit t= + co nicht 
explizit in die Gleichungssysteme eingehen. 

In der vorliegenden Arbeit soll deshalb ein drittes aquivalentes System 
angegeben werden, dessen Ableitung als wesentliches Hilfsmittel den Yang- 
Feldmanformalismus benutzt und keinen Zweifel dariber läBt, in welcher 
Weise die Randbedingungen fiir t= + co zu berücksichtigen sind. Die dabei 
erhaltenen analytischen Ausdriicke kônnen sehr leicht durch die Dyson’schen 
Doppelgraphen (8) veranschaulicht werden Der Formalismus wird nur in der 
untersten Näherung entwickelt, da eine strenge Fassung der verschiedenen 
Systeme einzeitiger Wellenfunktionen (in funktionaler Schreibweise) in einer 
Arbeit von SYMANZIK (*) erfolgen wird. 


2. — Analytische Behandlung. 


Wir leiten zunächst zwei Integralgleichungen ab, welche die einzeitige 
|2|-Wellenfunktion (!°) und die einzeitige 1 |2|-Wellenfunktion (!!) untereinan- 
der verknüpfen, falls die nàchsthòheren Wellenfunktionen vernachlässigt werden. 
Beide Gleichungen werden sodann in eine Integralgleichung für die einzeitige 
|2|-Wellenfunktion allein umgewandelt. 

In der Absicht die einzeitige |2|-Wellenfunktion 

Pil | X1X2|) = (Q, Lp (21) pes) Px 1 = (9, pla) p(@2)®) 2-4 (12) 
eines Zustands ® zweier gebundener oder gestreuter Nukleonen durch die 
einzeitige 1|2|-Wellenfunktion 


Ply | 152 |) = (Q, pas) p(@2) A (y) Poe = 10 = ¢ (43) 


auszudrücken, multiplizieren wir zunächst beide Seiten der Yang-Feldmanglei- 


) Zur Diskussion dieser Schwierigkeit vgl. A, Abschnitt IX. 
) F. J. Dyson: Phys. Rev., 82, 428 (1951). 

9) K. SyMANZIK: im Erscheinen. 
) 
) 
) 


gangselementen einzeitiger Operatorenprodukte vgl. Formel (7). 
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chungen 


t t } 
(La) (0) = y (as) + dg li S(a,— x,)y5A (a,)p(a7,) dar, + dm | Sla @,)y(a,) de, (1°), 


t t 
ie 


S(a—-0,)y,A (ata) dal 4- dm | Sla ah) pa) dx, 


— (28) 


(1.5) — y(@2) = y,, (#2) + ig 


8— 


gliedweise miteinander. Von der Richtigkeit der dann folgenden Beziehung . 


(2) (a) (2) =. = Win (V1) Yin (Le) 5 


t 
a i dé | dx, dx, S(a,— 2,)8 (2 — a, )y*{gylp(a,)p(w,)A (21) — tdmy(a,)p(a,)} + 
SIC 4! " 
t 


+] aw | da, dx, Sdi a) S(a— «yf gy? y(a,)y(a,) A (a,) — idmy(x,)p(a,)} 


— © €: 


überzeugt man sich am besten durch Einsetzen der Yang-Feldmangleichungen 
für y(€,) im ersten Integralterm und ¥:(@,) im zweiten Integralterm. Nach 
einer einfachen Umformung erhält man das Produkt der rechten Seiten von 
Gl. (1.4) und (1.5). Von (Gl. 2) bildet man die Matrixelemente Q + @: 


+ 
Polia) = (ia) + far fax dx, Sa(2,— #,) Sp(4,— &, )y?* 
— oo he 


)— idme(| LIA )}—symm. (14). 


{9yp(a | 21%, 
Hierin beschreibt 


(4) Py (lt) = (9, y, (my, (02) Die 


(13) S(z) gehôrt zur experimentellen Masse m. Man kann auch von der Yang- 


Feldmangleichung 
t 


p(a) = | Sm, (€ — a! y; A(x';y(x') da’ , Mo = M + om 
ausgehen, doch ist die Aufspaltung von m, in experimentelle und Zusatzmasse für die 
Diskussion der Randbedingungen praktischer. 

(4) Der mit symm. abgekürzte zum angegebenen Integralausdruck symmetrische 
Term entsteht durch Vertauschen der Indizes 1 und 2, so ist z.B. gy(x1|41x3|) durch 
Py (ay |x2x,|) zu ersetzen. 
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die beiden einlaufenden Nukleonen bzw. verschwindet, wenn ® Grundzustand 
zweier gebundener Nukleonen ist. 

Ganz ähnlich gehen wir vor, um die 1|2|-Wellenfunktion p(y |x,x,|) zu 
den nächsthôheren einzeitigen Wellenfunktionen in Beziehung zu setzen. Man 
drickt zunächst das Operatorenprodukt 


P(L1)P(V) A(Y) ea 


durch die einzeitigen Operatorenprodukte 


pla)p(@)A(YW)AY2), — pla)p(@)AY)AY2), TA p(y), p(y)] 
aus, indem man die drei Yang-Feldmangleichungen (1.a,b) und 
: 
(5) MEA | Avy—Yy LPY)s vs dy" (5) 


miteinander multipliziert. Die gesuchte Operatorbeziehung wird damit: 


(6) VICAUICAENUIF EEE = Win (01) Pin (2) A, (Y) su 


t 
ar if ae fax da, dy’ S(—2,)S (2 — 2,)A(y uni y{0x(x,, 525 DE 0.(x;; ot) 
— © Pi 


t » 
— if av’ fax: dx, dy’ S(a,— a, )S(a,— «,)A(y — y')Ox(x;, y', 1) 
— 0 +! 


0, = + igyly2y(a,)y(@,)A(@ JAY )oraveae + igviviu (a) pe) A (Yair ur + 
+ my + (AP) AY Jorn yee 
On = Di gviviv(a (e FU): PVT pane 


die dadurch nachzupriifen ist, da& im ersten Term von 0, y(æ,) und A(y'), im 
zweiten Term von 0, y(@,) und A(y’) und in 0, p(x,) und (ay) durch die rechten 
Seiten der Yang-Feldmangleichungen ersetzt werden. Von Gl. (6) hat man 
dann nur noch die Matrixelemente Q+ @ zu bilden, und die entstandenen 
Matrixelemente einzeitiger Operatorenprodukte nach der Wick’schen Regel in 


(18) Das Glied mit der Renormierungskonstanten dx? ist weggelassen, da sich diese 
in der behandelten Näherung ohnehin zu Null ergibt. 
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einzeitige Wellenfunktionen zu verwandeln, z.B.: 


(7) (9, (di) (Le) À (y1) A (Y2) D), Pa 


= lim (Q, Ty(a)y(as) A (y1)A (Yo) ®) 


CET ETES N 
0 40 
Lys Yj 


= Pil ¥1¥2 | ¥1%2 |) + A (y aii V2) V1 | 212% 


Ja 


Bei Vernachlassigung aller Wellenfunktionen aufer der |2|-Wellenfunktion 
auf der rechten Seite folgt so die Gleichung 


+ © 
(8) DAY | 2145 = p, °(y| aa: 1) + ig [ar | dx, dx, L(y TR 4; oy dv) 
i -— © ue 
Visa &)¥2Ru(| x, |) +- 
+ © 
+ ig far [ax: dx, Lay — 2,3 C:— 2) (tr — xy: |x,2, ) 
— © +’ 


[ Lp(213 2) = At(— &)Se(%e2) + Ax(€)S*(22)] 
welche zusammen mit Gl. (3) für Pil | %1%2 |) die Integralgleichung 


n 


max) = 97° (| eee) + | dt’ À ax ax, Kyle Sur — a) pe(læx |) + 


=n 1! 
O n 
+ | ar | da, dx, Kofi; £,2,)9(| xx, |) — symm. 
- © di 
159 i 
K,(æ, —- A) == — 10m8,(a,—- x.) +9? | dq 5, (œ — q)y;Lr(Q — lee q— y} 
+ © 


(9) 
| K,3(%1%2; 2,0) = g? | dg Sri — q)ysLr(9 — da La —- x, )¥2Sn(q— æ,)y} 
ergibt, wenn 9; °(¥|x,%.|) = 0 gesetzt wird (15). 


~ 


(6) Zur Diskussion der Randbedingungen für t= — co vgl. den letzten Abschnitt. 
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3. — Graphische Behandlung. 


Die Rechnungen des letzten Abschnitts lassen sich wesentlich schneller mit 
Hilfe der Dyson’schen Doppelgraphen (171%) durchfiihren. Im Gegensatz zu 
den gewohnlichen Feynmangraphen enthalten Doppelgraphen für jede Teilchen- 
sorte zwei Arten von Linien: Einfachlinien, denen S+, A+-Funktionen ent- 
sprechen, und Doppellinien, denen retardierte Funktionen entsprechen. Die 
in der vorliegenden Arbeit zur Ubersetzung von Doppelgraphen in analytische 
Ausdrücke gültigen Regeln sind folgende: Einer inneren Nukleonendoppel- 
linie (stark ausgezogen) und einer inneren Mesonendoppellinie (stark gestri- 
chelt) entsprechen bzw. die Faktoren iS,(#— x), —i1,(x— x), wenn die 
Linien von x’ nach x gerichtet sind. Einer inneren Nukleoneneinfachlinie 
(schwach ausgezogen) und einer inneren Mesoneneinfachlinie (schwach ge- 
strichelt) entsprechen bzw. die Faktoren — iS+(#— x'), î4+(x — x'), wenn die 
Linien von a’ nach x gerichtet sind. Ist eine von x nach +’ weisende Mesonen- 
linie bei x oder x’ durch einen Punkt markiert, so ist für sie — 1A, (4 — a') 
oder —-14,(%— a’) baw. i4+(@ — a') oder iA+(x— x") einzusetzen. Die Einse- 
tzungsregeln für 1|1|1-Eckpunkte, |1|1-Eckpunkte (durch 6m in der Zeichnung 
von den Nukleonenunterteilungsstellen unterschieden) und Nukleonen- oder 
Mesonenunterteilungsstellen sind dieselben wie bei Feynmangraphen: Einem 
1|1|1-Eckpunkt entspricht der Faktor —- gy;, einem 1|1|-Eckpunkt der 
Faktor idm, einer Nukleonenunterteilungsstelle der Faktor y,, einer Mesonen- 
unterteilungsstelle der Faktor + 1 oder — 1, je nachdem, ob die Stelle ober- 
halb oder unterhalb durch einen Punkt markiert ist. Die den auBeren Linien 
ensprechenden Operatoren oder Wellenfunktionen gehen aus der Beschriftung 
der Zeichnungen hervor. Sollen mehrere Punkte die gleiche Zeitkoordinate 
erhalten, so sind sie in gleicher Hohe gezeichnet. Sollen aber die Zeitinte- 
grationen in anderer Weise beschränkt werden, z.B. daB die Zeitkoordinate 
des Eckpunkts x, stets früher sein soll als die Zeitkoordinate des Eckpunkts 
æ,, so wird das in der Form 2° = ee ausdrucklich angegeben. 

Fig. 1 zeigt die Yang-Feldmangleichungen (1a), (1b) und (5). In Fig. 2a, 
wird durch Multiplikation von Fig. la mit Fig. 1b die graphische Form der 
Gl. (2) abgeleitet. Ubergang zu den Wellenfunktionen ergibt Fig. 2b [GI (3)]. 
Fig. 3a stellt die einzeitige Operatorgleichung (6) dar (1°), aus der bei Ver- 
nachlassigung der hòheren Wellenfunktionen Fig. 3b [Gl. (8)] folgt. Durch 
Einsetzen von Fig. 3b in 2a ergibt sich die Integralgleichung Fig. 4 [Gl. (9)] 
der einzeitigen |2|-Wellenfunktion. 


(1?) F. J. Dyson: Phys. Rev., 82, 428 (1951). 

(18) Uber die graphische Behandlung von Integralgleichungen der Feldphysik 
vel. B, Anhang. 

(1°) Die dm-Glieder sind in Fig. 3a der Einfachheit halber weggelassen. 


Le fa _- ~ — 7 A LE we Ae 2 ” F A 
x te | 4 4 i i 
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ie 
AX2 AX. Xo Ney 
Bio MI 4 À AS 
\ 
N 
Y Yin A Ÿ 
1Y 1Y Y 
! 
o di + © 
i 
1 I X 
D Ain 2l¥,¥] 
X4 Xe Ye Xe Xi X4 X2 
= ee - $m = SYMM. + 
ni 
da dé Yin Yin Ay A v Yin 
Xi 7 Xe Xo 
€ i 1 ° — SYMM. 
Fig. 9 cla SG + Xi{òm Vote x $m 
\ 10 70 \ 
ARE SAGESSE NY ¥ 
x, XS x, X2 Xi x x X2 
© = + ; + Em — SYMM. 
\ 
SN 5 SE LR 
P Po P S? 
Xi xa Xi x 1Y X ip AY X X2 AY 
1 i ' ' 
© i = ELLA SES 4 + 
i} 1 \ | 1 x I 
' ' \ ' x ' 
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hs da A Yin Yin A TOR ATX à ie sb i 
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Ganz analog verfàhrt man, wenn man allgemein, für eine einzeitige Wellen- 
funktion beliebig vieler Koordinaten, die zugehôrige Integralgleichung angeben 
will. Die durch Elimination der h6heren Wellenfunktionen aus dem unendlichen 
) allein hat die 


Gleichungssystem erhaltene Integralgleichung fiir ;(|2,. 
Form 


(10) plaire |) == py ara) + DANS) ara, asma (20) 


worin auf der rechten Seite tiber alle einzeitig irreduziblen |2|-Doppelgraphen D 
summiert wird und A, den D zugeordneten analytischen Ausdruck bedeutet. 
Einen Doppelgraphen D nennen wir dabei |2|-Doppelgraph, wenn er zwei 
auBere einlaufende Nukleonenlinien (?’) und zwei in den Punkten 2, und #, 
blind endende innere Nukleonenlinien besitzt. Er heiBt ferner einzeitig irre- 
duzibel, wenn er nicht durch einen Schnitt konstanter Zeit in zwei |2|-Doppel- 
graphen zerlegt werden kann. Gl. (9) und Fig. 4 sind gerade die g?-Naherung 
der strengen Gleichung (10). 


4. — Vergleich mit aquivalenten Methoden. 
Wir stellen zunàchst den Zusammenhang mit den in A angegebenen Be- 


ziehungen her. Dazu transformieren wir die Gl. (3) und (8) mit Hilfe der 
Formeln (??) 


CE, = Vp'+m, w=Vk?+ x, np) = E;" (ap + Bm)] 


(1 
S(x) == (27)? 


for exp [t(px — n(p)£E;t)ly: , 


E(x,;%,) = At 2,)S(#,) + A(e,)S*(a,) = 4t(,)8*@.)— A-(@,)S—(#,) = 


il eg ul I : 
= A dp dk vi; exp [t(kx, + px:)] exp [— în(p)(oxti + E,t)]ya 


1 } È + 
P1{| X1%X2 |) 22 go ign | dpi dp;q( pip» |) exp [1(pix, + p:%:)]; 


x, 0 (27) 


| 1 
Ply mx) = 7 x)? exp [ ipa] dk dp, dp.9(k|p,p.|) exp [i(ky + pix + prxi)], 


eae (A 


(2°) Ausgenommen g; ”(|x,x, |) sind alle Funktionen y, ® 


siehe den letzten Abschnitt. 
(21) Für die beiden 4uReren Nukleonenlinien ist in Gl. (10) die Wellenfunktion 
p,y(|x1x3|) selbst einzusetzen. 
(22) Wir machen die Annahme, daf ® Energieeigenvektor ist. 


gleich Null gesetzt worden, 


ry 
; 
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in den Impulsraum: 


Gane 
IE QC m(p)E,— n\q)E): 


ia, 
[ak (yh ph. 4d |) +o(k |p, 4 —k))} — 
— 2ni6m5,(E — n,( p)E, — m(q)E:)(yi + y2)p(\pq|) [= iviv], 


(12) p(k p—k, q)) = °(k| p —k, q\)— 


‘a Ox dE — mp RE: + Ox) — mq)E) y p(|1pq |) — 
ati | 
— PIA mp ME MD + 0) —k, q+) 
k 


Es folgen sofort die Gleichungen (A.66), (A.67) und (A.68), wie auch Gl. (A.69), 
(A.70) und (A.71), die statt der Energienenner genauer 6*-Funktionen enthal- 
ten (2). Inwieweit die hier zwangsläufig auftretenden 6+-Funktionen durch 
14(6+ — dò) ersetzt werden dürfen, wird im nächsten Abschnitt besprochen. 

Um die in B aus den invarianten Differentialgleichungen der Wellen- 
funktionen gewonnenen Beziehungen (B.VI) abzuleiten, benôtigen wir noch 
die Integralgleichung 


+0 | 
(13) PAY | 2122 |) = pa [xx |) + igfae | dr, de, Ly a i A 
— o (A 


Sue ye (| ata, |) — symm. 


die analog zu Gl. (8) aus den Yang-Feldmangleichungen (1) und 
(14) Ay) = An(9) —$ fay y PY), Py DIys dy’ 


erhalten wird. Durch Anwendung der Differentialoperatoren (24) (0/et+fi+f,), 
A-+-xj auf auf Gl. (3), (8) und (13) erhalt man (zusammen mit den Ener- 
gieeigenwertgleichungen von @i(|%:%2)), ply | xX.) und PUY | 1% )) die 
Gleichungen (B.VI) in der hier behandelten Näherung (?°). 


(23) Die erw ähnten Gleichungen in A und B enthalten allerdings statt m überall my 
und dafür keine ôm-Glieder. Vgl. zu diesem Punkt Anmerkung ("8). 
(24° f = a(0/0x) + ibm. 
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5. — Diskussion der Randbedingungen. 


Die inhomogenen Glieder der Gleichungen (3) und (8) 


pr ral) = (Q, pa). (#2) P) = (0, Yolo) (12) D'(— 00) 

Pe? (¥ | ra |) = (2; AMP (0, VLAD 00) (29) 
geben (bis auf eventuelle Beteiligung gebundener Zustande, siehe weiter unten) 
über die Randbedingungen zur Zeit t = — co Auskunft. Um auch den Hin- 


fluB der Randbedingungen zur Zeit t = + co diskutieren zu kénnen, ziehen 
wir die aus den Yang-Feldmangleichungen 


+ © + © 


y(æ) — Pout (1 es in | Sas — L'yA(a' )y(x x’) dx’ E è om | Sasa — x (e) da”, 
. Lo 
A (y) = Aowt(Y) + af Any (Y-Y{YY), vly')lysdy’ , 


abgeleiteten und zu den Gl. (3), (8) analogen Integralgleichungen 


dt’ 


dx, dx, STE — 2,)S_ (de — 2,)y}: 


(15) (| %,x2|) = (A) 


ne | x,x, |) — iòmp.(|x,x,|)} — symm. 


+0 
(16) 9.(y¥|2,x,|) =@F °(y | a.%2|) + ig ar [ax dx, L, (y— 1. ay) 
CTS 4 
Bra 2,)yipe(| 212,1) — symm. 
Le 2p) IVA pel Ea 21), (22) DE A ,, (21) S+(22) 


pt? (|xi%2|) = (Q), (0) (0) D(+00)), 


pi (ya) = (9, ve), A (4) D'(+00)) 


(5) 2; ist das wechselwirkungsfreie Vakuum, vol). Ally) bedeuten die wechsel- 
wirkungsfreien Feldoperatoren (zu den experimentellen Massen). Zur genauen Defi- 
nition von @/(— co) vgl. B, Gl. (II). 


+ 4 LA 


LPO ARE ET A PA ba + Re 2er J'TE, À ATEN FES US PE y 
4 - > b. A af ai PS 
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hinzu, die im Impulsraum folgende Gestalt haben: 


ig | 
(27)? d-(E —m(p)8, — n(9)E): 


(17) v(\pq|) =¢t?(|pq|)+ 


fort IP—k,q|) + y°olk|p, q—k|)} + 


+ 2xidmé_(E — m(p)E,— no(q) E) (vt + y2)9(\pq |) 


(18) o(k|p—k, q|) = o*?(k|p—k, q|) + 
ETS (EB — mp—k)(B, x + or) — m(g)E)y'g(|pal) + 


Wr 


ug . 
hes cd Ô_(E — mp —k)E,_x — olq)(B. + ox) y?9(\|p—k, q+k|). 


Wx 


Wir vergleichen die Beziehungen (11), (12) mit den Beziehungen (17), (18) 


zunachst ganz unabhängig von der speziellen Wahl der Randbedingungen. 
Die Funktionen 


Pi” (| Aix 


e) 


lassen sich stets in der Form 


(19) Pe” (| %1%2!) = Dd Ge {PRA | [PP (| 2 |) — 212) [ent La 


(20) p#(y|arza]) = Do ea] — Pt] PU] az asz ve 


nach vollständigen und orthonormierten Systemen {g{(|2])} und {y{{y|)} 
von Energieimpulseigenlosungen der wechselwirkungsfreien Feldgleichungen 
(mit experimentellen Massenkonstanten) entwickeln.. Nach Definition durfen 
in Gl. (19), (20) nur die Lôsungen g,(}æ|), g;(y|) positiver Frequenzen einen 
Beitrag geben. Daraus folgt nach Subtraktion der Gl. (11) und (17), bzw. 
Gl. (12) und (18) 


= (E — m(p) Ey — mg) Ea) : 


ig 
27 


gi Pi RE) 


fepretk]p, q|)+y*g(k|p,g-k|)}}— 


— 2ri 6m d(E — n:(P)E, — n:(QE) (vi + vi)e(|pq |), 
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gt? (k|p—k, q|)—9-?(k'p—k, q|) = 


= — 79 5(B — mp — RE, + + ox) — miq)Ed y'v(|pa |) — 


Dr 


i 


(42) 


d Ô(E M(p k)E 3 — nq) Es + ©x))y"9(|P — k, q+ k |) ; 
k 


daB auf den rechten Seiten alle Glieder verschwinden, die nicht d-Funktionen 


der Form 
0ET—- E E.) bzw. dE-E—- E,— 0) 


enthalten. In den Integralgleichungen (11), (12) dürfen also alle diejenigen 
ò+-Funktionen durch 1,4(6+—0-) ersetzt werden, die nicht Argumente der 
Form 

E-E-E., E — E, — E, — «w 
besitzen. Das bedeutet, daB in der durch Einsetzen von (12) in (11) entstan- 
denen Integralgleichung für g(|pq|) bei Integration über Energienenner der 
Form 


E-E+E,+@ oder HE—EH,+#H,+o oder E+E,+E,+% 
der Hauptwert genommen werden kann. 


Der Einfluf der Randbedingungen für t = + co oder — co soll nun im 
Einzelnen diskutiert werden. Wir unterscheiden drei Falle: 


1. ® sei der Grundzustand des Deuterons. Dann ist 
D'(— ©) = D'(+ co) = 0 (28) 


und es folgt 
PE” (|21X3 |) = PE” (y |x1%2 |) = 0 


die 6*-Funktion kann in diesem Fall also generell durch 14(6+— 6-) ersetzt 
werden. 


2. ® beschreibt die Streuung zweier Nukleonen der Wellenfunktionen 
A(læ)) und g;(|2|), d.h. ® hat die Form 


I 


D = U'(0, — c0)D'(— ow) 


DT co) = | ae, dx, D(a) (AP |)@ol| x |) DE Po(| Xr AE 2) . 


t 


(5) S. FuBINI und M. Cini: Nuovo Cimento, 10, 1695 (1953). 
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Ks ist: 


p, (| *%s )= P1(| #1 |)Pa(| Le |) — po x |)@1(| a> |) - PLY da ) = 0 
Im allgemeinen vereinfachen sich die Funktionen g*° nicht weiter, falls aber 
die Erzeugung eines reellen Mesons mit Energie- und Impulssatz nicht ver- 


traglich ist, folgt 
pi (ye) = 0. 


In diesem Sonderfall darf also in Gl. (12) d+ generell durch 15(d* — 6-) ersetzt 
werden. 

Als letzte Môglichkeit erwähnen wir drittens den Zustand ®, der die 
Streuung zweier Nukleonen mit den Wellenfunktionen @;(|4|), g:(|#|) und 
eines Mesons der Wellenfunktion ®3(4||) beschreibt. In diesem Fall ist 


Qi (y | X12) = Ps(y| ){qx(| di |)@o(| do |) — Po x. |)q1(| 22 |)} 


— © 


CRIE SE ON 


Die weitere Môglichkeit der Streuung eines Mesons an einem Deuteron im 
Grundzustand wird erst erhalten, wenn im Gleichungssystem der einzeitigen 
Wellenfunktionen auch y,(y:7:|x.x,|) berücksichtigt wird. Die allgemeinsteu 
Zustände, die in der behandelten Näherung des Zwei-Nukleonenproblems, 
nämlich Vernachlässigung aller einzeitigen Wellenfunktionen auBer der |2 
und 1|2|-Wellenfunktion, überhaupt in Betracht kommen, gehen durch Linear- 
kombination aus den besprochenen drei Moglichkeiten hervor. 


RIASSUNTO (*) 


Con l’ausilio del formalismo di Yang-Feldman si derivano equazioni integrali di 
funzioni d’onda monotemporali, i cui membri omogenei esprimono le condizioni al con- 
torno al tempo t= + co. 


(*) Traduzione a cura della Redazione. 
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Su di un’onda di pressione nei liquidi 


A. CARRELLI e G. BRANCA 


Istituto di Fisica dell’Università - Napoli 


(ricevuto il 18 Marzo 1954) 


Riassunto. — Se un liquido in cui si ha un sistema di onde ultrasonore 
stazionarie è attraversato da un fascio luminoso in direzione normale a 
quella lungo cui ha luogo la distribuzione di pressione, e quindi dell’indice 
di rifrazione, si osserva — come è ben noto — ad una certa distanza una 
serie di linee chiare ed oscure (HIEDEMANN). Sovrapponendo ora alla 
variazione di pressione periodica, prodotta dalle onde ultracustiche, 
un’altra variazione dipendente da una qualunque altra causa, ed avente 
la stessa direzione della prima, si osserva che tali righe si spostano paral- 
lelamente a loro stesse. Facendo scorrere una pellicola fotografica nella 
zona in cui si osservano le strie, si può misurare sui successivi fotogrammi 
lo spostamento di una generica riga luminosa (o oscura). Se le direzioni 
di propagazione degli ultrasuoni e dell’onda di pressione sono normali 
fra di loro e normali alla direzione del fascio luminoso non si osserva alcuno 
spostamento delle frange. Con tale dispositivo si può dunque rilevare la 
presenza di un’onda di pressione in seno a un liquido, mediante un me- 
todo ottico. 


Com’é noto, quando in un liquido si propaga un sistema di onde ultra- 
sonore stazionarie, si ha in esso una distribuzione di pressione variabile col 
tempo e con la posizione, distribuzione che produce a sua volta una variazione 
di indice di rifrazione, che segue la stessa legge. Se un fascio luminoso attra- 
versa il mezzo in. cui ha luogo tale distribuzione di indice, l’onda, all’uscita 
del mezzo, presenta una «modulazione di fase » che dà luogo ai ben cono- 
sciuti fenomeni di diffrazione; e cioè si hanno strati del mezzo alternativa- 
mente condensati e rarefatti che generano le note manifestazioni. 

Lungo l’asse ottico del sistema cioè, a partire da una certa distanza, sono 
visibili mediante una lente di ingrandimento, o anche ad occhio nudo, con 
l’ausilio di un vetro smerigliato, delle strie luminose che occupano tutto il 
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campo. Tali strie, parallele alla fenditura del collimatore, costituiscono una 
sorta di reticolo, che acquista la massima nettezza in una determinata posi- 
zione. Il processo formativo di tale reticolo è, com’é noto, il seguente: i raggi 
luminosi incidenti, supposti paralleli fra di loro e normali alla direzione x di 
propagazione del fascio ultrasonoro, penetrano nella regione che possiede una 
rifrangenza variabile ed il calcolo (!) dimostra che pur essendo il punto di 
convergenza variabile, essi convergono per un tempo più lungo durante tutto 
il periodo d’oscillazione in una certa zona, dove appunto si formano le strie 
luminose di cui si è parlato. 

Supponiamo ora che alla variazione di pressione periodica in t ed in x si 
aggiunga lungo la stessa direzione x 
una variazione di pressione dipendente 
da una qualunque altra causa. Si può pen- 
sare, ad esempio, che la nuova onda di 
pressione sia rappresentabile schematica- 
mente mediante la curva in fig. 1 a) e che 
essa si sposti lungo æ con la velocità di 
propagazione delle onde elastiche. Dunque 
in un certo istante la distribuzione di p Fig. 1. 
per effetto dei soli ultrasuoni è data dalla 
fig. 1b), e per la presenza contemporanea dell’onda di pressione, la distribu- 
zione della pressione sarà data dalla somma delle due curve riportate. 

Si suppone, come si ricava dalle due figure, che l’intervallo Aw in cui si ha 
la variazione di pressione sia dello stesso ordine di grandezza della lunghezza 
d’onda À degli ultrasuoni. 

Quando le due onde di pressione coesistono, durante il periodo di tempo At 
che l’onda Ap della fig. 1 impiega per propagarsi lungo l’asse x della vaschetta, 
è da prevedersi una variazione dell'andamento delle striscie chiare e oscure 
e più precisamente queste righe per effetto della variazione Ap debbono spo- 
starsi. In altre parole tutto avviene come se l’unica discontinuità nella distri- 
buzione dell’indice di rifrazione, dovuta all’onda di pressione, producesse local- 
mente un effetto analogo a quello che si è riconosciuto nel caso dell’onda ultra- 
sonora; vale a dire che i raggi luminosi modificano la loro traiettoria nella 
zona di discontinuità. Ma poichè le due onde si propagano simultaneamente, 
all’incurvamento dovuto all’onda ultrasonora e che produce il reticolo, si 
aggiunge un incurvamento dovuto all’onda di pressione e che determina in 
un certo istante lo spostamento delle frange. Poichè, come si è detto, la discon- 
tinuità dovuta all’onda di pressione viaggia con la velocità dell’onda elastica, 
in un tempo brevissimo e cioè nel tempo At dato dal quoziente del tratto { 


od 


(1) R. Lucas e P. Brquarp: Journ. de Phys., 3, 471 (1932). 
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lungo il quale si propaga l’onda e della velocità V dell’onda stessa, tutte le 
successive frange saranno interessate, per modo che alla fine il reticolo ap- 
parirà spostato. 

Noi abbiamo sperimentalmente messo in evidenza questo processo di pro- 
pagazione dell’onda di pressione nel liquido col dispositivo rappresentato sche- 
maticamente in fig. 2: V è la vaschetta 
che porta nella parte inferiore il quarzo @ 
e lateralmente due finestre per consentire il 
passaggio dei raggi luminosi. Il mezzo con- 
tenuto in V può porsi in comunicazione — 
mediante il rubinetto a due vie À — con il 
ramo 7, oppure con il ramo 7, di un tubo 
ad U. Il livello del liquido in 7, è uguale a 

Fig. 2. quello che esso raggiunge in V, per cui, sta- 
bilita la comunicazione tra V e T,, il li- 
quido stesso risulta sottoposto alla pressione atmosferica. In 7, invece vi è un 
dislivello di una quarantina di cm, rappresentato da una colonna di Hg, per 
cui la pressione addizionale che si propaga nel mezzo, una volta posti in co- 
municazione V e 7,, è di circa 1/2 atm. 
La velocità di tale propagazione è evi- 
dentemente la velocità di spostamento 
delle perturbazioni elastiche e cioè, nel 
caso specifico, trattandosi di acqua, tale 
velocità è di 1399 m/s. Tenuto conto 
che la zona di osservazione della va- 
schetta è di 3 cm, il tempo At di propa- 
gazione del Ap (e quindi della variazione 
di indice ad esso connessa) è dell’ordine 
di 10-5 s, mentre il periodo degli ul- 
trasuoni è dell’ordine di 10-* s. Questo 
significa che, mentre il Ap attraversa la zona di osservazione, il reticolo si 
è creato e distrutto circa 10 volte. 

Abbiamo potuto accertare, guardando il reticolo, che effettivamente nel 
periodo in cui si propaga l’onda di pressione è nettamente visibile il previsto 
spostamento delle frange. Lo spostamento si osserva anche nella fase di de- 
compressione ed ha luogo naturalmente in senso opposto a quello che si ha 
nella compressione. 

Si è voluto poi ottenere una prova obbiettiva che permettesse una misura, 
e si è proceduto nel modo seguente: si è fatto scorrere una pellicola cinemato- 
grafica, nella zona in cui il reticolo acquista il massimo di nettezza, alla 
velocità di 8 fotogrammi al secondo e si è poi misurato su di essa lo spo- 
stamento di una generica frangia: si ottiene in tal modo il grafico di fig. 3 


a 


Spostamenti di una frangia — 
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in cui sono chiaramente osservabili le due fasi di compressione e di decom- 
pressione. 
Acclusi si riportano gli ingrandimenti di tre fotogrammi, che mostrano le 
successive posizioni assunte da una frangia rispetto ad un segno di riferimento. 
A conferma di quanto precede si è eseguita una seconda esperienza in cui 
le direzioni di propagazione degli ultrasuoni e quella dell’onda di pressione 


sono normali fra di loro e normali alla direzione del fascio luminoso. Com'era 
da aspettarsi nessuno spostamento delle frange si verifica in tali condizioni. 

Infine si è modificato il dispositivo illustrato nella fig. 2 abolendo il rubi- 
netto À e producendo l’onda di pressione con l’innalzare un tubo flessibile 
pieno di Hg collegato al ramo 7, e ciò allo scopo di verificare se l’effetto potesse 
dipendere da uno scuotimento meccanico del sistema. Anche questa prova 
ha pienamente confermato lo spostamento del reticolo ed ha permesso inoltre 
di stabilire un punto fondamentale, che cioè l’effetto dipende non solo dal 
valore di Ap ma anche da At e cioè dall’intervallo di tempo durante il quale 
si ha la variazione di pressione. 

Questo dispositivo si presenta dunque particolarmente indicato quando si 
voglia misurare l’andamento di un’onda di pressione in un liquido, essendo 
sicuri con tale metodo che non si ha aleuna inerzia dipendente da dispositivi 
Meccanici (indici, leve, ecc.). 


APPENDICE 


Nella ipotesi fatta (fig. 1) la variazione della pressione è una funzione f(t) 
data in un intervallo infinito e si è supposto che l’andamento di Ap che è col- 
legato alla variazione dell’indice di rifrazione dalla relazione 

1 8n° An 


Eee *), 
De ORS eee Synge 


(*) R. Lucas e P. BIQUARD: op. cit. 


39 — Il Nuovo Cimento. 
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essendo « il coefficiente di compressibilita ed », l'indice di rifrazione, sia quello 
di fig. 5. Tale f(t) può esprimersi mediante integrale di Fourier. Definiremo una 
nuova funzione y(t) mediante le posizioni (fig. 5) 


PASO jit = ts PAU per eile aan 
La g(t) sviluppata in integrale di Fourier da 
Di x | C(w) sin wt do , 
0 


i cui coefficienti possono calcolarsi facilmente e sono: 


2K: 
(1) C(w) = — sino, 
O) 


essendo K=tg x. Può dunque dirsi che nell’onda di pressione discontinua ci 
siano onde progressive di tutte le pulsazioni ©. Più particolarmente le ampiezze 
C(w) sono date dalla (1) che è 
graficamente riportata in fig. 6: 
si vede che hanno importanza 
solo le onde di pulsazione piut- 


CW 


Fig. 5. Fig. 6. 


tosto basse e quanto più grande è K= Ap/At le più alte pulsazioni diven- 
tano più efficaci, 

Può dunque dirsi che nell’onda di pressione discontinua, e cioè del tipo 
riportato nelle figg. 1a) e 5 esistono onde progressive di tutte le pulsazioni @. 
In particolare le ampiezze C(w) sono date dalla (1) che è graficamente ripor- 
tata in fig. 6, la quale è stata ottenuta per 27 = 2-10-5s e nella ipotesi 
Re (a= 45°). 


Si vede che tali ampiezze assumono tutti i valori compresi tra -+-co e —co. 
Infatti 


lim C(0) = +o, lim C(w) = -— co. 


&—>0 +0 &@—>0— 0 


Inoltre i massimi e i minimi intermedi si hanno per 


(Qn + 1)m 


2T 
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ossia per valori del periodo che siano dello stesso ordine di grandezza del- 
Vintervallo di tempo in cui ha luogo la variazione di pressione (o multipli di 
esso). Le onde corrispondenti a tali pulsazioni si sommano con l’onda sinu- 
soidale generata dagli ultrasuoni, se, evidentemente le onde hanno lo stesso 
periodo. 

Questa semplice osservazione da ragione del fatto che producendo una 
variazione lenta di pressione non si osserva col nostro dispositivo lo sposta- 
mento delle frange; infatti in tal caso il periodo delle onde che compongono 
l’onda di pressione è più grande di quello degli ultrasuoni. 

Volendo osservare il fenomeno in tali condizioni occorrerebbe variare la 


frequenza del circuito oscillante sino ad «accordarla » con la frequenza delle 
onde costituenti l’onda di pressione. 


SUMMARY 


If a liquid in which there is a system of stationary ultrasonic waves is crossed by 
a light in direction normal to the one along which exists the variation of pressure and 
therefore of the reflection index, one observes, as it is well known, at a certain distance, 
a number of clear and dark lines (HTEDEMANN). Superposing now to the periodic 
variation of pressure due to ultrasonics another variation due to any other reason, and 
having the same direction as the first one, we observe that the lines deplace parallel 
to themselves. Taking a sequence of pictures on a film placed in the spot where we 
observe the lines, one can measure on the series of photograms the deplacement of 
any line (either a luminous or a dark one). If the direction of propagation of ultra- 
sonics and the direction of the pressure wave are perpendicular to each other, and 
both perpendicular to the direction ef light, one does not observe any deplacement of 
the lines. With such an experimental arrangement a pressure wave in a liquid can 
be detected with an optical method. 
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Sulla struttura elettronica del nitrobenzolo. 


M. SIMONETTA e A. VACIAGO 


Istituto di Chimica Generale e Analitica del Politecnico - Milano 


(ricevuto il 3 Aprile 1954) 


Riassunto. — È presentato un calcolo di alcune grandezze fisiche del 
nitrobenzolo mediante il metodo M.0. con una ragionevole scelta dei , 
AI parametri necessari. Per il momento elettrico, la frequenza del primo mas- 
Di‘ simo di assorbimento nell’ultravioletto, le densità elettroniche, gli ordini 
| di legame, gli indici di valenza libera e l’energia di risonanza si ottengono 
valori in ottimo accordo con le previsioni. Anche il comportamento cine- 
tico prevedibile nelle reazioni di sostituzione dei vari tipi concorda con 
i dati sperimentali: è valida la «chemical non-crossing rule ». 


1. — Introduzione. 


Nell’applicazione della teoria degli orbitali molecolari (M.O.) nella sua ap- 
prossimazione standard di combinazione lineare di orbitali atomici (L.C.A.0.) 

a molecole complesse contenenti eteroatomi, uno dei punti di maggior diffi- 

coltà è l’assegnazione di opportuni valori ai parametri necessari per tener conto 

della natura chimica degli eteroatomi presenti. Si è ricorso essenzialmente a 

metodi semi-empirici, ma anche questa via conduce a risultati incerti, soprat- 

| tutto per il fatto che il valore di alcune delle grandezze fisiche per le quali sono 
3 a disposizione dati sperimentali di confronto appare pressoché indipendente 
dai valori assegnati ai parametri necessari, purché scelti entro limiti ragio- 
nevoli (1). Recentemente JAFFE (?) ha sviluppato un metodo empirico per la 
determinazione dei migliori valori per questi parametri, nel caso di composti 
benzenici monosostituiti, basato sulla correlazione tra le costanti o di Ham- 


(1) H. C. Loneurt-Hiacins e G. W. WHELAND: Ann. Rev. Phys. Chem., 1, 138 
(1950). 
(2) H. H. Jarre: Journ. Chem. Phys., 20, 279, 778 (1952). 
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mett (*) e le densità elettroniche nei vari atomi di carbonio nel nucleo benze- 
nico. JAFFÉ ha preso in esame i seguenti derivati del benzolo: fluoro-, bromo-, 
cloro-, iodo-benzolo, fenolo, tiofenolo, anisolo, anilina e toluolo, ed ha stabi- 
lito per via empirica la seguente relazione: 


(1) Gr 16,6(4—1), 


dove o è la costante di Hammett per la posizione meta o para e q la corri- 
spondente densità elettronica. In tutti i casi considerati il gruppo sostituente 
contribuisce con un orbitale x al sistema in risonanza, mentre per il nitro- 
benzolo gli orbitali x sono complessivamente nove. Il nitrogruppo si distingue 
dai sostituenti considerati da JAFFE anche per avere un potere di attrazione | 
degli elettroni notevolmente superiore, in conseguenza del fatto che effetto 
induttivo ed effetto di risonanza si sommano. Sono quindi da aspettarsi varia- 
zioni di concentrazione elettronica sensibili, il che pone in dubbio la validità 
di una delle ipotesi di Jaffé, cioè che, ammesso che o sia una funzione mono- 
tona di Aq, questa funzione si possa sviluppare in serie e i termini con potenza 
di Aq superiore alla prima siano trascurabili. Inoltre, se Ag non è piccolo, 
l’integrale di Coulomb dell’atomo sede di tale variazione non è più indipen- 
dente dalla variazione stessa. Un altro fatto che fa apparire poco probabile 
l’applicabilità del metodo di Jaffé al nitrobenzolo appare evidente conside- 
rando i valori di g ottenuti dai valori sperimentali o, = + 0,710 e o, = 
= + 0,779 mediante l’equazione (1): qn = 0,954 e gq, = 0,950. Mentre cia- 
scun valore preso a sè potrebbe essere accettabile, la loro differenza sembra 
troppo piccola per render conto dell’effetto decisamente meta-orientante 
del gruppo NO, nelle reazioni di sostituzione elettrofila degli atomi di 
idrogeno del nucleo benzenico: da questo punto di vista miglior risultato ave- 
vano ottenuto WHELAND e PAULING (*): q, = 0,61; Qn = 0,95; q, = 0,79. 

Abbiamo comunque eseguito alcuni calcoli orientativi. I parametri in gioco 
sono l’integrale di Coulomb dell’atomo di azeto e di ossigeno (x, e x;) e l’inte- 
grale di scambio per i legami azoto-carbonio e azoto-ossigeno (yy, © Yxyo) 
Per applicare il metodo di Jaffé abbiamo in un primo tempo fissato due di 
tali parametri: 


go dee Ce | 


(in unità BP, dove f è Vintegrale di scambio C—C nel benzolo). 

Il primo di questi valori corrisponde a quello trovato da JAFFÉ nel suo 
trattamento del fenolo e dell’anisolo, mentre il secondo si uniforma all’uso 

(3) L. P. HammeTT: Physical Organie Chemistry (New York, 1941), cap. VII; 
H. H. Jarré: Ohem. Rev., 53, 191 (1953). 

(4) G. W. WHELAND e L. PAULING: J.A.C.S8., 57, 2086 (1935). 
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corrente. Inoltre, conformemente a JAFFE, abbiamo posto %, = Yn a,, dove 
do, è l’integrale di Coulomb per l’atomo di C separato da n legami dall’atomo 
di N. Naturalmente la nostra equazione secolare è risultata del nono ordine 
e fattorizzabile in una equazione di sesto e una di terzo grado. Le radici di 
queste equazioni rappresentano i livelli energetici corrispondenti ai vari orbi- 
tali molecolari y; = > ay, (dove: j, î = 1, 2,..., 9; e le y, rappresentano gli or- 
bitali atomici p,). ‘ 

Abbiamo calcolato successivamente i valori dei coefficienti a;; e le densità 


elettroniche ai vari atomi: qg;= 2 Y aj, (j = 1, 2, ...,5 in corrispondenza ai 
5 orbitali leganti). Con i valori q, e q, per varie coppie xx; y, abbiamo costruito 
i grafici riportati in fig. 1, 2, 3, dove è € = (1— q,)/(1— qm) = 0/0». Per il 
nitrobenzolo = 1,1. Pertanto, trasportando la linea £ = 1,1 dal grafico 3 


=1.01 


In 


-05 0:5 10 15 -05 05 1.0 15 -05 05 Non ane 
Fig. 1. — Linee di livello Fig. 2. — Linee di livello Fig. 3 — Linee di livello 
Gm = fil%xs Yxc): G = fe(%x: Yxc): E = fs(ax: Yxo)- 


sul grafico 1 o sul grafico 2, l’intersezione con la linea q,, = 0,954 o con la 
linea g, = 0,950 dovrebbe dare i valori di x, e y,,. 

Ma nel campo dei valori ragionevoli di x, e yy, non si ha intersezione; 
inoltre la linea € = 1,1 cade nel campo dei valori negativi di «,, mentre 
l’elettronegatività dell’azote relativa al carbonio, aumentata poi dalla pre- 
senza dell’ossigeno, e il fatto dimostrato sperimentalmente (5) che il momento 
elettrico è diretto dal nucleo benzenico verso il gruppo NO, permettono di 
prevedere un valore positivo per x; in corrispondenza poi alla linea £ = 1,1 
è dn > Am > 1, in disaccrodo con l’esperienza. 

Benchè apparisse poco probabile che l’insuccesso del metodo di Jaffé di- 
pendesse dalla scelta fatta per x, e y,,, abbiamo saggiato l’influenza di tali 


(5) L. G. Groves e S. SUGDEN: Journ. Chem. Soc., 1094 (1934). 
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parametri per valori di x, e y,, Compresi tra 0,5 e 2: in nessun caso abbiamo 
ottenuto risultati anche approssimativamente soddisfacenti. 


2. — Momento elettrico di risonanza e transizioni elettroniche N — y. 


Abbandonate pertanto le costanti o come grandezze sperimentali fonda- 
mentali per la scelta dei parametri nel caso del nitrobenzolo, abbiamo rivolto 
la nostra attenzione ad un’altra costante fisica direttamente collegata con le 
concentrazioni degli elettroni x: il momento elettrico di risonanza. Tale mo- 
mento, dovuto alla non uniforme distribuzione degli elettroni x, in una mole- 
cola aromaticà monosostituita C,H,X ha un valore determinabile come dif- 
ferenza tra il momento di tale molecola e quello di una molecola RX di rife- 
rimento in cui R sia un radicale alchilico che non contenga elettroni x. Per 
il nitrobenzolo abbiamo scelto come sostanza di riferimento il nitrometano. 
I valori sperimentali più attendibili per i momenti dipolari di queste sostanze, 
espressi in unità Debye, sono: 


per: il: nitrobenzolo< 3°. =. 4,200) 04,22 (8) 


peril nitrometanoz.-. ee n= 9:01 (3), (0 IAA). 


Da questi dati abbiamo dedotto un valore medio per il momento di riso- 
nanza del nitrobenzolo: u= 0,75 D. 

Per il calcolo di tale momento per via teorica è necessario conoscere oltre 
i valori di qg per tutti gli atomi, la struttura geometrica delle due molecole. 

Per quanto sappiamo, non esistono dati sperimentali per il nitrobenzolo 
(liquido a temperatura ordinaria); pertanto dai lavori di LLEWELLYN (°) e 
ABRAHAMS (!°) relativi allo studio roentgenografico dei cristalli di p-dinitro- 
benzolo, abbiamo dedotto il seguente modello approssimato per la molecola 
CH;NO, allo stato di gas: nitrogruppo e nucleo benzenico coplanari (*), rego- 
larità dell'anello benzenico esagonale, N—O = 1,23 À, N-C=141À, 
ONO == 124°. 

Abbiamo inoltre ammesso che la struttura geometrica del gruppo NO, 
nel nitrometano sia praticamente uguale a quella dello stesso gruppo nel nitro- 


K. B. McALPINE e C. P. SmytTH: Journ. Chem. Phys., 3, 55 (1935). 
L. G. Groves e S. SUGDEN: Journ. Chem. Soc., 158 (1937). 

. P. SmrTH e K. B.McALPINE: J.A.C.S., 56. 1697 (1934). 

F. J. LLEWELLYN: Journ. Chem. Soc., 884 (1947). 

) S. C. ABRAHAMS: Acta Cryst., 3, 194 (1950). 
*) Nel calcolo di y l'eventuale angolo fra i piani del nucleo benzenico e del nitro- 
gruppo, imputabile alla repulsione di Van der Waals tra gli atomi di ossigeno e gli 
atomi di idrogeno in posizione orto, non ha influenza. 
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benzolo: i dati sperimentali trovati in letteratura (11?) sono sostanzialmente 
in accordo con questa ipotesi. Per il calcolo del momento dipolare di meso- 
meria sono necessarie anche le concentrazioni elettroniche per gli atomi di 
azoto e ossigeno nel nitrometano, calcolate per gli stessi valori di &,, yy € % 
usati per il nitrobenzolo. Tale calcolo è stato pure eseguito secondo il metodo 
standard M.O.-L.C.A.O., trascurando l’effetto di iperconiugazione. Il calcolo 
di u è stato eseguito per tutti quei casi che avevano presentato densità elet- 
troniche ragionevoli (in particolare q, > 4, > 4): i valori trovati sono com- 
presi tra 0,8948 e 4,9290, cioè sono tutti superiori al dato sperimentale: però 
il minore di essi, corrispondente al gruppo di parametri: 


(2) dora Pe ti oe AUDE Veo = 0,5 


può essere considerato, nell’ambito delle nostre approssimazioni, in discreto 
accordo con l’esperienza. D'altra parte i valori così ottenuti rappresentano 
una scelta ragionevole, perchè sono in accordo con le elettronegatività del- 
l’ossigeno e dell’azoto relative al carbonio (1%) e con i presunti ordini dei 
legami N—O e N—C. Anche i valori per le concentrazioni elettroniche si pos- 
sono considerare soddisfacenti. 


G = 09707, 9g, = 099818,  qn=0,9939. 


Per tutte le scelte di integrali di Coulomb e di scambio da noi eseguite 
abbiamo determinato mediante il metodo semiempirico di MATSEN (14) un’altra 
grandezza confrontabile con l’esperienza, precisamente la variazione di energia 
corrispondente alle transizioni N->V, cioè al passaggio di elettroni x su livelli 
eccitati: esse corrispondono al caratteristico assorbimento nel primo ultra- 
violetto, osservato negli spettri di molecole contenenti sistemi di doppi legami 
coniugati. Le determinazioni dello spettro di assorbimento nell’ultravioletto 
del nitrobenzolo in soluzione sono numerosissime, ma più interessante per noi 
è lo spettro allo stato di vapore. L’unico lavoro in materia pare essere la deter- 
minazione fatta da BERTON per scopi analitici (5), dove per il primo massimo 
del benzolo e del nitrobenzolo sono dati i valori: 

lari = 2 880 Ad 2 6610 À, 


Ce 


(11) L. O. Brockway: Rev. Mod. Phys., 8, 231 (1936). 
(12) F. RoGowSKI: Ber. dtsch. chem. Ges., B 75, 244 (1942). 
(3) L. PAULING: The Nature of Chemical Bond (Ithaca, 1939). 
(14) F. A. Marsen, W. W. RoBERTSON e J. F. Music: J.A.0.8., 72, 5243, 5248 
5250, 5252, 5256, 5260 (1950). 
(5) A. Berton: Ann. de Chim., (11) 19, 394 (1944). 
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da cui si ricava per il nitrobenzolo un valore «sperimentale » AH = — 1,8 
(nelle solite unità f). Tra tutti i valori teorici il più vicino è ancora quello 
ottenuto mediante gli integrali (2), eguale a — 1,7880. 


3. — Energia di risonanza e diagramma molecolare. 


Adottati pertanto i valori (2), abbiamo calcolato alcune grandezze più 
facilmente determinabili con il metodo M.O. e precisamente l’energia di riso- 
nanza, gli ordini di legame mobile e gli indici di valenza libera. 

La prima grandezza è ricavabile come differenza dell’energia degli elettroni z 
nel caso di completa delocalizzazione e nel caso di un sistema fittizio in cui 
tutti i legami siano localizzati: 


È risultato E.R. = 2,9146 B = — 58,292 kcal/mole, se si assegna a f il 
valore corrente di — 20 kcal/mole. 

Il confronto con l’energia di risonanza « empirica », dedotta dai valori del- 
l’energia dei legami presenti nella molecola e dal calore di combustione del 
nitrobenzolo, presenta alcune difficoltà. 

Una difficoltà di carattere generale in problemi di questo tipo deriva dal 
fatto che la maggior parte dei dati termochimici dai quali si deducono le 
energie di legame a 0 °K sono invece riferiti alla temperatura di 20 0 25 °C: 
peraltro le imprecisioni così introdotte sono trascurabili nel nostro ordine di 
approssimazione. Una seconda difficoltà deriva dalla nota discordanza dei 
valori dell’entalpia di dissociazione della molecola N, e di sublimazione del 
carbonio riportati in letteratura. La questione è tuttora aperta: peraltro, 
soprattutto in base alle esperienze di KISTIAKOWSKY, KNIGHT e MALIN (16) 
abbiamo adottato i valori seguenti: 


LN, >N, AH = 113.0 kcal/mole , 


293.1 °K 


C (grafite) + C (stato 3P) AH = 171,6 keal/mole . 


293.1 °K 


La difficoltà più grave deriva dalla mancanza di precise determinazioni 
(CO), KISTIAKOWSKY, H. T. Knrgut e M. E. Matin: Journ. Chem. Phys., 
20, 876 (1952). 
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A0) 
per il gruppo nd : l’unico valore trovato in letteratura (SYRKIN e DyAT- 
N-0 
KINA (17)) è più che altro un ordine di grandezza: 169-186 kcal/mole. L’osta- 
colo può essere evitato calcolando Venergia degli elettroni x per il seguente 
sistema fittizio 


nel quale sono localizzati solo i doppi legami dell’anello benzenico. Si ottiene 
così l’energia di risonanza del nitrobenzolo dovuta al solo anello benzenico: 


E.R. = 2,1048 = — 42,08 kcal/mole. 


Facciamo notare che la differenza col valore precedente (~ 16 kcal/mole) 
è circa eguale all’energia di risonanza calcolata per lo ione allilico : OH = CHSGHS 

Per ottenere il valore «sperimentale » da confrontare, occorrerà l’energia di 
legame empirica del nuovo sistema fittizio: a tale scopo è necessario conoscere 
(ad esempio dall’entalpia di formazione del nitrometano gassoso) l’energia di 
legame del gruppo NO,. Abbiamo adottato per le entalpie di formazione del 
nitrobenzolo e del nitrometano liquidi dagli elementi nello stato standard e 
per le relative entalpie di evaporazione i dati riportati da D’ANS e Lax (18), 
e per le energie di legame i valori di COATES e SUTTON (1%), senza applicare 
correzioni, nonostante alcune lievi differenze nelle temperature di riferimento 
(comportanti peraltro errori certamente trascurabili in rapporto alle nostre 
approssimazioni). Abbiamo ottenuto: 


E.R. (empirica) = — 41,41 keal/mole , 


in ottimo accordo con il valore teorico. 
L’ordine di legame mobile tra gli atomi » e À, pn;, e Vindice di valenza libera 


(7) Y. K. SyREIN e M. E. DyATKINA: Structure of Molecules (London, 1950), p. 240. 

(18) J. D'ANS e E. LAnx: Taschenbuch für Chemiker und Physiker (Berlin, 1949), 
p.721 

(19) G. E. Coates e L. E. Sutton: Journ. Chem. Soc., 1187 (1948). 

(70) C. A. Courson: Proc. Roy. Soc., A 269, 413 (1939). 


x 
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dell’atomo n, F,, nel nitrobenzolo sono definiti dalle relazioni: 


Pri = > 24 jn ji (20) ’ Fa ii 2 DE Pri (3) 1 
4) 


dove gli a;, e 4;; sono i coefficienti degli orbitali atomici dell'atomo » e à nel- 
l’orbitale molecolare j e la sommatoria va estesa ai cinque orbitali molecolari 


leganti; la >. pri va estesa a tutti i le- 
gami partenti dall’atomo n e come atomo n 
si intendono gli atomi di carbonio in posi- 
zione orto, meta o para al nitrogruppo. I 
risultati numerici sono riportati nel dia- 
gramma molecolare di fig. 4, unitamente ai 
valori delle concentrazioni elettroniche e ai 
valori delle distanze interatomiche prevedi- 
bili. Questi ultimi sono stati ricavati per 
i legami C—C mediante la formula (?°): 


Res at ag 
1 +k(1— p)/p 
dove R è la lunghezza cercata in A, p è Vor- 
dine di legame mobile, s è la lunghezza di 
un legame semplice (1,540 A), d quella di 
un legame doppio (1,340 A) e k è una co- 
stante (0,8095) scelta in modo che la curva 
«lunghezza di legame/ordine di legame » 


(o) O 15969 
SE 


0.8356 


ou 
oO 
(©) 
N 
Le} 
u 
a 
10501 
0.093 0.9810 
SOLI 
ox 
AS 
D 
n 
0.080 0.9939 
ia = —T_—_m 
0.087 0.9707 
Indici di valenza libera Distribuzione degli elettroni TT 
Fig. 4. - Diagramma molecolare 


per il nitrobenzolo. 


passi per il punto corrispondente al legame triplo (1,204 À). 

I valori delle costanti qui riportati non sono quelli originali di COULSON, 
ma quelli recentemente pubblicati da MOSER (??). 

Per il calcolo delle lunghezze dei legami C—N e N—O abbiamo invece 
fatto ricorso alla formula di Gordy (?*), malgrado questa sia stata stabilita 
per atomi con cariche relative trascurabili (noi abbiamo invece 1,1644 cariche 
elettroniche formali sull’azoto e 0,5969 su ciascun atomo di ossigeno). Peraltro 
i valori ottenuti risultano in buon accordo con i valori «sperimentali » citati 


precedentemente. 


(21) C. A. Coutson: Faraday Soc. Discussion, 2, 9 (1947); Journ. de Chim. Phys., 


45, 243 (1948). . 


(22) C. M. Moser: Journ. Chem. Soc., 1073 (1953). 
(23) W. Gorpy: Journ. Chem. Phys., 15, 305 (1947). 
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4. — Effetto orientante del nitrogruppo. 


Conoscendo le densità elettroniche, gli indici di valenza libera e i livelli 
energetici degli elettroni è possibile calcolare per il nitrobenzolo l’ordine di 
reattività delle posizioni orto, meta e para nelle reazioni di sostituzione nucleo- 
fila, elettrofila e di tipo radicalico. I metodi approssimati a disposizione sono 
essenzialmente. due (24): l’approssimazione della molecola isolata e la teoria 
della localizzazione, entrambi basati sull’ipotesi che nel confronto tra la reat- 
tività delle varie posizioni in una molecola contenente elettroni x l’entropia 
di attivazione si possa considerare costante, e l’energia di attivazione dipenda 
essenzialmente dalla variazione dell’energia degli elettroni 7 nel passaggio 
dalla molecola iniziale al complesso attivato. Per l'applicazione della teoria 
della localizzazione abbiamo determinato i valori dei livelli energetici nei se- 
guenti modelli di complesso attivato: 


ASL 

| | 
prata 
Seed SE 


I valori così ottenuti per l’energia di attivazione nelle varie posizioni e 
per i vari tipi di reazione sono riportati nella tabella I (in unità B). 


TABELLA I. 
| 
| Orto | Meta Para Reazioni 
| 
— 2,4506 | — 2,5186 — 2,3988 nucleofile 
— 2,5572 ANSE | SS KG radicaliche 
—- 2,6628 — 2,5968 | 2 0702 elettrofile 


Pertanto gli ordini di velocità risultano i seguenti: 
per la teoria della localizzazione, 


e) m>(p,0) 
r) (p,o)>m 
n) . (p,0) > m ; 


(4) R. D. BROWN: Quart. Rev., 6, 63 (1952). 
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per il metodo della molecola isolata, 


e) m> (0, p) 
T) (590) SS Go 
N\-(p. 0) > m. 


Per quanto riguarda la differenza tra la posizione meta e le posizioni orto 
e para troviamo completo accordo tra i due metodi e l’esperienza: la « chemical 
non-crossing rule » di BROWN (?4) è verificata. Il comportamento relativo delle 
due posizioni orto e para non è ben definito nè sperimentalmente nè teorica- 
mente: in questo caso però gli effetti sterici rendono non applicabile l’ipotesi 
della costanza dell’entropia di attivazione. 

Tra le convenzioni di calcolo introdotte da JAFFÉ appare facilmente criti- 
cabile (24) quella relativa all’assegnazione degli integrali di Coulomb agli atomi 
di carbonio in posizione orto, meta e para. Abbiamo pertanto ripetuto il cal- 
colo con l’usuale scelta dei parametri, ponendo «, = 0 per tutti gli atomi di 
carbonio eccetto quello al quale è attaccato il nitrogruppo, per il quale RELAIS 

Abbiamo ottenuto per le varie grandezze una concordanza con i dati speri- 
mentali superiore a quella riscontrata con il calcolo precedente. Ecco i dati 
che interessano: 


q, = 0,9640 Im = 1,0025 n,— 0,913 
fe = 051425) AE = —1,7916 E.R. = — 58,308 kcal/mole. 


TABELLA II. — Energie di attivazione 


Orto Meta | Para | Reazioni | 
bet: | | | 
— DS _ 2440? | — 2,4064 nucleofile 
— 2,5866 — 2,4402 — 2,6076 elettrofile 


Poichè questo trattamento semplificato pare possa dare risultati altret- 
tanto buoni, esso potrà essere utilmente applicato a molecole più complesse 
contenenti il nitrogruppo. 

Complessivamente il nostro metodo ha permesso di determinare una scelta 
di integrali in accordo con le previsioni basate sull’elettronegatività degli atomi 
di azoto e ossigeno e con le caratteristiche dei legami N—-O e N—C nel nitro- 
benzolo; l’impiego di questi parametri ha portato a valori in buon accordo con 
l’esperienza per tutte le grandezze prese in considerazione. 
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SUMMARY 


A molecular orbital treatment of the nitrobenzene molecule is given: dipole moment, 
U.V. absorption frequency, electronic densities, bond orders, free valences, resonance 
energy have been calculated with a reasonable choice of the involved parameters and 
with good agreement with experimental data. The kinetic behavior of the molecule 
in all types of substitutions is also correctly predicted: the chemical non-crossing rule 
is verified. 
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Sulla cattura dei mesoni »- da parte dei nuclei leggeri. 
A. ALBERIGI-QUARANTA e E. PANCINI 


Istituti di Fisica delle Università di Roma e Genova 


Istituto Nazionale di Fisica Nucleare - Sezione di Roma 


(ricevuto il 5 Aprile 1954) 


Riassunto. — Viene descritta una esperienza per decidere la validità della 
legge di Wheeler sulla dipendenza della probabilità di cattura dei mesoni u= 
da parte dei nuclei leggeri (Z < 26). A tale scopo è stata misurata, con 
un apparato descritto in un precedente lavoro, la vita media apparente 
dei mesoni uw in Al e S. Risulta in Al +t = 1,02 + 0,04 us e in $ 
t=0,70+0,04 us. Utilizzando questi dati, quelli ottenuti in un prece- 
dente lavoro per il C, ed un risultato ottenuto da KEUFFEL e coll. per 
il Fe, si conclude che la teoria di Wheeler è sufficiente a spiegare, nei limiti 
delle precisioni oggi raggiungibili, la dipendenza della probabilità di 
cattura p, dei mesoni u= da parte dei nuclei con Z < 26. Si ricavano i 
valori k = 4,4 + 0,3 e Z, = 11,3 per le costanti introdotte da WHEELER 
nella legge p, = (1/70) (Ze/Zo)*. 


Introduzione. 


La dipendenza dal numero atomico Z della probabilità di cattura nucleare, 
pe, di un mesone pu, che si trovi sull’orbita AK di un certo atomo sembrava 
aver trovato soddisfacente spiegazione nella teoria di WHEELER (!). 

Recentemente, però, gli interessanti risultati delle misure delle probabi- 
lità della cattura dei mesoni u- da parte di nuclei pesanti (Z > 30), eseguite 
da KEUFFEL et al. (2), mostrano tali discrepanze con le previsioni della teoria 
di Wheeler da indurre a considerare opportuno un ulteriore studio sperimen- 


(1) J. WHEELER: Rev. Mod. Phys., 21, 133 (1949). 
(2) J. W. Keurret, F. B. HARRISON, T. N. K. Goprrey e G. T. REYNOLDS: Phys. 
Rev., 87, 942 (1952). 
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tale del fenomeno in questione anche nel campo dei nuclei leggeri per i quali 
la validità della teoria è basata soltanto sui risultati delle esperienze di 
TICHO (8). 

Tali risultati infatti, a causa dei forti errori statistici dai quali sono affetti, 
non permettono nè di confermare il valore % = 4 nella formula 


a) ee Gi 


TONZE 


(Ove To = 2,22 us è la vita media naturale del mesone u, Z, una funzione di Z 
calcolata da WHEELER e Z, una costante), nè di mettere in luce eventuali 
deviazioni da tale legge per alcuni nuclei di particolare struttura. Ci è sem- 
brato perciò opportuno verificare più rigorosamente, di quanto non fosse stato 
fatto fino ad ora, la legge di Wheeler. 

A tale scopo abbiamo eseguito misure più accurate delle vite medie appa- 
renti t dei mesoni u negativi condotti a riposo in elementi di basso numero 
atomico (Carbonio, Alluminio e Zolfo). Infatti dalla conoscenza delle vite 
medie apparenti si può risalire subito alla probabilità di cattura per micro- 
secondo. 


1. — Apparecchio. 


L'apparecchio da noi usato in questa ricerca è già stato descritto in un 
nostro precedente lavoro (*) a cui faremo costantemente riferimento anche per 
le notazioni. Durante le misure per la determinazione della vita media appa- 
rente del y- in Al e in S sono stati modificati, rispetto alla precedente ricerca, 
solamente i limiti dei gruppi di ritardi in cui vengono catalogati tutti gli 
eventi [(I, II, R) — AC]. 

Nelle misure di t in Al e Fe sono stati scelti i limiti: 


Gruppo da : a 
I — 0,075 us + 0,025 us 
II + 0,025 + 0,175 
III + 0,175 + 0,35 
IV + 0,35 i 0,65 
V + 0,65 + 1,00 
VI + 1,00 + 1,45 
VII + 1,45 + 2,25 
VIII + 2,20 E 10 


(3) H. K. TicHo: Phys. Rev., 74, 1337 (1948). 
(4) A. ALBERIGI-QUARANTA e E. PANCINI: Nuovo Cimento, 9, 959 (1952). 


SULLA CATTURA DEI MESONI uw DA PARTE DEI NUCLEI LEGGERI 609 


Nelle misure di 77 nelle S lo zero dei tempi è stato anticipato di 5-10-8 s, 
lasciando immutata l’ampiezza temporale di ciascun gruppo. 

Ricordiamo che l’istante dal quale si misurano i ritardi è determinato dal- 
l’arrivo del segnale dalla doppia D, e quindi dai contatori del gruppo II. 
Questo segnale segue il passaggio della particella attraverso i contatori II di 
un tempo che include un ritardo costante e un ritardo variabile statisticamente 
da 0 a 0,18 us (*). Da ciò segue che eventi fisicamente contemporanei in I, II, R 
possono essere registrati come eventi anticipati o ritardati a seconda del valore 
dei ritardi propri dei contatori II e R. 

Con la scelta dei limiti dei gruppi di ritardi esposta nella tabella sopra 
riportata si può calcolare (5) che il rapporto tra i numeri di eventi nel I e nel 
II gruppo deve essere circa 0,9 nel caso delle misure con Al e 1,75 nel caso 
delle misure con S. La costanza di questi valori nelle singole misure è stata 
usata come una delle indicazioni del buon funzionamento dell’apparecchio. 
Altre indicazioni sul regolare funzionamento dell’apparecchio sono tratte dai 
dati letti sui numeratori N;, N,, N,. Ogni misura parziale in cui uno qual- 
siasi di questi numeri differisse dal valore medio per più di quanto compa- 
tibile con le fluttuazioni statistiche è stata eliminata. 


2. — Risultati sperimentali. 


Con il dispositivo sperimentale già descritto sono state eseguite principal- 
mente misure per la determinazione della vita dei mesoni u- nell’ Al e nello S. 
Tali misure venivano fatte selezionando con lenti magnetiche mesoni uw e 
usando come assorbitori Alluminio metallico (densità = 2,7 g/cm*) o Zolfo 
(densità = 2 g/cm?) contenuto in una scatola di ottone con pareti di 0,5 mm. 
A queste venivano regolarmente intercalate misure senza assorbitore. 

Sono state eseguite per controllo anche serie di misure usando come assor- 
bitori Ferro in trucioli (densità = 1,2 g/em*) e serie di misure selezionando 
mesoni wr. 


a) Misure senza assorbitore. — I dati delle misure senza assorbitore sono 
riassunti nella tabella I. 

Nel primo gruppo vengono riportati i dati, le correzioni per le casuali e i 
dati corretti relativi alle misure eseguite senza assorbitore, selezionando x 
e con gli stessi limiti per i gruppi di ritardi usati per Ie misure con Al e Fe. 

Analogamente nel II gruppo vengono riportati i dati, le correzioni per 
casuali e i dati corretti, relativi alle misure eseguite senza assorbitore, selezio- 
nando u-, e con gli stessi limiti per i gruppi di ritardo usati per le misure con S. 

(CRAG ALBERIGI-QUARANTA, L. MezzeTTI, E. PANCINI e G. StoPPINI: Nuovo Ci- 
mento, 8, 618 (1951). 
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TABELLA I, — | 


: Gruppi 
Durata della misura Segno | di ritardi N pin Nam Nan 
come con 

ie LLY ett) ee ey ares = Al 56,7 + 0,6 156,7 +1 1084 4 
sf Correzioni casuali . . . . AI — — = 
È Valori COrrevtiE- RE _ AI — — = 
= LOS TRINO = S 46,7 + 0,7 159 + !,4 1 069 + 
Fe Correzioni casuali . . . . — S = = = 
Valonecorroutie nee an S — -- — 

DOCS ESS St media valori corretti — 
TABELLA II. — Mis 

Durata della misura Segno Nom | Nei 

IO nnt aes 62,4 + 0,8 163 + 

Correzione per le coincidenze casuali . . = sn 
Correzionemper tl) fondom RCE = es 

| Wotalestorrezioni ei utero) == = 

At COLL LUI Meee iru a eet Re ae, = LT 

Dati sprevwistiapet cee (020 = = 

SLOPE RE pee ae oes. Sree + 71,5 + 1,5 158 + 

Correzione per le coincidenze casuali . . = | cd 

Correzione per i fondo to = =a es 

Mo Caley COLreZlONI's a. ame ere ae = n 

Daticcorretty n RI ee P| = 32 
TABELLA III. — Misure con asse 

Materiale | Durata della misura | Segno Nom Nan Nam 

Al | 452.46 h | = 98,9 + 0,5 163,9 + 0,6 | 1108,3 + 

SA AI 75.75 h lee Sei 163,5 + 1,5 1160,7 + 
i ki Ba n 
‘ S 514.95 h = 78,1 + 0,4 170,0 + 0,6 1077,2 + 
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PI 
assorbitori. 
f Ny Nm My Ny Nyy Nyy Nvm 
153 1336 105 11 11 9 15 47 
= cs — 1 1 122 1,8 18 
153 1336 105 10 10 13,2 29 
461 284 36 7 1 —= 3 7 
= = — _ — — 0,6 6 
461 284 36 7 1 — 2,4 1 
— ca 35,2 4,3 3,3 2 3,9 159 
lo con assorbitore di Fe. 
Nam N; Ny Nin Mwy Ny Nyy Nyy Nvm | 
51 +.3,4 399 654 78 253 11 9 5 29 
= = == — = — — | 1 9 
a — — 33 4 3) 2 4 7,9 
= — — 33 4 3 2 5 16,5 
. — — 35 19 8 7 = 12,5 
— ae =a 35 24 4 0.8 nu = 
99 + 6 155 215 39 18 12 25 30 57 
— a — — — — — a 3 
E Ets = 9 2 4 4 7 12 
<# L= e 9 2 4 4 7 15 
e a = 30 16 8 2A 23 42 
e S. Dati corretti. 
Nu Nu Ny Ny Nyy 3. Nvyn Nein 
a 3 438 185 + 19 197 +15,4 220 + 16 191 + 15 133 + 13,5 155 + 17 
2 480 47 51 DS 89 WS 234 
2 2 066 116 +17 143 + 14 121 + 12 75 + 10 9 + 7,4 80 + 14,5 
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Il cambiamento dei valori di N,/h e N;/h nel passare dal primo al secondo 
gruppo di misure, è dovuto al fatto che la pressione atmosferica media du- 
rante il periodo in cui vennero eseguite le misure appartenenti al primo gruppo 
era alquanto minore (circa di 1 cm di Hg) di quella che si è avuta durante 
le misure del secondo gruppo. 

Come si vede chiaramente l’avere anticipato di 5-10-* s lo zero temporale, 
passando dal primo al secondo gruppo di misure, altera sostanzialmente sol- 
tanto le distribuzioni relative di eventi tra il I e il II gruppo di ritardi, mentre 
la frequenza di eventi nel III gruppo resta immutata nei due casi, come era 
da attendersi. 

Quest’ultima osservazione dimostra che nel caso dell’ Al agli eventi nel 
III canale non contribuiscono affatto i ritardi propri dei contatori del gruppo R. 
Infatti, se così fosse, lo spostamento dello zero dei tempi, e quindi del centro 
della curva di dispersione dei ritardi propri dei contatori, porterebbe, nel pas- 
saggio da un gruppo di misure all’altro, ad una notevole diminuzione della 
frequenza degli eventi nel III canale; il che non accade affatto. Accertato 
quindi che al III canale non contribuiscono i ritardi propri dei contatori (e se 
ciò è vero per le misure in Al lo sarà a maggior ragione per le misure in S) si 
possono senz’altro comprendere, nel caso delle misure con assorbitori, i dati 
relativi a questo canale nel calcolo delle vite medie. 

Nelle misure senza assorbitori gli eventi relativi ai canali dal III all’ VIII 
sono molto probabilmente dovuti a mesoni negativi (e a una piccola percen- 
tuale di positivi) disintegrantisi nei materiali di vario numero atomico costi- 
tuenti il sostegno e le pareti dei contatori con una vita media apparente abba- 
stanza breve (~ 1 us). 

Data la differenza trascurabile che il piccolo spostamento dello zero del 
tempo può portare sui risultati nel passare dalla prima alla seconda serie di 
misure, essi vengono senz’altro sommati considerandoli come facenti parte di 
uno stesso gruppo di misure. 


b) Misure con assorbitori di Fe. — Nella tabella II sono elencati i risul- 
tati di una serie di misure di controllo ottenuti usando come assorbitori due 
scatole contenenti Ferro in trucioli (densità = 1,2 g/cm?) di sezione trapezoi- 
dale e di cm 12 di spessore pari a 14,4 g/cm?. 

Nelle misure con assorbitore di Fe e mesoni u-, il numero di N,/h cresce, 
rispetto alle misure senza assorbitore, di 17 + 0,3/h; questo aumento cor- 
risponde ad un coefficiente medio di assorbimento, per i mesoni per questa 
regione di energia (intorno a 500 MeV), uguale a (0,12 + 0,03)% per g/cm? 
di Fe, in buon accordo con gli analoghi risultati da noi trovati in C, Al e S e 
con quanto si può dedurre dallo spettro di energia dei mesoni a 4300 m s.l.m. 
trovato da D. B. HALL (°). 


(5) D. B. HALL: Phys. Rev., 66, 321 (1944). 
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I corrispondenti dati per i mesoni positivi non sono stati calcolati per la 
povertà della statistica. 

Nella tabella II vengono elencate anche le correzioni effettuate sui dati 
sperimentali. Tali correzioni vengono calcolate per gli eventi registrati nel 
gruppo di ritardi dal terzo all’ottavo, poichè soltanto questi gruppi interes- 
sano per la determinazione della vita media 7. 

Le correzioni effettuate sono: 


1) Correzione per coincidenze casuali: da ogni gruppo di ritardi 
viene detratto il numero di coincidenze casuali calcolate. 


2) Correzione per il fondo: da ogni gruppo di ritardi viene detratto il 
numero di eventi registrati nello stesso gruppo durante le misure senza assor- 
bitore, ragguagliato allo stesso numero di ore. In questo modo si pensa dî 
eliminare dai dati il fondo di disintegrazioni avvenute nelle pareti e nei sostegni 
dei contatori, fondo che è ragionevole ritenere indipendente, in prima appros- 
simazione, dalla presenza dell’assorbitore. 

Nella tabella II è anche riportata la distribuzione dei ritardi calcolata per 
una vita media apparente nel Fe di 0,2 us, quale può essere prevista dalla 
teoria di Wheeler e in buon accordo con le misure di KEUFFEL et al. (?); la 
normalizzazione è stata eseguita in modo da avere nel terzo gruppo un numero 
di eventi pari a quello sperimentale. Come si vede non vi è accordo tra il 
risultato previsto e quello sperimentale almeno per quanto riguarda i gruppi 
di ritardi VI, VII, VIII, e questa discrepanza va ovviamente attribuita alla 
incompleta efficienza del dispositivo nel selezionare mesoni di segno diverso. 

In realtà l’efficienza delle lenti magnetiche può esser valutata con il cal- 
colo (7) e risulta praticamente totale, mentre i dati sperimentali da noi otte- 
nuti e riportati in tabella II sembrano mostrare che un numero di mesoni 
positivi, pari circa al 5% di quelli che si disintegrano nell’assorbitore con 
le lenti disposte per convergere mesoni positivi, si disintegrano nell’assorbitore 
anche quando le lenti sono disposte a selezionare mesoni negativi. 

Questo fatto riscontrato anche da H. K. TICHO (*) con un dispositivo ana- 
logo al nostro non è spiegabile se non ammettendo che tali mesoni positivi si 
arrestino nell’assorbitore senza avere attraversato entrambi i contatori del 
telescopio (I, II), e siano però associati a particelle che azionano almeno uno 
dei contatori I e II. In tale modo essi non risentono dell’azione delle lenti 
magnetiche, ma possono dare luogo ugualmente ad eventi [(I, II, R) — AC] 
che vengono quindi esaminati dal misuratore di ritardi 7. 

Questa ipotesi è anche confermata dal fatto che H. K. TIicHo trova per 
la frazione dei mesoni positivi incidenti non eliminata dalle lenti magnetiche, 
un valore del 4%, assai vicino a quello da noi trovato, mentre la geometria 


(7) M. Conversi, E. PancINI e 0. Piccioni: Nuovo Cimento, 3, 372 (1946). 
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del suo dispositivo, alquanto diversa dalla nostra (uguale lente magnetica, 
uguale campo, ma angolo solido del telescopio assai più stretto del nostro), 
indurrebbe ad attendere una disparità assai maggiore. 


3) Misure con assorbitore di Al. — Nella tabella III sono elencati i risul- 
tati ottenuti usando come assorbitori due blocchi di alluminio (densità 2,7 g/cm?) 
di sezione trapezoidale e di em 12 di spessore, equivalente a 32 g/em?. Nelle 
misure con assorbitore il numero di N,/h cresce, rispetto alle misure senza 
assorbitore, di (47,5+ 1,7)/h per i mesoni positivi e di (42,2+ 0,8)/h per i me- 
soni negativi. Questo corrisponde a un coefficiente medio di assorbimento per 
i mesoni in Al uguale a (0,16 + 0,01)% per g/em? di Al per i positivi e a 
(0,12 % 0,002)% per g/cm? per i negativi. 

Notiamo inoltre che gli incrementi di N,/h che si riscontrano passando 
dalle misure senza assorbitori a quelle con assorbitori di Al, e che sono dovuti 
ai mesoni incidenti nei telescopi (I, II) che vengono condotti a riposo negli 
assorbitori, presentano un eccesso positivo (*) del (11,8 + 4,2)%. Conside- 
rando invece gli eventi registrati nei gruppi di ritardi dal terzo all’ottavo 
dovuti a mesoni disintegrati nell’assorbitore, si può valutare (tenuto conto 
della probabilità di cattura dei u- e della diversa vita media dei + e dei 4°) 
che per i mesoni condotti a riposo negli assorbitori si ha un eccesso positivo del 
(20,7 + 5,8)%, non incompatibile con il valore trovato con le N,/h. 

I dati della tabella III sono dedotti dai risultati sperimentali effettuando 
le stesse correzioni che nel caso del Fe ed inoltre è stato tenuto conto, per i 
mesoni negativi, anche della correzione per l’incompleta selezione operata dal 
dispositivo. Abbiamo infatti calcolato, assumendo 7* = 2,22 us, il numero di 
eventi per ciascun gruppo dovuti a disintegrazione di questi mesoni positivi 
(5% del numero di u+ incidenti nell’angolo solido del telescopio) sfuggiti alla 
selezione delle lenti, e lo abbiamo detratto dal numero totale di eventi regi- 
strati in quel canale. 

Per il calcolo della vita media, non essendo applicabile nel nostro caso il 
metodo di Peierls (*) abbiamo usato le ben note formule 


5 irc VA Ci 
n dg CT Cle A 


dove r è il rapporto fra il numero totale N di eventi ritardati dopo il tempo ¢, 
e il numero N, di eventi ritardati dopo il tempo t. 


(*) Nel presente lavoro si intende con l’espressione eccesso positivo, la quantita 


Nee 
2 5 


Sa HOOF. 
Nea N= 


(8) R. PrrerRts: Proc. Roy. Soc., A 149, 467 (1937). 
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Dai nostri dati è possibile ottenere cinque valori indipendenti della vita 
media apparente del mesone u-: 


0,93 £ 0,26, 1,09+0,14, 1,01+0,08, 0,96+0,06, 1,07 +0,06 
da cui si ritrova per la vita media il valore 
Ty = 1,02 + 0,04 us. 


4) Misure con assorbitori di S. — Nella tabella III sono elencati i risul- 
tati ottenuti usando come assorbitori due blocchi di Zolfo (densità 2 g/cm8) 
di sezione trapezoidale e di cm 12 di spessore equivalente a 24 g/em?. 

Nelle misure con assorbitore il numero di N,/h cresce, rispetto alle misure 
senza assorbitore, di (31,4 + 0,8)/h per i mesoni negativi. Questo corrisponde 
ad un coefficiente medio di assorbimento per i mesoni in S uguale a 
(0,12 + 0,004)% per g/cm? di S. 

Usando le formule (2) si ottengono per la vita media i cinque valori indi- 
pendenti: 


0,73 0,25, 0,73 


Fit 222-0,69- 10,07, 0,70 = 0,07. 1,082 01 


come si vede l’ultimo di questi valori sembra scartare in maniera normale 
rispetto alla media pesata di essi che risulta 0,76 + 0,04. 

Dall’esame dei dati sperimentali e il confronto con le misure senza assor- 
bitori che sono state regolarmente alternate a quelle con assorbitori di S si 
può concludere che: 


a) gli eventi registrati nell’ottavo gruppo durante le misure con assor- 
bitore di S sono regolarmente distribuiti nel tempo; 


b) le misure senza assorbitore danno per l’ottavo gruppo il numero di 
eventi prevedibile in base a tutte le altre misure. 


Ciò indurrebbe a concludere che il numero anormale di eventi registrati, 
nel caso dello S, nell’ottavo gruppo non sia imputabile ad un effetto stru- | 
mentale. 

Le interpretazioni di questo effetto, se reale, possono essere diverse. Nel- 
l’attesa di indagare più a fondo il fenomeno dei mesoni y in $ con un nuovo 
dispositivo sperimentale, ci limiteremo ad assumere come valore più atten- 
dibile per 7; in S 


Ty = 0,70 + 0,04 ps, 
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che si può ottenere facendo la media pesata dei primi quattro valori sopra 
elencati. 


3.-— Conclusioni. 


Nella seguente tabella sono riportati i valori delle vite medie apparenti 77 
e della probabilità di cattura per microsecondo p, quali si possono desumere, 
dai risultati delle nostre misure. 

Insieme ai risultati da noi ottenuti per © (4), Al, S, sono stati elencati 
quelli relativi al Fe ottenuti da J. W. KEUFFEL et al (?). 


| | 
| Autori | Elemento | Vis ia De 
Ra 
| 

| NOR INO ai NATI ONE TS name. 18207 0,01 +0,01, | 
ME ree) C | 11,58 | 1,02 + 0,04 0,53 + 0,04 « | 
‘1 7 E II S je EB Te 070 50,04 0,97 + 0,06 

INEUPPED «220.04 Fe | 19,4 [Oo GIGS SE ULV BE |) Bites SE gay 


L’esame di questa tabella mostra anche con determinazioni alquanto accu- 
rate (almeno per misure eseguite con le piccole intensita dell’esperienze sui 
mesoni u della radiazione cosmica) dei valori di t i valori della p, siano affetti 
da errori che per i piccoli numeri atomici raggiungono il 100%. Si deve inoltre 
considerare che per la forma della dipendenza di % e Z, da p., Ventita degli 
errori relativi a queste ultime grandezze diventa in questo caso tale da non 
rendere più sicuramente applicabili quei risultati della teoria degli errori che 
sono validi quando gli errori sono piccoli rispetto alle grandezze cui si rife- 
riscono. 

Ciò nonostante abbiamo creduto opportuno calcolare i valori di % e Z 
applicando il metodo dei minimi quadrati non dimenticando per altro che 
i valori così ottenuti possone- non essere la migliore determinazione di tali 
grandezze. 

Utilizzando solo i nostri dati abbiamo ottenuto 


k=38+05, Z—111 (lg Z — 1,045+ 0,2) 


e tenendo conto insieme ai nostri risultati di quello di KEUFFEL relativo al 
Fe, abbiamo ottenuto 


k= 44403, Z—113 (lg Z =1,054+ 0,2). 


Si noti che la forte variazione tra i due valori di k è in gran parte dovuta 
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alla piccola variazione del valore di Z,. Infatti si ha dalla (1) 


k? 
dk = Beds 
Zo log RUE 
ni 


dalla quale si vede che se Z, ~ Z,, cioè per (to — t-)/t- prossimo ad 1 (che 
è il caso della nostra misura dotata di maggior peso) il coefficiente di dz, di- 
venta maggiore di 10. 

Date le grandi difficoltà di ottenere con i mesoni della radiazione cosmica 
dei valori della probabilità di cattura da parte di nuclei leggeri molto più 
precisi di quelli da noi riportati è prevedibile che la legge di dipendenza di 
queste probabilità da Z e la correttezza dei valori di Z, calcolati da WHEELER 
non potrà essere controllata con precisione molto maggiore di quella ottenuta 
nella presente ricerca almeno finchè non potranno essere utilizzati fasci 
di mesoni yw. artificiali. 

Si può dunque concludere che la teoria di Wheeler della cattura dei me- 
soni u- da parte dei nuclei leggeri è sufficiente per spiegare, nel limite delle 
precisioni oggi raggiungibili nelle esperienze sui raggi cosmici, la dipendenza 
delia vita media apparente dei mesoni t dal numero atomico, assumendo per 
le costanti Z, e ki valori 11 e 4 per i valori di Z non superiori a Z = 26. 


SUMMARY 


An experiment is described to test the validity of Wheeler’s law for the Z-dependence 
of the probability of capture of negative u-mesons by light nuclei (Z <26). With 
an apparatus described in a preceeding paper the apparent mean life of the y-mesons 
t- has been measured in Al and S, with the result 77 = 1.02 + 0,04 in Al and 
7= =0.70 + 0.04 in S. Using these values, those obtained previously for C, and 
a value obtained by Krurret et al. for Fe, it is concluded that Wheeler’s theory is 
sufficient to explain, within the limits of accuracy at present attainable, the depen- 
dence on Z of the capture probability p, in light nuclei (Z < 26). The values 
k = 4.4 + 0.3 and Z, = 11.3 are found for the constants introduced by WHEELER 


in the law p, = (1/79) (Z2/Z)F. 


en AR PT È 
j 


12 


Li 
pa 


ART 
SITO 


LEA 


+ 


Po EN RAR EURE) 


4 


L thy 
J 


IL NUOVO CIMENTO Vor Xr No 6 1° Giugno 1954 


High Polarization of Nuclei in Paramagnetic Substances. 


P. BROVETTO and G. CINI 


Istituto Nazionale di Fisica Nucleare - Sezione di Torino 


(ricevuto il 6 Aprile 1954) 


Summary. — In this paper a thermodynamical treatment of nuclear 
polarization in paramagnetic solids is developed. The general validity 
of the results obtained by A. W. OveRHAUSER for alkali metals, is con- 
firmed. We find that it is possible to achieve a degree of nuclear polarizat- 
ion equal to the one obtainable with the same magnetic field and tem- 
perature for electrons. 


1. — Introduction. 


A new method for producing polarization of nuclear spins in a metallic 
solid has been proposed recently by A. W. OVERHAUSER (1). 

This method is based on the effect produced by paramagnetic relaxation 
processes, due to coupling of electronic and nuclear spins (?): more exactly 
one considers the nuclear relaxation process in a metal placed in a static mag- 
netic field, which tends to restore the system of electronic spins, excited by 
the photons of a strong microwave source, to its equilibrium state, while 
inducing nuclear polarization. This polarization is calculated dynamically 
considering the evolution of the system towards the equilibrium state, viz. 
employing the knowledge of relaxation times. However, the polarization in 
the equilibrium state is found to depend only upon the values of the static 
magnetic field and of the temperature, and upon the degree to which the electron 
resonance can be saturated, hence not upon the nuclear relaxation time. 

In this work a very simple thermodynamical treatment of the problem is 


(*) A. W. OVERHAUSER: Phys. Rev., 92, 411 (1953). 
(2) A. W. OVERHAUSER: Phys. Rev., 89, 689 (1953). 
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developed, which enables us to obtain directly the results of A. W. OvER- 
HAUSER, extending their validity. 


2. — Thermodynamical treatment of the nuclear polarization. 


A number of relaxation processes brings to the equilibrium condition the 
difference of the population between the Zeeman levels of the electrons. The 
relaxation process due to the interaction between electrons and nuclei, depends 
upon the operator (*): 


87 
ae ee =: BeBn( 209402 = 20,-0+ 37 Spies) d(7) ? 


where f, is the electronic magnetic moment, f, is the ratio between the nu- 
clear magnetic moment and twice the nuclear spin I, 0,, 0, are the Pauli 
spin operators for the electron and for the nucleus, and 6(r) is the Dirac function 
of the distance between electron and nucleus. The first term of such an ope- 
rator is responsible for transitions of the type: 


(1) C,34 Mm: > Cy + Nm 
and the second term for the transition: 
(2) CES E Um > C14 37 Mm 3 


while the third term does not produce similar effects. 

We have denoted in (1) and (2), by e,, and e_, the electrons with com- 
ponents of the spin momentum respectively parallel and antiparallel to the 
static magnetic field, and analogously by 1, the nucleus with component m 
of the spin momentum J in the direction of the field. 

Let us call: N,,, N_, and N respectively the number of electrons ¢,,, ¢_, 
and the total number of electrons, M,, and M respectively the number of 
nuclei », and the total number of nuclei, that is to say: 


(3) Nee N ING; 


(4) a NM, =M. 


When the equilibrium between processes (1) and (2) is achieved, from the 


(3) E. FERMI: Zeits. f. Phys., 60, 320 (1930). 


620 P. BROVETTO and G. CINI 


law of mass action we have (*): 


NE VM, 
5 Bands eat ge 
) Nes 


where K, is obtainable by integration of the Van’t Hoff equation: 


dlog K, AU 
TONE 


(6) 


AU is the heat-content at constant volume, that is the variation of the in- 
ternal energy of the system, during the process (1). Therefore we have: 


(7) [e,1] ne [Ama] sé [e_,| as [nn] edili 5) 


where the square brackets mean the internal energy of electrons and nuclei. 
. From (7) it follows, remembering that f,, 6, have opposite signs, that: ’ 


AU=— 2H(6,,4+ |B.|)%, 


where H is the static magnetic field, 97 Avogadro’s number. 
Introducing the gyromagnetic ratio y = 26/h, the variation of the internal 
energy becomes: 


We determine the constant of integration in (6) by means of the conditions 
that at very high temperature we must have: 


Mo è VEST . 


From (5), (6) and (8) we then find: 


(9) NM, = exp pa nH A 


Nie ae kT 


(*) It is possible to obtain in statistical way the explicit expression of the mass 
action law (5) by assuming that the sums of the potentials of nuclei and electrons 
appearing in both sides of equation (1), expressed by means of their partition function, 
are equal. We have, in the Boltzmann quantum statistics, u=kT log (N/f(T)) where pu 
is the chemical potential and f(7) the partition function. However we feel that the 
thermodinamical considerations introduced here to evaluate HK, have a more direct 
physical interpretation. 
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Without microwaves, processes (1) and (2) cannot change the internal 
energy of the spin systems of electrons and nuclei, and so cannot modify the 
ratios of populations of electronic and nuclear spin levels in any way. These 
populations are due to the two independent conditions of the equilibrium: 


ty 63, 5 ips Sr o VIZ: 
NS [Vel M} y 
10 = ex 27) DR EE es ER À : 
(10) M, exp | ET h | Dan 06 exp FE nH 


If the system is strongly irradiated by microwaves with frequency 
o = 2|6.|H/h, the transitions e_ , + hw— e,, will increase the population of the 
level --1/2 and consequently a great number of transitions (1) occurs, inducing 
nuclear polarization. In other words, having reduced by means of the micro- 
waves the ratio N_,/N,,, when the equilibrium state is restored, the ratio 
M,,/ M, is found to be increased, as one can see by equation (9). 

Let us define: 


= 


11 pats ee ] 
( ) Ni i: exp [3] ’ 


where N_,, N; are the populations of the electronic spin levels when the altern- 
ating field of microwaves is turned on. At the saturation of paramagnetic 
resonance we shall have N_j= N,; and 6=0. 

Let us define, following OVERHAUSER a parameter of saturation s: 


ME SIE 
Nata 


(12) sil 


which is 0 when there is no alternating field and is 1 for complete saturation. 
The relation between s and Ô follows from (3), (10), (11), (12): i 


o == = ar, 
tgh (1 — ste OMT 


As the most interesting feature of the phenomenon arises when s ap- 
proaches 1, we make the following approximation: 


PAL: 
8) prec 


(13) 6=(1 
With the use of (13) and (11), (9) becomes: 


a ya tslvel à 
a = exp ZARE uz LI hH : 


ti 


(14) 
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The value of the ratio M,,/M,-, when s= 1, given by (14), is obtained as 
an exact result from (9). 
Let us define the nuclear polarization P in this way: 


TER M 
TM 2e M: 


m 
I 


P= 


By means of (4) and (14) the polarization becomes: 


I 
Da mexp [mu] 
(15) IP ze 


and therefore is given by the Brillouin function: 


21 +1, top et ali 
2I ae tae I, 


(16) Be cotgh 5 


This expression for P means that the polarization of nuclei, when the satu- 
ration is s, for given values of static field and temperature, is the same as the | 
polarization of a system of electrons would be, under the same conditions, 
if their gyromagnetic ratio were y, + s|y.|. 

For small values of w, viz. high temperature and low static field, (16) be- 
comes: 


Va + 8 Ye 


Beales) 


3kT hl 
For I — 1/2 (16) can be written: 
Vars || 
17 P = tgh ——_ AH. 
( ) 15 2kT i 


The equations (16) and (17) express in a more complete form the same 
results obtained by OVERHAUSER (1), as our formula (15) is formula (15) of 
his work. 

In Fig. 1 is plotted the value of P versus 7, calculated for a nucleus of 
spin 1/2 by means of (17) with s = 0.9, for different values of H and there- 
fore for different values of the wave lenght of microwaves (*). 


(*) The curves of Fig. 1 are valid for any nucleus with spin 1/2, because in every 
case Yn & Ye: 
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3. — Conclusion. 


The curves of Fig. 1 show that it is possible to obtain nearly total nuclear 
polarization at temperatures below 1 °K with centimeter waves, when the 
parameter s approaches unity.: 


= 21400 Gauss 
= 10700 Gouss 


= 3570 Gauss 
= 1070 Gouss 


n = n == == 
0 os i) 15 2 25 3 


Fig. 1. - Value of P calculated for a nucleus of spin 1/2. 


Having evaluated without paramagnetic saturation the. polarization of 
protons in a field of 21400 Gauss, at 0.5 °K, we find that under the same 
conditions of magnetic field and temperature, the polarization may be increased 
in this case by a factor of 230. 

Generally at temperatures of the order of 1 °K electronic relaxation times 
are quite long so that saturation may be achieved. 

For instance in the alkali metals saturation must be obtainable as the 
calculation of A. W. OVERHAUSER (?) predict values of 10-$ for the relaxation 
times already at room temperature. In metallic solids in order to avoid eddy 
current losses, it would be necessary to use finely divided metallic samples of 
particle size somewhat smaller than the depth of the skin-effect at low tem- 
perature, viz. smaller than a micron (*). 

We can conclude that the experimental conditions, under which the total 
nuclear polarization is predicted by (16), must be realizable at least in some 
cases, because they are the same conditions, apart from the paramagnetic 
saturation, encountered in the classical experiments of electron polarization (5). 

However even at room temperature and in a static field of 30 gauss the 
effect of nuclear polarization by means of this method, has been qualitatively 
checked ($). 

The effect has been shown by nuclear resonance, the intensity of which 
is proportional to the degree of nuclear alignment. 


(4) T. W. GriswoLp, A. F. Kip and C. Kirrer: Phys. Rev., 88, 951 (1952). 
(5) H. R. WoLTJER and H. K. ONNES: Comm. Phys. Lab. Leiden, no. 167c (1923). 
(6) T. R. Carver and C. P. SLICHTER: Phys. Rev., 92, 212 (1953). 
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The achievement of high nuclear polarization is obviously very important 
for performing nuclear scattering experiments. 

Only the alkali metals, in which electrons are distributed as in a Fermi 
gas, have been considered in Overhauser’s calculations. Our results are not 
based on any particular model for describing the electrons system, therefore 
it is possible to conclude that the nuclear polarization effect occurs in any 
paramagnetic substance. The reason is that the polarization is determined 
only by the fact that the heat-content of process (1) is positive and that 
of process (2) is strongly negative (*). 

Nothing instead we can tell, with thermodynamical ideas about the nuclear 
relaxation time which determines the velocity by means of which equilibrium 
is reached. We can however observe that in the case of paramagnetic non- 
metals, the process (1) requires a higher activation energy than the one re- 
quired for metals. This is because in the case of a metal a conduction electron 
in an unfilled band can always increase its kinetic energy by an amount equal 
to the magnetic energy lost in the transition (1); in a non-metal instead, where 
for the electrons only discrete energy states are available, transition (1) occurs 
. only when an energy is supplied to the electron-nucleus system in addition 
to that one released in the inversion of the magnetic moments so that the 
electron nucleus-system can jump to the next higher level. Obviously also 
the activation energy of the process (2) is increased by the same amount, so 
that the equilibrium constant of equation (5), and therefore the equilibrium 
state of nuclear polarization, remains always unchanged and not dependent 
on the type of solid. 

However the time to arrive to the equilibrium obviously depends on the 
type of solid, namely in non-metallic compounds this time must be much 
longer. 

The possibility of obtaining high nuclear polarization in non-metals is 
probably conditioned by the fact that other nuclear relaxation processes do 
not short-circuit the alignment partially; some difficulties may arise in obtaining 
paramagnetic saturation as electronic relaxation times are in some cases very 
short. 

In a future work this problem will be treated by means of statistical me- 
chanics methods. 


(*) One must consider only the heat-content of process (1), because in it the 
entropy of the system does not change. This appears from Boltzmann’s theorem, as 
the thermodynamic probability of the distribution of a certain number of electrons 
and nuclei with given spin components is found to be independent from the values 
of these components. 


_ We are SAAS to Prof. M. Cini for many helpfal discussions and to 
CIA A. GAMBA for his interest in QE work. 


RIASSUNTO 


; In questo ri si studia dal punto di vista termodinamico la polarizzazione 
nucleare in solidi paramagnetici. Si prova la generale validità dei risultati ottenuti Me 
da A. W. OVERHAUSER per i metalli alcalini. Si trova che è possibile raggiungere pola- 

izzazioni nucleari pari a quelle elettroniche nelle stesse condizioni di campo magne- 
tico e temperatura. i OA 
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Mouvement parasite vertical d’une suspension élastique symétrique 
à compensation et asservissement. 


J. E. PLAINEVAUX 


Université Libre de Bruxelles 


(ricevuto 18 Aprile 1954) 


Résumé. — Après un bref historique des dispositif proposés pour la réali- 
sation d’un guidage rectiligne sur lames élastiques, l’auteur expose une 
méthode approchée de caleul des suspensions symétriques et symétri- 
quement chargées, du type à asservissement. La méthode, valable pour 
des suspensions dont les lames sont de longueurs identiques, utilise les 
relations de Navier-Bresse, linéarisant ainsi le problème. La méthode 
suppose done de faibles charges de platine et de cadre raidisseur. Les 
formules sont données pour les deux cas principaux: lames élastiques 
sur toute leur hauteur et lames symétriques raidies dans leur partie centrale. 


Dans un précédent travail (1) nous avons décrit une méthode générale de 
calcul du mouvement parasite de dispositifs symétriques de guidage sur lames 
élastiques. 

Cette méthode s’applique aux montages dits 4 compensation qui utilisent 
la combinaison d’un montage en U et d’un montage en U renversé. Nous 
avons attiré attention sur le fait que dans le cas de dispositifs à compensation 
la précision des calculs intermédiaires devait étre fort poussée pour obtenir 
une précision acceptable dans le résultat final. Cette nécessité résulte de ce 
que les déplacements verticaux sont obtenus par la différence de valeurs trés 
voisines. De plus ces calculs doivent étre effectués au moyen de tables de 
fonctions elliptiques dont la précision est limitée en pratique. L’exécution de 


ces calculs est aussi assez longue. C’est pourquoi ils rebutent souvent ceux qui 
ont a les faire. 


(4) J. E. PLAINEVAUX: Nuovo Cimento, 10, 1451 (1953). 
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En présence de ces différents faits, il nous a semblé utile de chercher une 
autre voie d’attaque du probléme du mouvement parasite vertical. La méthode 
que nous présentons ici se limite au cas signalé à la fin de la note précitée: 
celui des systémes & compensation ot les déplacements du cadre raidisseur 
sont astreints a étre deux fois plus petits que ceux de la platine du microscope 
(Systèmes à asservissement) et-où les lames de la suspension en U sont de la 
même longueur que celles de la suspension en U renversé. 

Les dispositifs compensés, où les déplacements respectifs du cadre raidis- 
seur et de la platine résultent de l’équilibre statique des lames déformées, sont 
fort anciens. Il n’en est pas de méme du dispositif 4 asservissement géométrique 
des déplacements. 

Parmi les dispositifs compensés décrits depuis fort longtemps citons le type 
proposé par CLAY dans son adresse présidentielle à la Société Royale de Micro- 
scopie de Grande Bretagne, intitulée:-The mechanical development of the micro- 
scope (*) et en particulier le paragraphe: « A straight-line motion made with 
springs ». 

Citons aussi un article récent de J ONES (°) intéressant a divers points de vue. 
Cet article compare plusieurs dispositifs. Il attire l’attention sur les variations 
des erreurs par rapport à la ligne droite; variations dies à la non identité des 
lames. A ce propos il décrit un dispositif aussi simple qu’élégant permettant 
une égalisation des lames. Cet article confirme ce que nous pensons depuis 
longtemps: s’il est vrai que les suspensions à lames sont de réalisation écono- 
mique et ne nécessitent pas de machines-outils de haute précision pour leur 
fabrication, il n’en est pas moins vrai que le réglage de ces lames est long et 
pénible, et nécessite un personnel hautement qualifié et adroit si l’on souhaite 
obtenir une suspension de haute qualité. 

La solution consistant à asservir les déplacements de la platine 
cadre raidisseur semble intéressante à plusieurs points de vue. Si les lames 
sont remplacées par des bielles articulées à leurs extrémités, le déplacement 
de la platine est rigoureusement rectiligne si les articulations et le dispositif 
de duplication sont parfaits. Une erreur dans le dispositif de duplication des 
déplacements n’influence, comme nous l’avons montré, que faiblement le mou- 
vement vertical de la platine. Enfin, un tel dispositif est insensible à un charge- 
ment dissymétrique de la platine. Pour éviter les inconvénients inhérents aux 
articulations d’extrémités des bielles, on est conduit à remplacer les paliers 
par des articulations sur ressorts. Ces articulations se réalisent le plus sim- 
plement lorsqu'elles sont exécutées sous la forme d’une simple lame courte: 
d’où la réalisation de la bielle sous forme d’une lame de ressort raidie dans 
sa partie centrale. Enfin de toutes les exécutions possibles, la plus simple est 


x 


à ceux du 


(2) R. S. Cray: Journ. Roy. -Micr. Soc. (March 1937). 
(3) R. V. Jones: Journ. Sci. Instr., 28, 38 (1951). 
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celle où la lame de ressort n’est pas rigide en sa partie médiane. Notre étude 
portera sur les deux derniers cas: lame élastique travaillant en flexion sur 
toute sa hauteur et lame élastique rendue rigide dans sa partie centrale. 

Remarquons qu’au cours des simplifications de construction que nous venons 
d’énumérer nous ne nous sommes pas préoccupés de savoir si les avantages 
de la suspension à asservissement se conservaient. Aussi est-ce bien un peu 
sentimentalement que l’on réclame des avantages pour cette suspension. Nous 
étudierons cette question de plus près dans une note ultérieure. Bornons-nous 
à dire ici que cette note montrera que la suspension à asservissement présente 
des avantages indéniables sur la suspension non asservie. 

Disons encore ici que les dispositifs à asservissement sont plus faciles à 
étudier théoriquement. Cette simplification théorique résulte de ce que les 
lames de la suspension sont étudiables séparément, ce qui n’est pas le cas des 
suspensions non asservies. 

Signalons aussi que, a égalité de rigidité de suspension, l’effort subi par 
la vis est sensiblement le méme quel que soit le type de suspension. Cet effort 
est d’ailleurs essentiellement variable avec la position de la platine et change 
de sens au passage de cette dernière par la position moyenne. Il faut recon- 
naitre que cet effort variable est particulièrement gênant, d’autant plus que 
l’on a tendance à augmenter la rigidité pour réduire importance du mou- 
vement parasite dans le plan vertical. On est loin de travailler a pression faible 
et constante sur la vis comme il serait, souhaitable de le faire. A notre con- 
naissance il n’a pas encore été trouvé de solution satisfaisante a ce probleme 

trés important pour la lon- 
Platine gévité de la précision de la 
= vis ainsi que pour la pré- 
cision du guidage fourni par 

la suspension. 

Considérons (fig. 1), un 
dispositif à asservissement à 
lames alternées et où l’as- 
servissement est obtenu, 
pour fixer les idées, par 
deux vis dont les pas sont 
dans le rapport de un à deux. Déplacons la platine à partir de sa position 
médiane d’une longueur 2x,. Chaque lame de suspension présente alors une 
flèche horizontale x,, les tangentes aux extrémités de la lame restant verticales. 
Pour calculer le mouvement parasite dans le plan vertical, calculons les défor- 
mations en deux étapes. 

Supposons d’abord toutes les forces verticales nulles et déformons les lames 
comme il vient d’être dit. Pour obtenir cette déformation (fig. 2) nous devons 
introduire aux extrémités deux forces F, égales et de signes contraires et deux 
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moments d'encastrement M, égaux et de sens opposés. Dans cette déformation 
les lames se sont raccourcies dans le sens vertical d’une certaine quantité qui 
nous est indifférente puisque la platine n’est pas descendue par rapport au 
bâti dans cette déformation. 

Etudions une lame déformée (fig. 2). Nous supposons que les déformations 
sont dies uniquement au moment fléchissant. Admettons l'hypothèse des 
petites déformations, habituelle en résistance des matériaux. On assimile done 
la courbure de la lame déformée à la dérivée seconde de la flèche par rapport 
à l’axe de la pièce non déformée. L’arc de lame déformée est alors assimilé 
à sa projection sur l’axe Oz. Pour ce calcul, la lame ne se 
raccourcit done pas suivant son axe non déformé Oz. 
Comme nous l’avons déjà dit plus haut, ce raccourcis- 
sement ne présente pas d’intérét dans le cas qui nous 
occupe car il n’influence pas le mouvement vertical de la 
platine grace à la compensation du système de suspension. 

Tenant compte que M,— + FL, l’équation différen- 
tielle de la courbe déformée est: 


dx F, (52 : 
de ET\2 


Fig. 2. 


I étant le moment d’inertie géométrique de la lame par 
rapport à sa fibre neutre, et £ le module d’élasticité du 
matériau dont la lame est faite. Au besoin, dans cette équation ainsi que dans 
les équations suivantes, il faudra remplacer Æ par le module réduit 
E'=E(1— »?), v étant le coefficient de Poisson. 

Tenant compte des conditions aux limites: 


da 
en 0: pesa ET, aa 
da 
en A: HZ Ue pr 
il vient, tous calculs effectués: 
BR I 1 
(1) sei). 


Appelons x, la flèche horizontale à l’extrémité A, on en déduit: 


PL 


5) FADO 


L’élimination du terme 7,/EI entre les équations (1) et (2) nous donne l’équa- 


ah Ne 
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tion de la lame déformée en fonction de sa flèche maximum 2: 


\ 6a È 1 I 
(3) Ci T3 CÒ (3 L 3 2) a 


C’est cette équation que nous utiliserons dans la suite des calculs. Remarquons 
en passant que la déformée de la lame est indépendante de sa rigidité. 

C’est à partir de la connaissance de cette déformée que nous allons cal- 
culer l’abaissement de la platine et du cadre raidisseur lorsque des forces ver- 
ticales sont présentes. Nous nous bornerons à calculer des suspensions symé- 
triques et symétriquement chargées. Etudions d’abord le rapprochement de 
la platine du cadre raidisseur. 

Nous calculerons ce rapprochement en utilisant une idée de NAVIER rela- 
tive au principe de la superposition des petites déformations. Il est bien évident 
que dans le cas qui nous occupe, ce principe ne s’applique pas en toute rigueur 
et qu’il faudrait calculer les déformations avec toutes les forces en présence. 
Mais comme nous l’avons montré, ce calcul se heurte a certaines difficultés. 

Pour calculer la suspension en question avec l’idée de Navier, considérons 
que la déformée (équation (3)) sous l’action des charges horizontales est la 
forme d’une lame à l’état neutre, c’est-à-dire que la déformée que nous 
avons déterminée est une courbe géométrique et que nous ne nous préoccupons 
pas des forces et couples qwil a été nécessaire de lui appliquer pour l’obtenir. 

Supposons done qu’à cette lame «initialement » courbe, nous appliquions 
les forces verticale P et horizontale Q ainsi que des moments d’encastrement 
destinés à maintenir les extrémités verticales. Sous l'influence des moments 
fléchissants ainsi créés, tous les éléments de la courbe se déforment. Le pro- 
blème du déplacement de chaque point est généralement inextricable, puisque 
la déformation de chaque élément modifie la position et la direction des autres 
éléments; or, c’est à partir de cette position et de cette direction que nous 
calculons les déplacements. Au contraire, le problème est simple et se ramène 
à des quadratures, si nous calculons les déplacement dûs à la déformation 
d’un élément avec les moments calculés comme si l'élément en question restait 
dans sa position initiale. Sans nous occuper du déplacement que subit un élé- 
ment, nous déterminons ce qui résulte de sa flexion pour le déplacement du 
point considéré. Nous additionnons tous les déplacements ainsi calculés; autre- 
ment dit, nous prenons pour bras de levier définitifs, les bras de levier avant 
la déformation. Ce n’est légitime que si les déplacements sont petits. 

Ce procédé de calcul est en fait une linéarisation du problème. Il n’est done 
acceptable que si les charges verticales n’ont qu’une faible influence sur les 
flèches, ou d’une manière un peu plus précise, si les charges verticales sont 
petites par rapport à la charge critique d’Euler. Nous comptons préciser cette 
limite dans un travail ultérieur, en étudiant le même problème par une toute 


> 
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autre méthode. Remarquons en passant, qu’en réappliquant le principe de 
Navier a la déformée déja obtenue par cette méthode, on créerait une méthode 
d’approximations successives permettant d’améliorer les valeurs obtenues pour 
les fléches sous l’action des charges verticales. 

L’application du principe de Navier à une pièce dont l’axe est courbe 
conduit aux formules de Bresse (+). Celles-ci, lorsqu’on ne conserve que les 
déformations dùes au moment fléchissant se réduisent pour l’extrémité supé- 
rieure de la lame è: 


M 
= | —(L—z)dz 
Aw ni (l )dz, 
0 
I 
° M 
de =— | Fe (m8, 


lorsqu’on assimile l’arc de lame à sa projection sur l’axe Oz et que l’autre 
extrémité de la lame est immobile et de direction d’encastrement fixe. Re- 
marquons qu’on n’éprouverait pas de difficulté à tenir compte de l’influence 
de l’effort normal et de l’effort tranchant. Les formules seraient seulement 
plus lourdes. 

Considérons (fig. 3) la déformée sous l’influence des forces horizontales, 
sollicitée par des charges verticale P et horizontale @ 
ainsi que par des moments d’encastrement M,. 

Maintenons l’origine O fixe et calculons les déplace- 
ments Ax et Az de A sous l’influence de cette sollicita- 
tion. Tenant compte que 


1 1 
M, = ; Pr +50L, 


les formules de Bresse donnent: 


L 
P /1 AO Al 5 
4o= {|i —2) + mr] Bae, 
0 


(4) J. A. BRESSE: Cours de Mécanique appliquée a l'École impériale des Ponts et 
Chaussées, 2° édit. (Paris, 1866), 1ère partie, p. 97-111. 
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qui, tenant compte de l’équation de la déformée dont on part: 


donnent, tous calculs effectués: 


rate ap Qu: 
10.81: OVO BT? 
dal te 


14051 10EI : 


Mais le dispositif étant asservi, il s’en suit que: 


c’est-à-dire, que: 


portant cette valeur de Q dans Az, il vient: 


Pai L 
700 EI ’ 


le signe négatif indiquant que le point À descend, ce qui est évident à priori. 
Dans ce calcul nous avons négligé l’influence des forces de frottement lors 
de la descente du point A. Un simple coup d’ceil à la fig. 1 montre que cette 
simplification est parfaitement licite. 
Cela étant, soit 2P le poids de la platine; d’après ce que l’on vient de cal- 
culer la platine se rapproche du cadre raidisseur de 


PEL 
700 EI ? 


lorsqu'elle est écartée de sa position médiane d’une longueur 2x,. Si 2p repré- 
sente le poids du cadre raidisseur, la charge sur une lame supportant le cadre 
raidisseur est (P +») et sous cette influence (stabilisatrice d’ailleurs), le cadre 
raidisseur descend de 


(P + p)aiL 
700EI ? 


cone 
+ 
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donnant lieu à une descente de la platine par rapport au bâti de: 


(2P + phiL 


DÉOURE ST 


pour un déplacement 2x, de la platine, compté à partir de la position médiane. 
Les formules que nous venons d'établir permettent de calculer le mouve- 
ment parasite d’une telle platine. Aux approximations que nous avons faites 


x 


la trajectoire parasite est une parabole à axe vertical d’équation: 


_ GP + peed 
2800 HI ? 


si È représente le déplacement de la platine a partir de la position médiane; 
n représentant la descente de la platine à partir de sa position la plus haute 
ayant lieu par ailleurs à mi-course. 

Si les moments d’inertie des lames supportant platine et cadre raidisseur 
different, soit: 


I le moment d’inertie des lames supportant la platine; 


I' le moment d’inertie des lames supportant le cadre raidisseur. 


La trajectoire parasite a alors pour équation: 


e ao A 
iii ly REISER DEA TO o 


Pour terminer la question des lames de section constante sur toute la hauteur, 
il reste encore quelques mots a dire sur l’effort de poussée de la vis et sur les 
tensions qui règnent dans les lames. 

L’effort de poussée se calcule aisément puisqu’il est la résultante des forces 
F, et Q calculées précédemment. 

Quant à la tension maximum dans les lames, elle se calcule immédiatement 
par la connaissance du moment fléchissant maximum, qui se produit aux en- 
castrements. Ce moment maximum vaut, pour une lame supportant le cadre 
raidisseur : 


F 1 I 
Mn = 5 Pols 10 (P + p)a,, 


ou encore: 


634 : J. E. PLAINEVAUX 


Tandis que pour une lame supportant la platine, il vaut: 


6EI I 
ME, fa ES FE 10 Pri 


Il nous reste à étudier le cas où la lame est raidie dans sa partie centrale, 
Les lames que nous considérons possè- 
dent une partie centrale de longueur 24 
et de rigidité infinie, encadrée par deux 
parties flexibles de longueur L ayant I 
comme moment d'inertie géométrique 
par rapport à la fibre neutre. Une lame 
de suspension a done 22+2L pour lon- 
gueur totale. 
= L’étude d’une telle lame fléchie se | 
fait, en principe, comme celle d’une | 
lame de section constante. Pourtant il | 
x, xprsin? 5) est plus facile en pratique d’étudier une | 
demi-lame. 
Sous lVinfluence des seules forces | 
horizontales et couples d’encastrements | 
(fig. 4), l'équation différentielle de la déformée de la partie flexible de la lame est: 


2 


SX cos Pile ESTA 


dy Fs 
de? EI 


(Acosp+L— 2), 


en se limitant toujours à l’hypothèse des petites déformations de la partie 
flexible. 
Tenant compte des conditions aux limites: 


daz 

en Ore = 07 se = ON a = 0, 
dz 
da 

Ona By aly a, (89: 


il vient pour la partie flexible, tous calculs effectués: 


__tgy 3(4 cosp+ L)e®—- 28 
Boe 2A cos@ + Li 


ou: 
3(A cosy + Le? — 2 
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z P 
L+ Xc05 9 Q A 


avec: 


Ltgg 3A cosg + 2L 
Li — = is 
3 2%icosp+L 


Dans ces diverses expressions nous lais- 
serons l’angle g comme paramètre. 

Appliquons maintenant (fig. 5) les 
forces verticale P et horizontale @Q 
comme pour les lames d'épaisseur con- 
stante. 

L'utilisation des formules de Bresse 
conduit aux expressions suivantes du 
déplacement du point A: 
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L 
Ax, = zi] [Pte — @+Asing) + Q(L—2 + 7 cos p)|(L— 2 + A cosq) de, 
: 0 


L 
Ag, — à] ea; —x+ A4sing) + Q(L—2+7cosg)](e — x + Asing) dz, 
0 


qui, en tenant compte de l’équation de l’élastique calculé plus haut donne, 


tous calculs effectués: 


4 
62° sing cos? @ + 84?L sing cosy + 4412 sing + 5 L° te 9 


E: PE 
ESS ETE, 3(22 cos gp + L) | 
L 
ae o (A COS? | ab: 31L cos @ Li I») 
PL 
Agi Sa 


204 : 68 68 
364‘sin*?g cos?@ +604*L sing cos g + a 2 L?sin?o + 3 AL? singtgg + 35 L' tg? 


(24 cosp + L)? 


QL 


È 4 
645 sing cos? pg + 84?L sing cosy + 4AL? sing + 5 I tg 


SEL 2? cosy + L 


Mais le dispositif étant asservi il s’en suit que 4x, = 0, ce qui donne la 


636 J. E. PLAINEVAUX 
relation suivante entre les forces horizontale @ et verticale P: 


4 
Gus + 8u?L + 4uL? + 2 L3 


+ QuiL + bul? Ts’ 
où l’on a posé u= À cos y. 
Portant dans Az, la valeur de Q ainsi trouvée, on obtient, tous calculs 


effectués: 


o) 4e 
J2u+ + 12u3L 4 = wi = pls IA 


PL te? 
45EI = 1244 + 24u°L + 19u2L? + Tu L$ + Le 


AZ, = 


Si la longueur L de la partie flexible de la lame est petite par rapport a la 
longueur u — Àcosg on aura intérêt à utiliser le développement suivant: 


1 (PE) I° __25(1)° __1823(1)* 
nu) | 84\L 42\L} È 5040\w 


Le calcul du déplacement horizontal de la platine se fait sans difficulté. 


PIS tg *p 


ira 


Ù 


Il vaut en effet quatre fois l’abscisse x, du point A. Or, on a: 


62? cos? g + 64L cosp + 25° — ie 61° + 6uL + 20? 


a 
non Ier 64 cosy + 3L NACRE 


Si le rapport L/u est petit, on préférera utiliser le développement: 


io ee 
ANR ‘2\w 12\u De ANT Reta, 


Appelant encore 2p le poids du cadre raidisseur, on peut étudier de la méme 


- manière les lames supportant le cadre raidisseur et combiner les expressions 


ainsi obtenues. 

Cela va nous conduire aux expressions ci-dessous, où É représente le dépla- 
cement de la platine à partir de sa position médiane, 7 la descente de la 
platine 4 partir du point qu’elle occupait dans la position médiane: 


et Ey Ose 
MONT APN da ri 


36 
12u* + 12u3L + 7 u?L? + = ul + x LA 


Gu? + 6uL + 20? ee 
6u + 3L INA 


é—4tgp 


Re 2(2P+ PL”, È; 
TT TAGRI © Ÿ 1206 SausL + 19941 + Ual + Li 


i NOIA TA ge sea _-1823/1\! 
Seamer = ee a] 84\u)  42\u) 7 B040\n) "|" 


g étant un paramètre, et u = À cos y. 


i) ee ees see 3 n AR 6) eat Te ke ee ee ke na De (tela por” ta À AE - ln by gill 
I Se it 14 = = j + È ru M > 
SES y ' à TA i - 
- * 
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On constate que la descente de la platine est sensiblement proportionnelle 
au cube L* de la longueur libre de la lame. Il y a donc intérêt à réduire cette 
longueur le plus possible. Remarquons encore que l’angle g n’est jamais très 
grand dans les dispositifs réalisés; il vaut tout au plus 10 degrés. Pour beau- 
coup de dispositifs réalisés on constatera que la seconde fraction de l’expres- 
sion de 7 reste sensiblement: constante. 

Dans les expressions qui précédent on a supposé que les lames de support 
du cadre raidisseur étaient identiques 4 celles du support de platine. Au cas 
ou les lames supportant la platine ont I pour moment d’inertie, alors que 
celui des lames du cadre raidisseur est I’, on remplacera dans l’expression 
précédente 


par 


Ici, comme précédemment, le calcul de la poussée axiale de la vis se fait 
sans difficulté. Il suffit de se souvenir que cet effort est la résultante des 
forces Q, calculée précédemment et F, que l’on montre facilement être égale à: 


2EI tg @ 


1 = - . 
>  L(2Acosp + L) 


Quant à la tension régnant dans les lames, elle se calcule immédiatement par 
la connaissance du moment fléchissant maximum qui se produit aux encas- 
trements. 

Ce moment vaut, pour une lame supportant la platine: 


12 2 
Gu + 6u?L + — ul? + = Is 
SS 2EI(u tJ) te 9 oa fl 5 JU 5 L sia 5 
ae L(2u + L) BF Gus + 97D + Sul? + DB 


Pour une lame supportant le cadre raidisseur,-il suffira, dans cette expres- 
sion de changer I en I’ et P en —(P +9). 


L'auteur tient à remercier tout particulièrement le professeur G. P. S. 
OCCHIALINI de l’intérét qu’il n’a cessé de porter à ce travail. Il lui exprime toute 
sa gratitude pour l’avoir conduit à préciser sa pensée lors de nombreuses 
discussions aussi amicales qu’instructives. 
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RIASSUNTO (*) 


Dopo una breve rassegna storica dei dispositivi proposti per la realizzazione di 
una guidatura rettilinea su guide elastiche, l’autore espone un metodo approssimato 
per il calcolo delle sospensioni simmetriche e simmetricamente caricate di tipo asservito. 
Il metodo, valido anche per sospensioni le cui guide non sono di lunghezze identiche, 
utilizza le relazioni di Navier-Bresse, rendendo il problema lineare. Il metodo presup- 
pone quindi carichi deboli della platina e del telaio di irrigidimento. Le formule sono 
date pei due casi principali: guide elastiche per tutta la loro altezza e guide simme- 
triche irrigidite nella parte centrale. 


(*) A cura della Redazione. 
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On the Properties of a Spectrometer 
for Linear Electron Accelerators (*). 


C. C. GROSJEAN and V. J. VANHUYSE 


Interuniversitair Instituut voor Kernwetenschappen 
Centrum van de Rijksuniversiteit - Gent, Belgio 


(ricevuto il 12 Aprile 1954) 


Summary. — In connection with a specific problem concerning an electron 
spectrometer, we have solved the relativistic equations of motion of an 
electron in an arbitrary transverse magnetic field with radial symmetry. 
A detailed discussion of the exact formulae describing the motion of 
the particles on their trajectories, leads to the following conclusions: 
1) The invariance of the relation between coil current and detected 
electron energy. 2) The detected beam current does not directly provide 
the true energy spectrum, since it is proportional to the product of the 
beam intensity per unit energy interval times an energy dependent factor. 
3) A particular property of the trajectories may possibly lead to a new 
method for a precision measurement of the electron energy in a mono- 
kinetic beam. 


1. — Introduction. 


_ It is well-known that the energy spectrum of fast electrons produced by 
a linear accelerator can be determined by deviation of the particles in a 
magnetic field. When two identical air-cored coils are placed coaxially at a 
certain distance from each other, the magnetic field in each point of the sym- 
metry plane between them is perpendicular to that plane and only function 
of the distance 7 to the centre of the system. Directing the electrons in the 
symmetry plane towards the centre, they will describe trajectories depending 
on their kinetic energy £ and the electric current J through the coils. Putting 


(*) This paper is a part of the theoretical research done in connection with a linear 
electron accelerator project. 
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two slits in an appropriate position, it is possible to detect only a narrow 
beam of deviated particles having an energy in a small region (E — AH/2, 
E + AE/2) around a well-defined energy E. When I varies linearly with time 
from zero to a maximum 1,, the average £ of the detected energy region 
changes practically from zero to a maximum %,,, according to a certain law. 
Synchronizing the sweep of a cathode ray tube with the voltage from the 
field coils, and making the vertical deflection of the spot proportional to the 
amplified voltage produced by the detected beam current, one obtains a graph 
which is generally regarded as the energy spectrum. Of course, the /-scale 
has to be replaced by the H-scale and this requires the knowledge of the exact 
relationship between coil current and energy. This problem is usually sim- 
plified by replacement of the existing magnetic field by an idealized one being 
homogeneous in a circular region and zero elsewhere (1). After some simple 
calculations, one finds 


(1) dee CVE(E + 2m ec?) + 2m00*) , 


in which C is a proportionality factor. In reality, however, the exact field in 
the symmetry plane of the coils is a rather complicated function of r even 
changing sign at a certain value of r not very different from the average radius 
of the coils. Therefore, it seems reasonable to examine in how far (1) will 
be modified when the theory is worked out taking the existing field into account. 
The reason for writing this article is that the exact solution of the problem 
leads to some interesting conclusions which are, to the authors’ knowledge, 
not widely known. 


2. — Theoretical Development. 


Let H(r) describe a magnetic field with radial symmetry around a point 
O and perpendicular to the Oxy-plane. H > 0 corresponds to a field pointing 
in the upward direction. Suppose an electron comes in with an initial velo- 
city %, along the x-axis moving in the negative direction (Fig. 1). The rela- 
tivistic equations of motion are 


(2) eH 


(1) J. VASTEL: Ann. de Radioélectr., 7, n. 27 (Janv. 1952); A. SARAZIN: Ann. de 
Radioélectr., 8, n. 33 (Juill. 1953); J. G. THomasoN: A spectrometer for the 4-MeV lin. 
electr. accelerator, A.E.R.E. Report N. G/R217 (Harwell, 1948). 
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in which v represents the electron velocity, m the electron mass, e the electron 


charge and c the velocity of light. Multiplying on both sides by v, one can 
deduce the energy equation 


(3) “PES 
dé \/ v 2 
ne (2) 
Gs 
proving that 
(4) Vv=% Fig. 1. 


must be satisfied during the motion. Introducing polar coordinates, (2) gives 
rise to | 


(5) TER 270) = —<rbH(r) , 
? Mo Là nes __ Vo 
(6) FG (70 + ri) =< FH () (4 x aE 


which we have to solve with the initial conditions: 
(7) pw Lodi “6=0, (sind, =0, few 6=0. 


In these equations, e represents the positive value of the electron charge. 
The solution of the problem can further be based on eq. (6) only, since a 
suitable combination of (5) and (6) again leads to (4). Multiplying (6) by 7, 
we can write 


dre EVE DE dr 

SEE Hi mss 
(8) dt i Moe =i dt 
and after integration 

VI 
(9) TUE i È fear, 


r 


since (r20),_.—0. The velocity component v, can therefore be calculated 
in every point of the plane by the use of 


(10) | fi ro = — evi Ri J rate ar ; 


mcr 


r 


42 — Il Nuovo Cimento, 
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Now there are two cases to be considered: 


1) |r6|, calculated jrom (10), is always smaller than vp. 
In this case, it can be shown that the electron passes through the centre 0. 
Indeed, let » converge to zero in (10). If we had 


foo} 


CLL) fra di .0/, 


0 


we see that rd would become almost proportional to r-1 in the neighbourhood 
of the origin. Hence there would exist a value of r for which |rf| would be 
equal to v in contradiction to our supposition. Therefore, in the present 
case, we must have 


(12) | Yr Hr’) Gr =0 
0 
and 
esa 
(13) Ti ea 4, 
T->0 MoC 70 


since H(0) is always finite, in reality. 

As the particle progresses on its path, the radial velocity component v, 
diminishes in favor of v,, but it remains directed towards the centre O. In 
the neighbourhood of O, v, increases again as v, tends to zero and the particle 
passes exactly through the centre. From that instant on, the electron describes 
the second half of its trajectory which is symmetric with the first half with 
respect to the normal through 0. 


2) |rÒ| becomes equal to or larger than v%. 
Let 7, represent the largest value of r for which 


(14) rd == 
so that 
(15) |r |< % forny = 1 


Ym is the minimum distance of approach of the particle towards the origin 0. 


A) If r=r, is an ordinary point of the curve v, vs. r, the Taylor 
expansion of |r0| in a small region around r, is of the form 


(16) rd] = vo — Cr —1'm) + O[(r rm], 
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=, 


in which C, is positive according to (15). As the electron approaches Tm, We 
may write 


(17) = Vt — 76° = — Van Ci(r — rn) + (= 77] 
Cra 
; dr ; 
(18) | ie eee ie ee 
| J V 200 (rm) + OL r= Fn) à à 


r, being the distance of the electron at time t,. Choosing r, sufficiently near 
to rm, the terms represented by O[(r —r,,)?] can be made negligible, so that 
the approximate integration leads to 


BD oe Ne 
7 (19) Via la Tm Vr telat. 
Defining 
4 ce 2(7, — Tn) 
(20) tin a t, SE | VC, ? 


-(19) may be transformed into 


(21) PA 


9 (t—tn)? (Tm TETI). 
At r—7r, we have exactly 
d?r ; 
(22) di? = Voli = 0 . 


This proves that the particle reaches r,, in a finite time and that from that 
point on, it starts to describe the second half of its trajectory again sym- 
metrically with the first half. 


B) If r=r, happens to be a maximum or a point of inflexion with 
zero slope, the Taylor expansion of |r6 | will generally be of the form: 


(23) . 70 |= % — Cr — ru) + 


in which « can be any integer larger than 1 and C, > 0. In the same way as 
above, the approximate form of the trajectory in a small region (r,,, 73) of 7 
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will be given by 
(24) r— tn (Tr — 1m) exp [—V 20,0, (t — t,)] for a=2, 
(24') (7 — fm)! 4 * — (1, — Pn) 1% — V20,0, G — i (t—t,) fon 225 


Hence when «= 2, the time before the particle reaches r, is theoretically | 
infinite and the electron converges asymptotically towards r, describing an 
infinite number of revolutions around the origin 0. 

Considering a particle with initial conditions (7), it is clear that only in 
the cases [1] and [2A], one can talk about an angle of deviation namely the | 
angle between the two asymptots of the trajectories. Therefore, let us complete | 
the solution of the problem in these two cases by calculating the first half 
of the particle path. Combining (10) with 


(25) v(dt)? = (dr)? + r?(d0)? 
and eliminating dé, one obtains 
pr | 
| r'H(r')dr' | 
9 Mi Cpr? | 
(26) : me ; 
Le LS os 214 
fev B r'H(r') dr} | 


de | mcpr 


Therefore, the first half of the trajectory is described by 


di to r"H(r")ar" 
esatti È, 
= a) = % a | 3 
"OV « > roy 9 
r __ fevi—p? MT (y! de ; 
L magno r"H(r°)dr 


from which the final deviation angle @ can be calculated. 
For case [1], we find 


O = 2|6(0) |— 2ka if 2ka < 2|0(0)|< (2k + 1)x 
(28) 
O = 2kn— 2|6(0)| if (2k—1)x < 2|6(0)|< 2kha 
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and for case [2A]: 


| O = 2|0(rn)|— (2k — 1) if (2k — 1) < 2[0(r,)|< Qha 
(29) 
| O = (2k + 1)x— 2|6(rn) | if 2kn < 2|0(rn)|< (2k + 1)x 


These results are applicable to an infinitely thin electron beam which is 
perfectly directed towards the centre of the spectrometer. In reality however, 
these conditions are never exactly satisfied and in order to take this into 
account, we may try to extend the calculations to the case of a particle 
incoming with a non-zero impact parameter. This means electrons satisfying 
the initial conditions 


(30) = m= +00, %=0, (rsin0h=p, ty=—%, d=0. 
In exactly the same way as before, the problem can be solved with an 


arbitrary field H(r) and it turns out that (27) has to be replaced by the more 
general formula 


{SE fermen may 7 


DATE Mo? 
(31) 07) = Î 
Da fev1— PB? — f? 


r " Pp 21 
f | mer! È LL sti 


rl 


for the first half of the path. 


3. — Returning to the original problem, where the magnetic field is produced 
by the coil current J, we may say that 


(32) H(r) = If(r), 


in which f(r) characterizes the field distribution. Introducing this in (27) 
and (31), it is clear that a trajectory remains unchanged if the relation 


B 


—S02 
se Vip 


is satisfied. Expressed in terms of the kinetic energy Æ, (33) becomes equi- 
valent to (1). 

Now, let us consider an electron beam having a uniformly distributed energy 
spectrum throughout its cross-section. Suppose that with a current J, and 
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a particular position of the two slits, we are able to detect the electrons in 
the narrow interval (E, — AE,/2, E.+4£,/2). When the current changes to I, 
we see that the same set of trajectories passing through the slits will now be 
described by particles in a new interval (E — AE/2, E+AH/2) determined by 


TO VEN EME) 


RE i Ve a VE (E, + 2m,¢?) 
and i 

E(E + 2m,c?) 
(35) ADS (E + me?) 


AE, sw E,(E, + 2me?) 
(E, + mye?) 


This is not only valid for two slits placed at large distance from the centre 
where the trajectories become practically straight lines, but also for any 
arbitrary position of the slits. Therefore, we may conclude that the relation (1) 
between coil current and electron energy is rigorously valid for each considered 
transverse magnetic field with radial symmetry. In addition to this, (35) 
shows that the detected energy interval AH is not constant, but varies sensi- 
tively with the average energy £, proving that the detected beam is not 
directly proportional to the particle current per unit energy interval. If i(E) 
represents this current density, the vertical deviation of the cathode ray is 
proportional to 


E(E + 2mye? 


The factor by which 2(£) is multiplied causes a deformation of the real 
spectrum. Indeed, it is formed by the product of E times a factor varying 
from 2 to 1 as the energy increases from zero to infinity. Hence, a peak of 
accelerated electrons in a relatively high energy region has the appearance 
to correspond to a much larger intensity than an equivalent peak at lower 
energies. This greater efficiency of the spectrometer at higher energies is in 
fact an essential property of the apparatus. But it means that in practice 
when the graph on the cathode ray screen consists of a small peak due to 
electrons with initial energy plus a large peak of accelerated electrons, one may 
not directly conclude that nearly all the particles have been accelerated. Other 
effects due to the factor are the broadening of peaks and the slight shift of 
their maximum towards higher energies. 


È 
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4. — Application of the theory. 


We wish to apply the results of $ 2 to the case of two equal and coaxial 
air-cored coils. Let us first consider the simplified case of two circular con- 
ductors with radius À placed at a distance 27 from each other. One can easily 
calculate the field H(r) in the symmetry plane, when a current J is used: 


(37) Hm 4RI È (R —r cos y) dy 


STO ì r? — 2rR cos y + R°+ L?} 


0 


in which H is measured in oerstedt and J in ampere. The integral 


(38) | r'H(r')dr' 


can be calculated from 


r 


; FA 4RI — r' cosy) dy 
38) fran di La la 2r'E cos p+ B+ Lh 


0 


which can be reduced to 


0/2 


: 41 dy 
I nal nh ( 2 2 02 È. 5 
(40) IL H(r ) dy 5 R P+ 7 aa + L? —- 4rR sin? pl 
0 0 


RS five r)? + L? — 4rR sin? p}} dy 


by using carthesian coordinates (£, 7) instead of (r’, y) in the integration plane 
of (39) and by putting 


(41) 7 =r sin 29 , 


after integration with respect to £. Introducing the complete elliptic integrals 
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K(k) and E(k), we have 


I 


(42) | VE) dr = “TUR + n° + (1-5) K(k) — E(k) 


i 2aI TRA la 
i AOS ER Se pp 


in which K=2VrR/V(R+r?+I?. 
For r — ©, we find 


(co) 


(43) | e Air’ \dr’ = 0. 


0 


Defining the dimensionless parameters 


we obtain 


k2 2at 
es Bi) —(1—F Kb) (+ A+ ay) 


and from (10), we deduce 


5 Si nr 2 
TREIA ( k 


69) 10me  kyva 


In fig. 2, we have plotted 


(47) y = fla, a) o 


1 ke 
=i (1 = ) xu) 206) 


as a function of « for a few values of « The curves reach a maximum in the 
vicinity of x = 1 as it was to be expected. Let us define the current 


A0 ME ob I. : 
Se a fae a) ee Vip à 


(48) I.(a) 


Keeping in mind our previous discussion, three cases have to be considered: 
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1) I<I,. Then we know that v, is always smaller than v and the 
particles pass through the origin. The final deviation angle can be derived 
from 


fla, at) 7 | da 


(49) = nef om 
[(25mzc4p' /16e21?( (1 — py) — f(a, a) |? # 


which can easily be reduced to an angle between 0 and x (cfr. (28)). 


2) I>1I,. The distance of closest approach to the centre can be de- 
termined by intersection of 


= 10m,c’. ives 
Bel V1 8 


(50) 


with f(x, x), taking the point with the largest x-coordinate. For every value 
of a, this distance is at least of the order of the radius R, for it is always 
larger than the abscissa of the point where f(x, x) reaches its maximum. 


3)I =1I,. This is precisely the case [2B] in which r —r, happens 
to be a maximum. Theoretically the particle will describe a trajectory ap- 
proximately given by (24). In practice, however, some small perturbations 
may cause the particle to behave as in one of the preceding cases. 


Anyway, it seems a remarkable fact that the distance of closest approach 
to the centre undergoes a sharp discontinuity at the «critical» current I., 
as it is exactly zero when J is smaller than J, and at least of the order of È 
when J is larger than I,. In principle, this could lead to the possibility of 
an accurate determination of the electron velocity in a monoenergetic beam. 
An electron detector placed at the centre of our system of coaxial coils would 
first transmit the full beam intensity as long as the increasing J is lower than 
the critical current, and then drop to zero for I exceeding J,. In the case 
of a beam in which the energies are distributed as a peak with a finite width, 
the method would at least permit the measurement of the average energy. 
The relation between J, and v/c is relatively simple to calculate in each special 
case using (48), which shows the same energy dependence as (33). For instance, 
applying (48) in the case of L= R/2 or «x — 1/2, we find 


LOTO (B10) TUE 2 Lee 
8 x 4.80 10-10 0.443 V1 — p "Vip 


(Ol) Ve 


It is clear that the simplicity of this calculation would make the method 
particularly interesting for the absolute calibration of a monoenergetic beam. 
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All these considerations apply to the case of an infinitely thin and perfectly 
centered electron beam deviated in the magnetic field of two coaxial circular 
conductors. But in practice, the problem is somewhat more involved due to 
the following factors: 


1) The field is produced by two coils. This makes the application of 

the formulae slightly more complicated but the curves replacing those of Fig. 2 
will show essentially the same behaviour 

and all our conclusions will remain unal- 


"04 tered. In the case of two coils at an ave- 
rage distance equal to half their average 
03 diameter, (51) provides approximately 


the number of ampere-turns required. 


0.2 2) The particle beam is not ideally 
thin. This necessitates an analysis of the 


0.1 behaviour of particles incoming with a 
Paie small impact parameter p. It can easily 
be shown that the distance of closest ap- 


RE Soe On a SE GA ye proach to the centre remains practically 

Fig 2. equal to p as long as I is smaller than a 

critical current and shows the same dis- 

continuity as before when J exceeds -this current. The latter is slightly 

dependent on p and will therefore not be rigorously equal to I, calculated 

for p= 0. This will only have the effect of making the discontinuity in the 
detected current less abrupt, but it must still be very clearly observable. 


It is a pleasure to thank the Institut Interuniversitaire des Sciences Nuclé- 
aires and Prof. Dr. J. L. VERHAEGHE for his interest in this work. 


RIASSUNTO (*) 


In connessione con un problema specifico concernente uno spettrometro elettronico, 
abbiamo risolto le equazioni relativistiche del moto di un elettrone in un campo magne- 
tico trasversale arbitrario a simmetria radiale. La discussione dettagliata delle formole 
esatte che descrivono il moto delle particelle lungo le loro traiettorie porta alle conclu- 
sioni seguenti: 1) la relazione tra la corrente di campo e l’energia elettronica rivelata 
è invariante; 2) la corrente elettronica non fornisce direttamente il vero spettro ener- 
getico, poichè è proporzionale al prodotto dell’intensità della corrente per intervallo 
unitario dell’energia moltiplicata per un fattore dipendente dall’energia; 3) una pro- 
prietà particolare delle traiettorie può forse condurre a un nuovo metodo per la misura 
di precisione dell’energia elettronica in un fascio monocinetico. 


(*) Traduzione a cura della Redazione. 
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An Amplifier with Low Dissipation and Short 
Rise-Time for C.R.T. (*). 


D. BRINI, L. PELI, O. RimonpI and P. VERONESI 


Istituto di Fisica dell’ Università - Bologna 


(ricevuto il 15 Aprile 1954) 


Summary. — A description is made here of a C.R.T. Amplifier for ne- 
gative input, set up with secondary-emission tubes, whose dissipation 
was 55 mA, gain 400 and rise time not greater than 7 mus. 


In the last few years scintillation counters have reached a high degree 
of development and have been used in many fields of nuclear physics. The 
counters with organic liquids or crystals give pulses with a short rise time. 
This has necessitated a study of the appropriate electronic techniques. 

The presentation of pulses to a C.R.T. requires amplifiers whose funda- 
mental requirements are as follows: 


a) short rise time; 
b) high level output. 


Point a) is in practice limited, for conventional types of linear amplifiers, 
by the characteristics of used tubes, namely by their figure of merit F' which, 
as is known, is represented by 


(where C, is usually the sum of the input capacity with the tube output ca- 
pacity). 


(*) This work has been partially supported by the Padua Section of the I.N.F.N.. 
(1) VALLEY WALMAN: Vacuum Tube Amplifiers (M.I.T.), p. 73. 
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For the interstage tubes which are more currently used in wide-band am- 
plifiers, the figures of merit are as shown in Table I. 


TABLE I. 
g (mA/V) OC, (pF) F (1/us) Dissipation (mA) 
6AC7 9 . 19 230 13 
6CB6 6.2 10 280 125 
6AH6 9 14 290 13 
6AK5 012 8 300 10 


The figures of merit of the tubes which are commonly used as the push-pull 
output stage are as shown in Table II. 


TABLE II. ' 
g (mA/V) | C, (pF) F (1/us) Dissipation (mA) 
832 A 20 | 16 120 150 
two 807 6 21 130 200 
829 B 8.5 | LE 170 200 


The conventional amplifiers utilizing about 10 tubes, which dissipate 200 
to 300 mA.and have a gain of the order of 100 with a reasonable output for 
presentation to a C.R.T., can reach rise times about. 100 mus. 

With critical quadrupole compensation, such amplifiers can reach rise times 
not smaller than 20 to 30 mus (?). 

A substantial improvement can be obtained with the use of distributed 
amplifiers, which have been developed in the last few yers (*). Rise times 
which can be reached with such amplifiers are about 2 to 3 mus. This neces- 
sitates, however, the use of a large number of tubes, which includes a high 
dissipation of current and high costs. 

It is often required to have good rise times without arriving at those which 
are realized with distributed amplifiers. This can be obtained by using the 
secondary-emission tubes (4) such as the EFP60, recently put on the market 
by Philips. 


2 


(2) As (1), page 110. 

(3) G. F. Myers: Electronics, Jan., 128 (1952); G. C. KELLEY: Rev. Scient. Instr., 
71 (1950). , 

(4) C. W. MULLER: Proc. I.R.E., 39, 159 (1950). 
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By reasons of their secondary emission, these tubes possess great trans- 
conductance, which greatly increases their figure of merit. Table III reports 
the EF P60 figures of merit, where the outputs taken on the plate and on the 
multiplier dynode respectively were used. | 


TABLE III. 
| 
| g (mA/V) | ©, (pF) | F (1/us) Dissipation (mA) 
Plate | 25 | 17 670 22 
Dynode | 20 | 22 | 410 Î 22 


For signals of single polarity, the EFP60, however, for sufficiently high 
peak current, can also be used as output tubes, obtaining two outputs in 
push-pull, one on the plate and one on the multiplier dynode. 

On the basis of these observations we have realized an amplifier which 
satisfies following requirements: 


) input for negative pulses; 

) gain of about 400; 

) total rise time of 5 to 7 mus; 

) output of about 45 V (in order to obtain on a 5XP C.R.T. a deflexion 
of about 1 inch); 

e) low current dissipation (55 mA). 


a 
b 
e 


l 


J 
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Fig. 1. — All the resistors are in kQ. All the capacitors are in yF. All the induc- 
tors are in pH. 


The scheme is shown in Fig. 1. As will be seen it is a conventional scheme 
with quadruple critical compensation. The dynode load in the output stage 
is less than the plate one, notwithstanding the smaller transconductance of 


TE ad 47 at 7 dm 
Ca > 
654 D. BRINI, L. PELI, O, RIMONDI and P. VERONESI 


the tubes relative to the dynode. This is with the aim of obtaining equal 
rise times for the two outputs in push-pull. For this reason we obtained as 
maximum outputs about 30 V on the plate and about 15 on the dynode. 

The rise-time measurements were made by the oscillographic method, with 
pulses obtained by a mercury relay (of the Westinghouse 275B type). Pulses 
of a known length were used for the calibration of the syncroscope sweep. 
Various photographs were taken with different pulse-shaping times, the rise 
times being calculated therefrom. The measurements show that the rise time 
is certainly between 5 and 7 mys. 

The gain, measured with slow calibrated pulses, was about 400, and the 
output (linear) 45 V. 

These data proved to be in agreement with those calculated. 

In the practical setting-up of the amplifier great care had to be taken in 
the search for earth points for each stage, in order to avoid self-oscillations. 

Both the filtering of the alternating current heater and the electrolitical 
condensers by-passed with ceramic capacities were found to be convenient. 


Thanks are due to Mr. R. DEGLI Esposti for his collaboration in the setting 
up of the circuit. 


RIASSUNTO (*) 


Si descrive un amplificatore C.R.T. per impulsi di entrata negativi, costruito con 
tubi ad emissione secondaria, con dissipazione di 55 mA, guadagno 400 e tempo di 
salita non superiore a 7 mus. 


(*) Traduzione a cura della Redazione. 


Study of a Liquid Scintillation Counter (*). 


D. Brini, L, PELI, O. RIMONDI and P. VERONESI 


Istituto di Fisica dell’ Università - Bologna 


(ricevuto il 15 Aprile 1954) 


Summary. — This paper deals with a liquid scintillation counter of small 
thickness (1 cm) and of circular cross-section (10 cm in diameter). The 
measurements have shown that its efficiency in revealing relativistic 
particles is practically 100%; and that its real diameter corresponds to 
the geometric one. The measurements permit the conclusion that the 
collection of light is about 30% of that emitted in scintillation, with a 
disuniformity in the light collected not greater than 30%. 


To count the relativistic charged particles a liquid scintillation counter 
was constructed, which, in addition to the known property of rapidity, pre- 
sented the following characteristics: 


a) a high efficiency over the whole sensible area; 

b) a cylindrical geometrical form of small thickness; 

c) a minimum thickness in the direction of the path of the particles 
(that of the cylinder axis). 


The design of the apparatus used is shown in Fig. 1. 

A vessel (1), containing a solution (3 g lof p-terphenyl+10 mg of diphenyl- 
hexatriene in 1 litre of phenylcyclohexane) was constructed from a sheet of 
plexiglass, 12 mm in thickness, and closed in a cylindrical box (2) (3) made 
from a sheet of brass of 0.2 mm in thickness, the inside walls of which were 
varnished with MgO. On one side of the box there was a cylindrical tube (4), 
which contained a photomultiplier 5819 (6), in optical contact with the side 
wall of the vessel, which, in such a position, was flat. Optical contact was 


(*) This work was partially supported by the Padua Section of the I.N.F.N.. 
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realized with silicon oil D.C. 200 placed between the above-mentioned wall 
and the plexiglass dise (5); (9), between this disc and the photomultiplier 
head we put vaseline oil held in by rubber gaskets (9) between the photo- 
multiplier and the metal tube (4). 

The photomultiplier was attached 
laterally in order to satisfy condit- 
ion c) of the project. The type of 
optical collection used was in accor- 
dance with the prime necessity for 
obtaining the greatest collection of 
light, even at the risk of a loss of 
uniformity and rapidity, which would 
suggest the use of a different optical 
collection (1). 

Using the following hypotheses 
for purposes of simplification, ,we 
may calculate light- collection as 
under: 


beer LU: x 6 > 


LIZ, 


a) the cylindrical box was conceived as a white sphere containing the 
vessel in its maximum section; 


b) the index of refraction of the material, of which the vessel was cons- 
tructed (n,—1.49) was the same as that of the scintillator (n,—1.52); 


c) the opening made on the sphere was small in respect to its diameter; 
d) scintillation was in point form; 


e) the body was perfectly transparent. 
Total light collection was due to the following contributions: 


1) the fraction of light which arrived on the photomultiplier directly 
from the point of scintillation; 


2) the fraction of light which was entrapped in the vessel by total re- 
flection, and which partly reached the photomultiplier; 


3) the fraction of light partly diffused by the white vall. 


The calculation (?) (*) of the second contribution showed that the quantity 
of light which remained entrapped in the vessel was a growing function of 


(1) R. L. GARWIN: Rev. Scient. Instr., 23, 255 (1952). 
(2) W. A. SHURCLIFF and R. C. Jones: Journ. Opt. Soc. of Am., 39, 912 (1949). 
(*) See Appendix. 
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the distance from the centre of the point at which scintillation took place. 
In fact, the rays making an angle (see Fig. 2) between the values y, and ¢, 


defined by 


Va T. 


1 


1 
COS Pi == ; COS Y, = 
cd i n 


Vi &? sen? 0 


with k— d/r, 


remained entrapped, since they made an angle greater than the critical angle. 
This light, after successive reflections, reached the anterior face of the light 


conductor 5 (Fig. 1). Not all of it, 
however, reached the photomulti- 
plier, because part of if escaped la- 
terally, owing to the angles of inci- 
dence. 

In Fig. 3 the curve A, represents 
contributions 1) and 2) for the dia- 
meter perpendicular to the photo- 
multiplier axis. The light diffused 
by the white wall is not present. 

The curve A, represents, on the 
other hand, the percentage of the 
total flux emitted in scintillation 
and collected by the photomultiplier, 
with the addition of contribution 3). 
The calculation of the light diffused 
was made with the formula of Ul- 
bricht’s sphere with a diffusion co- 
efficient {= 0.9. 


Fig. 2. 


The other two curves refer to the diameter B (photomultiplier axis). The 
greater disuniformity in the collection depends evidently on the notable in- 


—— Diameter A with black walls(An) 
i UD RTS) 


» A » white » (Ab) 
» Bs » » (Bb) 


Fig. 3. 


43 — Il Nuovo Cimento. 


crease of contribution 1) when 
scintillation took place near the 
photomultiplier. 

The substitution of the white 
sphere with a cylinder, leads, for 
reasons of symmetry and im- 
perfect diffusion, to a superval- 
uation in the calculations made, 
of the contribution of the diffu- 
sion for those points far from the 
cylinder axis: that leads, in 
practice, to a better uniformity in 
light collection. As will be seen 
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from Fig. 3, the diffusion increases light collection in the central zone 
to about 70%. 

Experimental measurements of light collection were effected with a source 
of y (50 mC of Co), collimating a bundle with a lead screen of 10 cm thick, 
in which a hole of 6 mm in diameter was bored. The counter was hit by such 
a bundle at several distances from the centre, along a diameter, measurements 


o — Diameter A with black walls (An) 
» » (Bn) 


A » white » (Ab) 
B » » » (Bb) 


Fig. 4. 


being taken of the corresponding anodic current. The intensity of the current 
versus the distance of scintillation from the centre was then set down in diagram 
form (Fig. 4). 

The decrease of the intensity of the current near the edge was not due 
to a minor light collection, but to diffusion of y and to the fact than the bundle 
of y had finite dimensions. 

These contributions can be taken into account from the values of the current 
beyond 5 cm from the centre of the counter. In so doing a good agreement 
with the curves of Fig. 3 which give light collection can be noted. Regarding 
the percentual increase of light collection due to the white walls it was ob- 
served that the agreement was excellent. 

It is to be noted that the counter is not uniform; yet it can be affirmed 
that the disuniformity is not greater than 30%, although such a maximum 
disuniformity takes place only in a very limited area near the photomultiplier 
cathode. 


The efficiency proofs regarding the revelation of the relativistic particles 
were effected with the apparatus shown in Fig. 5. A telescope, formed by 
the GM counter plates 1, 2 and 3, defined a solid angle covered by 10 cm of 
lead. The photomultiplier 5819 functioned at 1000 V, the tensions being 
distributed in the proportions of 2, 1, 1, etc.. 

The threefold coincidence noted the number of particles crossing a 4 x6 cm? 
area of the scintillation counter. The fourfold coincidence showed the number 


Se IRR 
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of particles revealed by the scintillation counter too. A registration of the 
events of the threefold and fourfold coincidences made over some days gave 
the following results: threefold 1369; fourfold 1359. The efficiency was then 


0.992 + 0.0024. 

The three and fourfold chance 
coincidences were calculated, but 
proved negligible, since the fre- 
quency of each plate was rather low 
and the resolving power of the co- 
incidence about 5.107 s. 

Such an efficiency value repre- 
sents, however, a limit inferior to 
the true one, since the side showers 
produced an apparent inefficiency. 
In our case we hold that 1% of 


inefficiency may be attributed to a cause of this nature, rather than to a real 


inefficiency. 


In order to measure the minimum height of the counter pulses, the three 
to fourfold rate as a function of the bias of the integral discriminator, placed 


0 0,05 0,10 0,15 0,20 0,25 


Fig. 6. 


on the scintillation counter channel after 
a proportional amplifier, was determined. 
The result obtained is shown in Fig. 6. 
Point P on the curve corresponds instead 
to the measurement made without lead. As 
can be seen, the efficiency without a screen 
is lowered to about 98%, which justifies the 
affirmation that 1% of inefficiency may be 
attributed to side showers. 

The curve in Fig. 6 shows the integral 
pulses of the photomultiplier, due to part- 
icles collimated by the GM-counter tele- 
scope. It is noted that the maximum effi- 
ciency is obtained with a height not superior 
to about 0.03 V. 

Bearing in mind that light collection was 


about 30% of that emitted in scintillation, and considering that, in this 
case, the photomultiplier gain was 2:10°, the size of the pulses, according to 
Morton’s calculus (3), is comparable to those very small pulses found exper- 


(3) G. A. Morton: RCA Rev., 10, 525 (1949). 
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imentally, considering also the statistical fluctuations in the minimum 
energy loss (4) (1.8 MeV/g/cm?). 

In conclsuion, the measurements obtained show that the area which is 
effectively sensible corresponds to the geometrical area of the counter. In 
fact, even if the zone explored had a surface of 4x6 cm? round the centre 
of the counter, the form of the curve of light collection at the edges does not 
prejudice efficiency; the points of maximum disuniformity refer to positive 
differences in respect to the central area, which leads to an increase in the 
pulses. For all these characteristics, the counter proves useful in cases where 
well-defined geometry and a high-resolving power and minimum thickness of 
the counter are required. 


Thanks are due to Mr. A. COLLINA for his construction of the counter, and 
Mr. A. MICHELINI for help given in the measurements made. 


APPENDIX 


Tt is required to calculate the percentage of the luminous flux which remains 
entrapped in a transparent cylindrical body, having a refraction index # and 
radius r, when the source is in point form and placed inside it. 

The entrapped light is due to: 


1) a constant contribution 
independent of the source po- 
sition ; 

2) a contribution which var- 
ies with the source position. 


Given the distance d of the 
source from the cylinder axis, a 
generic ray emitted at point P, 
where the source is placed, is fixed 
by the angles œ and 0 (Fig. 2). 
The first angle is that made by the ray with the plane z passing through P 
and perpendicular to the cylinder axis; the second is the angle made by 
the plane © containing the luminous ray and perpendicular to x with the 
plane À passing through P and containing the cylinder axis. 


(4) T. BowEN and F. X. Roser: Phys. Rev., 85, 992 (1952); R. Baskin and J. R. 
WINCKLER: Phys. Rev., 92, 454 (1953); T. E. GRAUSHAW: Progress in Nuclear Physics 
(ol: 22), 2 ie 
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The contribution 1) is that due to the rays which have the angle between 
P and g, (see Fig. 7), which satisfie the equations 


si 1 ae al 
ingo ==; cosg, ==. 
O 3 n 


This is true because the rays having these values of œ fall independently 
of d and 0 on the surface of the cylinder with angles greater than the critical 
one. The value of the solid angle where these rays are included is obviously 


veti 


n 


Wo = 40 


If we assume that n = 1.5, the percentage of light which remains en- 
trapped is 7.87%. When n< 42 this contribution is nul. 

To calculate contribution 2) it is necessary to know the value of g because 
when a ray leaves P, it undergoes a first total reflection and falls thereafter 
on the surface at angles still greater than the 
critical one. This can never happen if g > m1. 
In fact, the critical angle being (7/2) -—-g,, the 
rays go out whether they are incident on the 
cylinder bases before or after total reflexions on 
the side walls. A part only of the rays which 
made the angle <q will remain entrapped, * 
since, when P departs from the cylinder axis, 
the angle of incidence on the side walls can 
become greater than the critical angle. 

Applying the sinus theorem (Fig. 8) to the 
triangle DCP, eliminating sin à and holding 
that the angle « is equal to the critical angle, we find the angle g which fixes 
the rays coming from P and incident on the surface of the cylinder with the 
critical angle. The value of this angle satisfies the equation 
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Fig. 8. 
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where % is equal to d/r. 

For this reason the rays with g:<@< % are incident with angles greater 
than the critical angle and remain entrapped. ‘The solid angle defined by 
the preceding condition is given by 


N Po n 
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The calculation of the percentage of light entrapped, for different distances, 
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reported in Table I, has been carried out with numerical integrations in the 
particular case where n = 1.5. 


TABLE I. 
= djr | 0 0.33 | 0.66 0.834 
eg le aut ip ea a see | 


RIASSUNTO (*) 


Il lavoro descrive un contatore a scintillazione a liquido di piccolo spessore (1 cm) 
e di sezione circolare (10 cm di diametro). Le misure eseguite hanno dimostrato che la 
sua efficienza alla rivelazione di particelle relativistiche è praticamente del 100%, 
e che il suo diametro reale corrisponde a quello geometrico. Le misure permettono di 
concludere che la luce raccolta è circa il 30% di quella emessa nella scintillazione, con 
una disuniformità nella raccolta della luce non superiore al 30%. 


(*) Traduzione a cura della Redazione. 
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Isotopenzusammensetzung und Radioaktivitàt 
von rezentem Vesuvblei. 


F. BEGEMANN, J. GEISS, F. G. HOUTERMANS 


Physikalisches Institut der Universität - Bern 


W. BUSER 


Anorganisch-Chemisches Institut der Universitdt - Bern 


(ricevuto il 16 Aprile 1954) 


Zusammenfassung. — Das Blei aus drei Cotunnitproben vom Vesuv 
wurde radioaktiv und massenspektrometrisch untersucht. Die Radio- 
aktivitàt wurde als zum RaD und seinen Folgeprodukten gehòrig identi- 
fiziert und ergab — auf den Auffindungszeitpunkt extrapoliert — 
21,2 + 0,5 Zerfälle/s-mg Blei. Das Cotunnitblei lift sich als ein bisher 
unbekannter «intermediärer » Bleityp auffassen, der zu etwa 0,7% aus 
radiogenem Zerfallsblei und im übrigen aus rezentem « gewdhnlichem 
Blei» besteht. Das Zerfallsblei entstammt einem « Urankonzentrat », 
dessen Maximalalter sich nach der RaD-Methode zu 17-10° a ermitteln 
lieB. Die Isotopenzusammensetzung (vgl. Tab. I) zeigt eine ähnliche 
Anomalie wie die der Bleierzlagerstatten von Joplin, Sudbury u.a. Es 
wird ein Modell fiir den Entstehungsmechanismus dieser « J-Typen » 
angegeben. 


1. — Einleitung. 


Für die Frage der Bildung von Bleimineralien und für die Theorie des 
« gewohnlichen Bleis » ist es von besonderem Interesse, die Isotopenzusam- 
mensetzung eines Bleiminerals zu messen, dessen rezente Bildung aus einem 
definierten Magma sichergestellt ist. Als solches bot sich Cotunnit, ein Subli- 
mationsprodukt vulkanischer Fumarolen. Uns standen 3 Cotunnit-Proben vom 
Vesuv zur Verfiigung. Diese zeigten eine gut meBbare Radioaktivitàt, wie sie 
fiir dieses Mineral schon von ZAMBONINI [Za 07] festgestellt worden war. 
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2. — Messungen. 


Die drei genannten Proben wurden zu verschiedenen Zeiten an drei ver- 
schiedenen Orten am Vesuv gefunden (vgl. Tabelle I). Der schon kristallisierte, 
rein weiBe Cotunnit lieB sich gut mechanisch vom roten Tuffstein ablôsen. 
Von jeder Probe standen einige hundert mg zur Verfiigung. Debye-Scherrer- 
Aufnahmen der drei Proben zeigten voéllige Identitat mit dem Diagramm von 
synthetischem PbCl,. Die chemische Analyse ergab, daB es sich bei den 
vorliegenden Cotunnitproben um sehr reines PbCl, handelte. (Der stôchio- 
metrische Bleigehalt betragt 74,5%, gefunden wurden 74,4%). 

Zur Untersuchung wurden etwa 100 mg jeder Probe aufgelôst, Beimen- 
gungen des Trägermaterials abfiltriert und aus der Lésung das Blei als Sulfid 
gefallt. Dieses diente als Ausgangsmaterial fiir die massenspektrometrische 
Untersuchung und fiir die Aktivitàtsbestimmung. Auf diese Weise wurde 
erreicht, daB die Bestimmung der Isotopenverhältnisse und die radioaktiven 
Messungen an genau gleichem Material durchgeführt wurden, da beim un- 
veränderten natürlichen Mineral mit der Môglichkeit zu rechnen war, an ver- 
schiedenen Stellen des Sublimats oder selbst innerhalb einzelner Kristalle 
Unterschiede der Radioaktivitàt oder der Bleiisotopenverhaltnisse zu finden. 
Die spezifische Aktivitàt würde z.B. variieren, wenn die Dauer der Sublimat- 
bildung vergleichbar mit der Halbwertszeit des RaD ware. 

Die Bleiisotopenverhaltnisse wurden mit einem einfach-fokussierenden Mas- 
senspektrometer mit 60°-Sektorfeld und Elektronensto8-Ionenquelle gemessen. 
Das Massenspektrometer wird a.a. O. [Eh 54] genauer beschrieben. Das 
Bleisulfid wird aus einem Ofen bei ca. 300-400 °C verdampft und die 
Isotopenverhältnisse des Bleies durch Messung der Ionenstréme der Pb+- und 
PbS+-Ionen bestimmt [Eh 53, Gei 54]. Die Ergebnisse der Messungen sind in 
Tabelle I zusammengestellt, die Isotopenverhältnisse für die drei Cotunnite 
sind innerhalb der Fehlergrenzen gleich. 


TABELLE I. 
Pb206 Pb?207 Pk 208 
Probe Funddatum oS = ea A == + 
Pb204 Pb24 Pb 04 
MECOLUNNICA RO 12.8.52 19,08 LOST 39,24 
Cotunnit Bas 3.7.49 19,14 15,75 39,37 
sacouunmit i. 24.8.52 19,19 15,80 39,62 


Da die Vermutung nahe lag, daf der grôfite Teil der Radioaktivitat vom 
RaD(Pb?!°) und seinen Folgeprodukten herrührt [Ro 07], wurde eine quan- 
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titative Abtrennung des Bleis von allen anderen Elementen nach der Dithizon- 
Methode vorgenommen (*). Die Aktivitàt der Bleifraktionen wurde in Becher- 
zahlrohren durch Vergleich mit einer RaD-Standard-Lésung [Be 52] bestimmt. 
(Gemessen wird in beiden Fallen die B-Strahlung des RaE). 

Der Bleigehalt der Lésungen wurde gravimetrisch als PbSO, bestimmt. 
Die Ergebnisse der Aktivitàtsmessungen und der Bleibestimmungen zeigt 
Tabelle II. 


(RABELLEOHE 


E Zerfalle 
Fund Pb-Gehalt Zerfalle ag, mg 
Probe a i der Losung| s-mg Pb i | Th-Aquivalent Tn 
atum | (mg/em3) al (extrapol. auf | — DD 
(Marz 54) | Punddatum) | = 
Cot. A 12.8.52 1,720 20,5 21,6 + 0,5 
Cot. B 3.7.49 1,212 18,4 21,4+0,5 | ZA 
le Cota 24.8.52 1.163 19,7 20,7 +0,5 | | 


Schon die Untersuchung der bleifreien Fraktion nach der Dithizonaus- 
schüttelung ergab, daf bei allen drei Cotunniten mehr als 99% der gesamten 
B-Aktivität dem RaD und seinen Folgeprodukten zuzuschreiben sind. Darüber 
hinaus wurden die Cotunnit-Lésungen in einer speziellen Proportional-Zahl- 
rohrapparatur, mit der die Tn- und Rn-Aktivitàten getrennt gemessen werden 
kônnen auf ihren Gehalt an RdTh bzw. Ra untersucht (Einzelheiten der Meb- 
anordnung sowie die genauen Resultate siehe [Be 54]). Die vorläufigen Er- 
gebnisse sind ebenfalls in Tabelle II aufgeführt. 


3. — Diskussion der Ergebnisse. 


Sämtliche bisher untersuchten Typen von natürlichem Blei lassen sich 
entweder dem Typus des «gewühnlichen Bleis » oder dem des «radiogenen 
Bleis » zuordnen. Der erstere Typus ist dadurch charakterisiert, daf sein gegen- 
wärtiger Uran- und Thoriumgehalt von verschwindender GròBenordnung ist 
und keinesfalls fiir das Vorhandensein eines nennenswerten Anteils der Blei- 
isotope Pb?°, Pb? und Pb? verantwortlich gemacht werden kann. Nach 
dem von verschiedenen Autoren [Ger 42, Hol 46, Hou 46, Bu 49, Co 52] gege- 
benen Modell der Entwicklung von Bleimineralien lassen sich diese auffassen 


(*) Vgl. z.B. F. J. WELCHER: Organic Analytical Reagents, Vol. 3. 
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als ein Gemisch aus einem « Urblei » einer zur Zeit der Bildung der Erdkruste 
gegebenen Isotopenzusammensetzung und Zerfallsblei, das in der Zeit von der 
Differenzierung der Erdkruste bis zur Bildung des Bleiminerals aus Uran und 
Thorium entstanden ist. 

Radiogenes Blei, das bisher ausschlieBlich in Uran- und Thoriummineralien 
oder in engster Assoziation mit diesen gefunden wurde, besteht vorwiegend 
aus den durch radioaktiven Zerfall entstandenen Isotopen Pb?°, Pb? und 
Pb?°s, kann aber einen mehr oder minder groBen Anteil an sog. « accessorischen 
Blei » enthalten, dessen Isotopenzusammensetzung einem «gewòhnlichen » ent- 
spricht. In dem hier untersuchten Cotunnit haben wir es offensichtlich mit 
einem neuen «intermediären » Bleitypus zu tun, der durch seinen RaD-Gehalt 
seine recente Bildung beweist. Uberdies zeigt die noch vorhandene Aktivitat 


15,8 


15,7 


15,6 


15,5 


15,4 


176 17,8 


Fig. 1 zeigt einen Ausschnitt des «-8-Diagramms für « gewohnliches Blei ». Ein- 
getragen sind die Isotopenverhältnisse aller bisher gemessenen [Ni 38, Vi 52, Eh 53, 
Gei 54, Eb 54] europäischen Proben von « gewohvlichem Blei », mit Ausnahme solcher, 
deren Isotopenverhältnisse auBerhalb dieses Diagrammausschnittes liegen. Diese Proben 
haben jedoch alle ein « Modellalter » p von mehr als 700-108 a. Die Werte (®) A, B, C 
geben die Isotopenverhaltnisse der drei Cotunnitproben wieder, die Parallelogramme 
zeigen die Fehlergrenzen der massenspektrometrischen Messungen. (Die Form dieser 
Fehlergrenzen ist bedingt durch die kleinere MeBgenauigkeit des Pb). Eingezeichnet 
sind auBerdem die Isochrone fiir ein « Modellalter » von 500-105 a und die rezente 
Isochrone. 
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des RaD, daB dieses Blei noch vor einer mit der Halbwertszeit des RaD ver- 
gleichbaren Zeit mit Uran oder einem seiner langlebigen Folgeprodukte aufs 
engste assoziiert war. 

Andererseits wiirde ein Blei der gefundenen Isotopenzusammensetzung 


-— ware nur diese bekannt — als ein «gewéhnliches Blei » erscheinen, dessen 
Isotopenhaufigkeiten freilich em wenig von denen anderer junger — speziell 
europäischer — Bleisorten abweichen. Auf diesen Punkt soll weiter unten 


noch ausführlich eingegangen werden. 

Infolge seines intermediären Charakters kann also das hier vorliegende Blei 
von zwei Seiten her aufgefaBt werden: einmal als ein durch natiirliche Vor- 
gange separiertes radiogenes Blei, das einen groBen Anteil accessorischen Bleis 
enthalt, zum andern als ein in seiner Isotopenzusammensetzung von den iibli- 
chen europàischen Typen ein wenig abweichendes « gewòhnliches Blei » (vgl. 
Fig. 1). 

Diskutieren wir zunächst die erste Auffassung des Cotunnits als radiogenes 
Blei, so ist die wichtige Frage nach der Geschichte des Bleies in den letzten 
Jahrzehnten vor der Auffindung zu beantworten. Bekannt ist der Zeitpunkt 
der Auffindung und damit auch der des Aufhòrens der Nachentstehung. Als 
obere Grenze fir das Alter der Kristalle kann jederfalls der Zeitpunkt der 
Erstarrung der Lava, bzw. der der Entstehung der Tuffschicht, angenommen 
werden, auf welchen sich die Kristalle gebildet haben. Maf£fgebend für die 
hier im wesentlichen interessierende spezifische Anfangsaktivität, d.h. für das 
anfängliche Verhältnis des RaD zum Gesamtblei ist der Zeitpunkt der Se- 
paration des RaD von seinen langlebigen Muttersubstanzen. Hierbei kommen | 
infrage die Abtrennung des RaD vom Ra, des Ra vom Io und schlieBlich die 
des Io vom Uran. Die Frage, ob und um wieviel sich diese Zeitpunkte von- 
einander unterscheiden, wird durch detaillierte radioaktive Untersuchungen 
der verschiedenen Sublimationsprodukte der Fumarolen mit verschiedenen 
Entstehungsdaten zu entscheiden sein. Es handelt sich hier um einen ex- 
tremen Fall von radioaktivem Nichtgleichgewicht, und die Tatsache, daf in 
den natiirlichen Zerfallsreihen Folgeprodukte sehr verschiedener Halbwerts- 
zeiten existieren, verspricht mittels langwieriger radiochemischer Untersuchun- - 
gen interessante Riickschliisse auf den Mechanismus und den zeitlichen Ablauf 
vulkanischer Vorgänge ziehen zu kònnen. Leider sind seit der Auffindung der 
untersuchten Proben schon einige Jahre verganden, soda das laufende Po?!0- 
Gleichgewicht nahezu erreicht ist und eine Feststellung des Kristallisations- 
datums aus dem Anwachsen des Po?!° heute nur mehr ungenau môglich ist. 

Vergleicht man die in Tabelle II angegebenen spezifischen Aktivitaten mit- 
einander, so zeigt sich immerhin, daB diese — umgerechnet auf den Auffindungs- 
zeitpunkt — innerhalb der Fehlergrenzen gleichen Wert haben. (Hierbei ist 
zu bedenken, daB die Differenz der am weitesten auseinanderliegenden Auf- 
findungszeitpunkte nur etwas mehr als drei Jahre ausmacht, was einem Ab- 
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fall der RaD-Aktivitàt von nur etwa 10% entspricht). Wenngleich eine Nach- 
prifung dieses Ergebnisses an Proben mit genau definierten und môglichst 
weit auseinanderliegenden Bildungszeitpunkten erforderlich ist, so ist dieser 
Befund nicht im Widerspruch mit der Hypothese, daB der Auffindungszeit- 
punkt, der Bildungszeitpunkt und der Zeitpunkt der Abtrennung des RaD 
vom Ra nur um eine Zeitspanne auseinanderliegen, die klein ist im Vergleich 
zu der Differenz der Auffindungszeitpunkte von drei Jahren. 

Das Cotunnit-Blei laBt sich als ein radiogenes Blei auffassen, das aus einem 
« Urankonzentrat » stammt, fiir welches sich nach der RaD-Methode [Hou 51, 
Be 53] ein Maximalalter angeben läfBt. Eine solche Altersbestimmung erlaubt 
jedenfalls nur, einen Maximalwert ftir das Alter des im Vesuvmagma vorhan- 
denen « Urankonzentrats » anzugeben, da nicht bekannt ist, ob der Zeitpunkt 
der Abtrennung des RaD--RaG vom Ra mit der Stôrung des Gleichgewichts 
der Elemente vom Uran bis zum RaD zusammenfällt, und nur die spezifische 
Aktivitàt zum Auffindungszeitpunkt der Altersbestimmung zugrundegelegt 
werden kann. Die Berechnung dieses Maximalalters beruht auf der Tatsache, 
daB bei der Mitnahme des RaD durch den zur Cotunnitbildung führenden 
Sublimationsvorgang mit Sicherheit auch RaG in einem Verhältnis zu RaD 
mitgenommen worden ist, wie es in dem « Urankonzentrat » zum Zeitpunkt 
dieser Abtrennung des Bleis vorhanden war. Die Bestimmung dieses Maximal- 
alters wird naturgemäB sehr ungenau, da im Gegensatz zu allen anderen be- 
kannten Fallen radiogener Bleisorten, ftir welche das Verhältnis des accesso- 
rischen Bleis zum Gesamtblei hôchstens ca. 20-30 % ist, hier durch das Vor- 
handensein accessorischen Bleis in Héhe von iber 99% (*) eine genaue Kennt- 
nis der Isotopenzusammensetzung des accessorischen Bleis erforderlich ist, um 
die zum RaD gehòrige Menge RaG berechnen zu kônnen. Die sonst für ein 
thoriumfreies Uranmineral ibliche Methode, unter den fraglichen Werten fiir 
das Verhaltnis x = (Pb?°*)/(Pb?°) des accessorischen Bleis einen Anhaltspunkt 
aus dem Verhältnis y = (Pb?°)/(Pb?°) zu gewinnen, versagt hier, da der ur- 
sprüngliche ThB-Gehalt der Cotunnite und damit der Th-Gehalt des « Uran- 
konzentrats » nicht bekannt ist. Auch die naheliegende Methode, aus Unter- 
. Schieden der Isotopenzusammensetzung durch verschiedene spezifische Akti- 
vitàt und damit verschiedene Anteile an accessorischem Blei dessen Wert 
von « zu ermitteln, 148t sich hier nicht anwenden, da die bisher gemessenen 
Werte von « innerhalb der Fehlergrenzen übereinstimmen. Jedoch zeigen die 
gefundenen Werte von « und £ = (Pb?°)/(Pb?%) der Cotunnite nur eine so 
geringe Abweichung von denen der bekannten jungen europàischen Bleitypen, 
da8 ohne weiteres ersichtlich ist, daB nur ein geringer Prozentsatz, nämlich 


(*) A. Prurri und D. Miexiacci [Pi 22] fanden Abweichungen des Atomgewichtes 
des Cotunnitbleis von dem gewôhnlicher Bleie und schlossen auf eine Beimengung 
von Uranblei (uraniopiombo) zu gewohulichem Blei in Hühe von 12%. 


aes = 
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hôchstens 4% der gefundenen Menge Pb?°* aus dem « Urankonzentrat » stam- 
mendes RaG sein kann, d.h. da8 das fragliche « Urankonzentrat » ein Alter 
von héchstens 50-105 a haben kann. Da für junge radiogene Bleisorten das 
Verhaltnis AcD/RaG praktisch gleich 4,6% ist, haben wir zur genaueren Be- 
rechnung des Alters ein accessorisches Blei mit den Isotopenverhältnissen 
% = 18,85, Ba = 15,74 verwendet. Diese ergeben sich als Schnittpunkt einer 
durch die «-$-Werte des Cotunnit B gezogenen Geraden mit einer Neigung 
von 4,6% und der rezenten Isochrone des in [Hou 53] gegebenen Isochronen- 
systems (vel. Fig. 1). Damit erhält man für das « Urankonzentrat » im Vesuv 
ein Maximalalter von 17:10°a (Entsprechend ergibt sich aus den Proben 
Cotunnit A bzw. Cotunnit C ein Maximalalter von 17 bzw. 14-10%a). Die 
hier benutzten Isotopenverhältnisse des accessorischen Bleis fiir ein Gesteins- 
blei des Alters p = 0 entsprechen einem Wert von u = (U**),,,,/(Pb204) & 
= 9,7, was durchaus den aus europäischen Bleisorten ermittelten u-Werten 
fiir deren Muttergestein entspricht [Gei 54]. Diese Ubereinstimmung in den 
u-Werten bedingt, daB die Altersberechnung für das « Urankonzentrat » im 
Vesuv keineswegs von der speziellen Wahl des Isochronensystem abhangt, 
sondern daB eine Korrektur mit einem rezenten oder jungen europàischen 
« gewohnlichen Blei » auf ein ähnliches Alter führt. 

Betrachten wir nun das Cotunnit-Blei als ein «gewôühnliches Blei », als 
welches es dem Beobachter z.B. erschiene, wenn man es schon vor einigen 
hundert Jahren gefunden und erst heute untersucht hatte, so daB schon seine 
RaD-Aktivitàt vollstandig abgeklungen ware. Nur nach seiner Isotopenzu- 
sammensetzung beurteilt erscheint das Cotunnit-Blei als ein ausgesprochen 
« anormales » Blei, das demjenigen Typus zuzurechnen ist, für welchen das aus 
dem Isochronensystem [Hou 53, Gei 54] berechnete « Modellalter » geringer 
ist als das geologische Alter. Das « Modellalter » ergibt sich namlich zu 
p = —200-10%a. Das Auftreten solcher anormaler Typen von Bleierzen, wie 
wir sie z.B. in den Vorkommen von Joplin und am Colorado-Plateau und noch 
extremer in einigen der Bleiproben von Sudbury kennen, ist vielfach diskutiert 
worden. Diese etwa als «J-Typen» zu bezeichnenden Vorkommen lieBen es 
einigen Autoren [Je 48, Vi 52] überhaupt als zweifelhaft erscheinen, das er- 
wähnte Bildungsmodell fiir Bleierze für die Berechung des Erdalters zu ver- 
wenden. Es war eingewandt worden, daB zur Berechung des Erdalteis und 
der Isotopenzusammensetzung des « Urbleis» nicht alle analysierten Blei- 
proben verwendet wurden. So ergibt sich z.B. nach diesem Modell für die 
Probe Worthington Mine, Sudbury (Ontario) [Co 52] ein Modellalter von 
p =—4-10°a. Demgegenüber ist die Môglichkeit eròrtert worden [Vi 52, 
Gei 54], daB Bleitypen scheinbar negativen Alters sich ohne weiteres durch 
eine Erhôhung des u-Wertes im Muttergestein oder eine Mitnahme radiogenen 
Bleis bei der Ausscheidung erklàren lassen. Durch die Analyse des heute noch 
RaD-haltigen Cotunnitbleis scheint uns diese Erklarung fur das Zustande- 
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kommen anormaler Isotopenzusammensetzung des J-Typs, der hier gewisser- 
massen im statu nascendi beobachtet wird, auf die erwahnte Art experimentell 
erwiesen, da sich aus dem RaD-Gehalt des Cotunnitbleis ein Wert von ca. 
4 = 110 ergibt, während die Isotopenzusammensetzung zeigt, daB der Wert 
von u während des gròften Teils der Bildungsgeschichte dieses Bleis etwa 
10 gewesen sein mufi. Wir wollen damit keineswegs als bewiesen erachten, 
daB die geologisch-genetischen Bedingungen sonstiger anormaler Bleivorkommen 
des J-Typus eine besondere Analogie zu den hier vorliegenden Bedingungen 
zeigen, doch scheint uns auch fur diese als Erklärung ein Anstieg des u-Wertes 
im Muttergestein oder eine Mitnahme radiogenen Bleis nahezuliegen. Diese 
Auffassung wird dadurch gestützt, daf im Falle der Sudbury-Vorkommen 
auBer der oben erwähnten anormalen Proben [Co 52 Nr. 25 und 26, Ru 54 
Nr. 1, 2, 5, 6, 7,9, 11 und 13] auch drei Proben [Ru 54 Nr. 3,4 und 8] existieren, 
die als durchaus normal gelten künnen, da sich für diese « Modellalter » von 
1800 bis 2000-10 a ergeben. 

Auch im Falle der Proben gewühnlichen Bleis vom Colorado-Plateau [Geo 54] 
scheinen ähnliche Verhältnisse insofern vorzuliegen, als auch hier Isotopen- 
verhältnisse gemessen wurden, die von stark negativen « Modellaltern » bis zu 
solchen reichen (Nr. 6 und 7), fiir die das « Modellalter » mit dem geologischen 
Alter übereinstimmt. Im Falle des Joplin-Vorkommens zeigen alle gemessenen 
Proben negatives « Modellalter ». Hier ist also das zugehòrige Blei des « Normal- 
typs » nicht bekannt, doch fallt auch hier wie in den beiden anderen angefùhrten 
Fallen eine starke Streuung der Isotopenverhältnisse auf. Môglicherweise wird 
sich gerade die starke Dispersion der Isotopenverhältnisse eines Vorkommens 
als Kriterium für «anormale » Bleilager erweisen. 

Zuletzt môgen noch einige Punkte unserer physikalischen MeBergebnisse 
hervorgehoben werden, die vielleicht eine gewisse allgemeinere Bedeutung für 
die Vulkanologie haben kénnten. Der heutige u-Wert der Muttersubstanz, 
aus dem die Cotunnite stammen, ergab sich experimentell zu w= (U?**)/(Pb?°4) > 
>110. Dies ist ein extrem hoher Wert, der weit über dem aller bekannten 
Eruptivgesteine liegt. Dieser Umstand hat uns dazu veranlaft, oben von 
einem « Urankonzentrat » als der Muttersubstanz eines Teils des Cotunnitbleis 
zu sprechen. Zunächst kann hierbei keinesfalls behauptet werden, daf eine 
Erhôhung der absoluten Urankonzentration gegenüber der gewôhnlicher Ge- 
steine vorliegt. Jede Erhéhung des chemischen Verhältnisses Uran/Blei muB 
in dem hier gebrauchten Sinne als Urankonzentration aufgefaBt werden, auch 

dann, wenn diese nur durch einen Bleiverlust zustande gekommen ist. Der 
ÜberschuB an Pb?98 gegenüber « normalen » jungen europàischen Bleisorten 
lässt vermuten, daB eine der Urankonzentration vergleichbare Thorium- 
konzentration stattgefunden hat. Eine exakte Deutung des Mechanismus, der 
die hier gefundene Erhéhung des u-Wertes und damit des Verhältnisses 
Uran/Blei bewirkt, muf einer eingehenden geologisch-petrographischen Ana- 
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lyse der vulkanischen Vorgänge vorbehalten bleiben. Wir môchten im fol- 
genden nur einige Modelle angeben, durch die das Zustandekommen des 
hohen u-Wertes und der gegeniber anderen jungen europàischen Bleierzen 
erhòhten Anteile an Pb?% und Pb?’ im Cotunnitblei zu erklären wire. 

Eine mégliche Deutung ist das Zusammentreffen von Gesteinsblei mit 
weniger als 1% radiogenem Blei aus Mineralanteilen, die Uran und Thorium 
in sehr viel héherem MaBe enthalten als gewéhnliche Eruptivgesteine. In 
diesem Falle bedeutet der oben berechnete Wert von 17:10%a das Maximal- 
alter der uran- und thoriumreichen Minerale. 

Als andere Deutungsmôglichkeit für das Zustandekommen des hohen 
u-Wertes und der gefundenen Isotopenzusammensetzung bietet sich die An- 
nahme, daB der hohe u-Wert durch eine Verarmung des Magmas an Blei be- 
dingt ist. Das starke Auftreten von HCl in den Dämpfen der Fumarolen und 
das Auftreten von Chloriden in deren Sublimationsprodukten ergibt die Méglich- 
keit, daB wegen der leichten Verdampfbarkeit von PbCl, allmàhlich eine Blei- 
verarmung des Magmas eintritt. Eine Nachpriifung dieser Hypothese lieBe 
sich durchführen, wenn eine Abschätzung der bei der Vulkantätigkeit abge- 
gebenen Bleimenge môglich ware. Auch ist von diesem Standpunkt aus das 
Studium der spezifischen Anfangsaktivitàt frischgebildeter Bleiverbindungen 
und anderer Sublimationsprodukte in ihrer zeitlichen Abhängigkeit von der 
Vulkantàtigkeit von bedeutendem Interesse. Die Angabe des Maximalalters 
von 17-106a für den radiogenen Anteil ist in diesem Falle ein Mañ für die Zeit, 
seit der diese Bleiverarmungsvorgange begonnen haben. Eine Erhôhung des 
u-Wertes in dem beobachteten MaBe kénnte freilich auch als Folge von Aus- 
brüchen in jüngster Zeit eingetreten sein, doch würde dieser Mechanismus 
nicht die gemessene Erhôhung des Gehaltes an Pb?% im Blei erklären und 
macht die Annahme einer Dauer der genannten Vorgänge von der angegebenen 
GròBenordnung erforderlich. 

SchlieBlich laBt sich einwenden, daB der gemessene u-Wert aus dem Ver- 
hältnis RaD/Pb2% bestimmt wurde und damit nur das Verhältnis U?**-Aqui- 
valent des RaD/Pb?%, nicht aber das direkt gemessene Verhältnis U?38/Pb?2% 
reprasentiert. Es wire somit denkbar, daB die hohe RaD-Aktivitat lediglich 
dadurch zustande kommt, daB die die Sublimationsprodukte erzeugenden 
Dämpfe Radon aus einer sehr groBen Gesteinsmenge aufgenommen haben. 
So wiirde z.B. eine stark radonhaltige Quelle mit minimalem Bleigehalt gleich- 
falls ein sehr hohes Verhältnis RaD/Pb?° aufweisen, wenn die Zeitdauer 
zwischen der Aufnahme der Emanation und dem Austritt der Quelle vergleich- 
bar mit der Halbwertszeit des Radon ist. Diese Môglichkeit würde im Falle 
des Cotunnitbleis zwar das hohe Verhältnis RaD/Pb?% erklären, nicht aber 
den aus der Isotopenzusammensetzung ermittelten UberschuB an Pb2%. Ein 
weiteres Argument gegen diese Auffassung ist der gemessene ÜberschuB an 
Pb2% im Vergleich zu anderen jungen europäischen Erzbleien, denn wegen 
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der kurzen Halbwertszeit des Thoron ist es besonders unwahrscheinlich, daf 
der erhôhte Gehalt des Cotunnits an Pb? sich allein durch eine Mitnahme von 
Emanation erklaren 1aBt. 


Zum SchluB méchten wir unseren besonderen Dank für die Beschaffung 
der Cotunnite und fiir interessante Auskiinfte tiber diese den Herren Dr. BEVI- 
LACQUA, Prof. IMBo und ganz besonders Prof. PARASCANDOLA, Neapel, der die 
Proben gesammelt hat, aussprechen. 

Fur anregende Diskussionen danken wir den Herren Dr. HùGI, Bern, und 
Dr. PicctoTro, Brüssel. 

Der Schweizerischen Kommission fiir Atomenergie danken wir fir die 
Bereitstellung von Mitteln fiir die massenspektrometrischen und radioaktiven 
Messungen. 
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RIASSUNTO (*) 


\ 


Il piombo di tre campioni di Cotunnite del Vesuvio è stato esaminato per la sua 
radioattività e allo spettrometro di massa. La radioattività fu identificata come dipen- 
dente dal RaD e dai suoi prodotti di decadimento e risultò dopo estrapolazione al tempo 
di ritrovamento 21,2+0,5 decadimenti/s-mg Pb. Il Piombo della Cotunnite è da 
ritenere un tipo «intermedio », finora sconosciuto, consistente per circa il 0,7% in 
piombo di decadimento radiogeno e per il resto in piombo «ordinario » recente. Il 
piombo di decadimento ha origine da un « concentrato di uranio » la cui antichità mas- 
sima, col procedimento del metodo RaD, fu potuta fissare a 17-105 a. La composi- 
zione isotopica (v. Tab. I) rivela una anomalia simile a quella dei giacimenti piombi- 
feri di Joplin, Sudbury, ecc.. Si propone un modello per il meccanismo di formazione 


di questi « Tipi J ». 


(*) Traduzione a cura della Redazione, 
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A Non-linear Generalization of the Schrôdinger 
and Dirac Equations. 


M. ScHONBERG 


Faculdade de Filosofia, Ciencias e Letras da Universidade de Sao Paulo (Brasil) 


(ricevuto il 23 Aprile 1954) 


Summary. — It is shown that it is possible to modify the Hamilton-Jacobi 
equation by the introduction of generalized potentials that lead to a 
modification of the total momentum, without introducing new forces. 
By means of the solutions of the modified Hamilton-Jacobi equations 
it is possible to describe the motions of more general classical ensembles 
of particles than those corresponding to the solutions of the ordinary 
Hamilton-Jacobi equation. The extension to quantum mechanics is 
straightforward and leads to a generalization of the Schrodinger equation 
by means of a system of non-linear equations. The same procedure leads 
to a generalization of the Dirac equation. 


1. - The Schrédinger equation for a spinless particle of mass m and charge e 


moving in an electromagnetic field described by the potentials A,, A 


(1) 


ih 


Da al = 24) B+ eA, 


ct 2m\i dx ce 


corresponds to the classical Hamilton-Jacobi equation: 


(2) 


OS 1 /08 € È 
Sii epr, E ed 0. 
aa Er c | ae 


Each solution of the Hamilton-Jacobi equation (2) defines a set of co® traject- 
ories of the particle determined by the first order differential equations: 


(3) 


p=mv + “ = = £(v = velocity of the particle). 


-% 
di 
| 


A NON-LINEAR GENERALIZATION OF THE/SCHRÒDINGER AND DIRAC EQUATIONS 675 


This association can be rendered more pictorial, by considering an ensemble | 


of identical particles distributed in space with a density (of points) o(f, x) 
satisfying the partial differential equation: 


(4) 2 + div (gv) = 2° | div £ (Se ca)}-c. 
Thus we get a kind of pressureless fluid, whose Euler equations are the Newton 
equations of motion of the particle and whose equation of continuity is (4). 
Each quantal state corresponds to a state of our pressureless fluid. 

The motions of the pressureless fluid associated to the solutions $ of the 
Hamilton-Jacobi equation are not the most general ones. This is seen more 
clearly when there is no magnetic field, since it follows from equation (3) that, 


in this case, the motion of the fluid is irrotational: rot v = 0. When there ~ 


is a magnetic field, there it still a restriction on the vorticity: 


(5) | rot» = — _ H (H = magnetic field). 


In order to get the general motion of the fluid, we shall introduce Clebsch 
parameters : 


e 
(6) mv + A i 


It is easily seen that the Hamilton-Jacobi equation (2) can now be replaced 
by the system of three equations: 


oS ESTE RCA 
À LG, = 0 
(7) ot es Che Saal Ga "a LF ee ; 
2 0A 0A du du 
et NS E ULI eg eed ES I 
(8) a o a 


By differentiating both sides of (6) with respect to #, we get, as a consequence. 


of (7) and (8), the equations of motion of the particle: 


dv v x 
EE = \H E — electric field). 
(9) | map =e(E+ nu) Re Soe 


Equation (7) may be regarded as the Hamilton-Jacobi equation for a part- 
icle moving under the action of the electromagnetic field and an extra force f 
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derived from the generalized potentials (04/0) and — A(du/ ex): 
ef, Ou Of, Ou (eu 
Ra | es 
Ft f dx (2 ca) SC]; (4 =) LS iz 2) 


In such a condition, the total momentum is given by the following equation 


(11) p=mv PA 


which replaces (3). Equation (11) is obviously equivalent to (6). It is easily 
seen that 


__ Oufoà. OA 0A (Ou op 
oe LE a By ale ele 3x)? 


so that the force f will vanish when 2 and wu satisfy equations (8). Thus we 
get an explanation of (9). Equation (4) is now replaced by the following 
equation: 


00 Si ON EOS EC ci 
CIS) % divi (2 : 4 Ax =i), 


It is interesting to notice that equations (7), (8) and (13) can be derived from 
a variational principle: 


A 2 Ù 
(14) felt (e SA +42) + ea, dtd = 0. 


Ox C 


2. — We shall now discuss a generalization of the Schrédinger equation 
which corresponds to the above generalization of the Hamilton-Jacobi equation. 
The Schrédinger equation (1) can be derived from the well known variational 
principle: 


A 2 
(15) o (finita Page Vars —0 
ot 2m 1 Ox C 


We shall now consider the Schrédinger variational principle for a particle 
which is under the action of an electromagnetic field and of a force derived 
from the generalized potential 7(0u/0t) and — A(du/ex): 


(16) offer: d A: 12) P Cr ! Ae) ea | 


2m è Ox Ox 


o” 3 - x "i - ” È a - i 
ae 2% mn ; 1 ] - i aA 
Fy, È ao 2 
È | 
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By variating only ¥* and Ÿ we get the equations: 


(17) in OF _ DI (È 0 e 
2m 


n du\? Ou 
AR oe 22) P+( + a2 ly, 


ot 
; SES 1/70, Ou\? © 
17 Sue bes SE 
(17a) ih ie (; Pa + — -A— A =) pri (ea, + À à jes ; 


By variating 4 and u, we get from (16) two more equations: 


Ou Ou _ o) : . 
(13) On, Soot Le ct a We) + div (4j) = 
o=P*P, 
OE (Aloe dre al 


It follows from (17) and (17a) that: 
00 AE 
(20) e ad AE 


By means of (20), we can transform the equations (18), as follows: 


SIAT GÀ ° Mi 
(21) orata NE 5; (ue) + div (uj) = 


The first equation (18) and the first equation (21) correspond to the two classical 
equations (8). The analogy with the classical theory of section 1 becomes 
even clearer by putting in evidence the amplitude R and the phase S/h of W: 


i 
(22) P = Rexp PA] : 
We have 
: 17/08 ee , du 
@) den, je, = 5 (Epa), 


ESE du\? h? AR 
en re + = 0 
we athe nl = +14) ne 2m E : 
5 de i ay) e fes 2 4 4 se) =o 
>) A tar) (E Hes | Lie À 


LP STR PI NO, PAT e 
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The equations for À and w can be written as follows: 


oA” -\ en Ou Ou 
eA eee. slo nr 
ce) pra CRA 


The only difference between the above equations and those of section 1 arises 
from the term —(#?/2m)(AR/R) in the left hand side of (24). This term 
vanishes in the classical limit. 

. It is interesting to notice that the second equation (18) and the second 
equation (21) express the conservation of quantities whose spatial densities 
are Ao and uo, the current densities being Aj and wj, respectively. A similar 
circumstance appears in the classical theory, because it is possible to replace 
the equations (8) by others similar to the second equation (18) and the second 
equation (21). Moreover any function F of À and y satisfies the equation 
(0/ét) (Fo) + div (Fj) == 0, so that: 


iY 


+a 


(27) a F(A, ujo(t, x) dx = 0. 


— © 


It follows from equations (19) that @ and j are multiplied by a factor |a |? 
when ¥ is replaced by a¥, a being an arbitrary constant. Therefore v remains 
constant and the equations for À and w are not changed. We may say that 
the sets Y, A, w and aŸ, A, u describe the same-condition of the particle. 


WY may be normalized as usual and |7|*dx can be taken as the probability 


of finding a charge e in the element of volume dx at the point in consideration. 
This choice is actually required by the form of the interaction terms with the 
electromagnetic field. It corresponds also to the meaning of the classical 
density o of section 1. The physical interpretation of 4 and w is still unknown. 


3. — We shall now see that it is possible to find solutions of (17) and (18) 
by taking A= G(4), G being an arbitrary function. Indeed the function Y% 


(28) | v= exp |; | away), 


Ho 


satisfies the usual Schrédinger equation (1). Since o and j can be expressed 


= 


_ in terms of ¥Y’ and ¥’’* only, the determination of w depends on the solution 


of a linear partial differential equation of the first order, whose coefficients 
are known, once (1) is solved. The generalization of the Schrédinger equation 
corresponding to À — G(u) appears as trivial, being essentially equivalent to 
a gauge transformation of the potentials Ay, A. 


iM 
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We shall now prove that 4 and w are independent when n 


m 0x’ 6x’ 


(29) | n= rot (» sie A) = 1 04% 


UD, 


is not 0 everywhere in space, and only in this case. Indeed, it follows from 
the form of the third member of (29) that À is a function G of w and t, when 
n=0. By taking into account that 4 is a solution of the first equation (21), 
we get: 


Of Ode OG” OG fou du\ . 36 
30 meet na 
(40) ot uc Ox 04 du (+0 v= 


It is easily seen that: 


© 
(31) DM _ rot (v An). 
ot 
Therefore, a necessary and sufficient condition for n to vanish everywhere for 


all ¢ is that y(t, x) = 0 for all x, at some instant of tone ie 
It follows from (23), (17) and (18) that: 


© 0 lat E GI AR 
Sleek HIER à 
(22) mv Ox Je 1 > gli | Sin da È 
This equation corresponds to the classical equation (9). It has the same form 
both for the Schrédinger equation and for the above generalization of the 
Schrodinger equation. By considering the solutions of the differential equation 


(33) | a = v(t, x), 


we get a family of mappings of the space into itself, in which the points x(t) 
correspond, for different values of #, to the point x, = x(0). This flow con- 


serves the measure li o(t, x)dx of the measurable sets Q(t), as a consequence — 


24) 
of (20). It has also some other properties, that are very similar to these of ‘the 


motion of an inviscid barotropic charged fluid, whose charge and mass den- 
sities have everywhere the ratio e/m. These properties can be proven by 
modifying very slightly the proof of the corresponding hydr Signal: theorems, 
Thus we get the Cauchy equation 


\ 0, x) © 
CA nis} rai fa, 


o(t, x) 


a ee ER) SR TA fe Me a oe D x LA yea tk pe er, De PP ER ey oe PRE Lote POE Re, ON à i oe Le 1 de 77 Lie 
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which expresses the vector y(t, x) in terms of the vector (0, x,), x, being the 
point which is mapped on x at the time ¢. The well known Helmholtz vortex 
theorems have their analogues in our flow. The vortex lines are here replaced 
by n-lines, tangert at all their points to the respective y vectors at an instant 
of time, and the vortex tubes by n tubes. 

The equations (17) and (18) may be replaced by (32) and (00/0) + 
+ div (ov) =0.. In order to get the ordinary Schrodinger equation we must 
impose the supplementary condition n = 0. This way of writing the equat- 
ions of our generalization of the Schrodinger theory has the advantage of 
showing that the above treatment of the Clebsch parameters is not the most 
general one, since from those equations and the third equation (23) it does 
not follow that Ÿ, 4 and w satisfy necessarily equations (17) and (18). It can 
be shown that the most general choice of the Clebsch parameters leads to the 
following equations 


nov i(k ae du? du n 
0 lo 
CRIS e | 4 OK 
(36) ot Ox mM\0x € dx Ou 
du 1/08 e du eK 
dui Ae a =F 


K(t, À, u) being an arbitrary function of t, À and w. The equations (35) and 
(36) can be derived from the following variational principle: 


(37) o [| see ur a+ E a) 


2m i Ox 


—- (ea, + 2 = —- K(t, A, 7) var | dida =05 


4. — The Dirac equation can be generalized by replacing the Dirac variat- 
ional principle by the following principle: 


ei ae CO du 
38 lap nets “élue 
(38) 0 | wy = Won = A+ A) 


— W*Bme?Y — P* (ea. + 12) | dtdx = 0. 
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The equations for Y, 2 and uw are now: 


Oe ORO: e cu 
(39) Uda (; da lA +e + per? + (cd, +af)v_o, 
oA oA Cli Role 
4 CR ES ale. perg 
GO DE EE O ordigno 
(41) A Hdi LES TC ae, 


We have now also the three conservation equations: 


Ô 
a; (Me) + div (uj) = 


00 + 
9 ce ke 
(42) > +divj=0, > ° (10) + div (Aj) = > 


It is possible to generalize, in à different way, the Dirac equation by means | 
of the following variational principle 


OP hire e h 0D  0D* 
LETO A EEE pgs ee ce hee == RD ea ee = 
(43) [te ih À P*ca (; 3 At ilo Ae = pol) 
* ” 
— me —¥* (ed, + — ion: do Li Ao} )v}aax 0, 


in which both Y and ® are four component spinors. The equations for Y 
and ® are: 


OP A ari opt 


or am OOTY 


© © 
—+j—|D=0. 
(45) (e ce) | d 
We have now the following conservation equations 
00 ere 0 ‘ È 
(46) Dr +divj=0, si (D.,0) + div (D.j)=0, 


the ®, being the components of ®. 

The generalization of the Dirac equation corresponding to the variational 
principle (43) seems to be of a different nature than those involving the two 
scalars À and uw. It is however possible to build with 4 and y a kind of wave 
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dA (47) O=Viexp 


es 
when 4 is taken dimensionless. It is easily seen that: 


* 
NC (or 20 00 0). 


i iy Ot ot 
: Se où K(p,0O 20% 
Lae * 
| 1 0x 2i ie, da dx e) ; 


| Given a solution y of the ordinary Dirac equation, we can form a solution 
_of the system (39)-(40) by taking 


(49) : Zip; VASO 


and a solution of the system (44)-(45) by taking: 


(50) + W—y, D=0. 


Thus it is seen that the generalized states described by Ÿ, 2, uw or by ¥, D | 
include the ordinary states of the Dirac particle as particular cases. 

A more detailed discussion of the generalized: wave equations will be pub- 
lished later in this journal. | 


M po DREAM S SION ON) 

fe) Si dimostra che è possibile modificare l’equazione di Hamilton-Jacobi introducendo 
dei potenziali generalizzati che conducono-a-una modificazione dell’impulso totale senza 
richiedere l’introduzione di nuove forze. Per mezzo delle soluzioni delle equazioni di 
Be Hamilton-Jacobi modificate è possibile descrivere i moti di un numero di insiemi gene- 
3 rici di particelle classiche maggiore di quello corrispondente alle soluzioni dell’equa- 
«> °°». zione ordinaria di Hamilton-Jacobi. L’estensione alla meccanica quantistica è imme- 
diata e porta a una generalizzazione dell'equazione di Schrodinger per mezzo di un 
sistema di equazioni non lineari. Lo stesso procedimento conduce a una generalizza- 
dr zione dell'equazione di Dirac. 


TLT 
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(*) Traduzione a cura della Redazione. 


Energetic Nuclear Collisions in the Upper Atmosphere. 
II. The Muon Component. 


U. HABER-SCHAIM 


Physikalisches Institut der Universität - Bern, Schweiz 


G. YEKUTIELI 


The Weizmann Institute of Science. - Rehovoth, Israel 


(ricevuto il 29 Aprile 1954) 


Summary. — The flux of pions in the upper atmosphere calculated in 
Part I is used to determine the flux of muons in the atmosphere. In 
particular, the energy spectrum of muons at sea level in the energy band 
30-1000 Mc? is evaluated and compared with experiment. 


In Part I of this paper (+) the flux of charged pions in the upper atmosphere 
was calculated starting from a given spectrum of primary protons and «-part- 
icles. Here the flux of the muons created by the decay of pions is calculated 
down to sea level. pane a 

The transport equation for muons of energy between E and E+d£ at 
atmospheric pressure p is 


au* (EH, p) 1.1u*(E, p) | 1267*(1.23, p) 


(1) » ip 


Ep 


Here 1.1/Ep and 126/Ep are the mean free paths for decay of muons and 
pions respectively (E is expressed in nuclear mass units). Energy losses of 
muons due to ionisation in the atmosphere and the energy spread due to their 
angular distribution in the rest system of the pions are neglected. 


(1) U. Haser-ScHarm and G. YegKutIELI: Nuovo Cimento, 11, 172 (1954). 
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The solution of (1) is given by 


(2) 


126 


uX(B, p) = = pre gusts a*(1.23 B, t)dt . 


0 


Following part I [eq. (13), (20)] we divide the total flux of pions x* in two 
parts x and x initiated by primary protons and «-particles respectively: 


(3) 


a*(E, p) = 


n(E, p) + (E, p) = 


= exp p/d] À (p/2)kas(B) + exp [— pla Ÿ (/A)*%s(B) 


where À and A, are the collision mean free paths for nucleons and «-particles. 


(5) 


Dividing the muon flux in a similar fashion: 


ee) = 2 pl, p) > HE, p) ; 
pi =1 
JA 
1267,(1.23 E) 2 
RARE RUE (DA) | Oe ee AT, 
0 
ola, 


E 2 


EE, 2) = (ay 126%, (1.232) Ore I di pig os 


(0 


We approximate u* by 


. (4a) 


2 4 


u*(E, p) = 2 HE, D) + 2 HE, PD) + (È, p). 


TABLE I. 

È E R kg/em? pb M* Mae 2.1-M* 
30 9.4 1.8-1075 2.5-10-4 4.4-10-4 5.2-1074 

50 16 3.9-10-8 8.8-1075 1.6-1074 1.8-10-4 

100 31 4.9-10-6 1.8-1075 3.8-1075 3.8-10-5 
200 59 3.3-1078 3.2-1078 8.6-1075 6.7-10-8 
500 117 1.6-10-9 321110 1.2-10-8 6.7-10-? 
1000 180 1.3-10-8 4,2-10-8 2.0-107? 8.8-10-8 


| eee à ER ele es Ot ee oe ee ~ ay eee O° Le +? J cS IR I I e ayo 
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By a calculation similar to that used in the appendicies of Part I it can be 
shown that 7, (H > 30, p) is less than 2.5% of the total flux u*. Table I 
summarises our results for the differential and integral fluxes at sea level. 

The range in column two is taken from GEORGE (2) and the experimental 
values in column five from BARRET et al. (*). 

Our calculated results given in columns three and four were obtained from 
a primary differential spectrum of protons aH-* with a = 0.57 and s = 2.56 (4). 
Evidently the calculated values are too small. 

The factor a enters throughout the calculation as a proportionality factor 
only. We think our model yields more reliable results for muons of energy 
about E=100; therefore by normalizing the calculated flux to the observed 
flux at that energy we find a = 1.2, i.e. an increase by a factor of 2.1 as 
compared to the old values. The corrected results given in column 6 are 
higher than the observed ones for E< 100 and lower for E > 100. This can 
be further corrected by a slight decrease in the exponent s. However the 
comparison between theory and experiment is strongly dependent on the range 
energy relation and since this relation has not yet been verified at our energies 
we shall not make any attempt to improve the agreement with experiment 
by charging the exponent of the primary spectrum. 

The corresponding integral spectrum of protons is 0.77-E-15 and has a 
value of 1.66-10-$ at 4800 Mc? as compared to 4-10 given by KAPLON et al. (5). 

In conclusion it seems that the nucleon-pion cascade described in Part I 
together with the above primary proton spectrum yields a muon spetcrum 
comparable to the one obtained from experiments underground. 


(2) E. P. GroRGE: Progress in Cosmic Ray Physics (1952), p. 413. 

(3) P. H. Barret, L. M. BoLLINGER, G. Cocconi, Y. EISENBERG and K. GREISEN: 
Rev. Mod. Phys., 24, 133 (1952). 

(4) This primary spectrum is based on older data of the muon flux underground 
which are considerably lower than those of BARRET et al. Compare U. HABER-SCHAIM: 
Phys. Rev., 84, 1199 (1951). 

(5) M. F. KapLoNn, D. M. Rirson and E. P. Wooprurr: Phys. Rev., 85, 933 (1952). 


RIASSUNTO (*) 


Il flusso dei pioni nell’atmosfera superiore calcolato nella Parte I si utilizza per 
determinare il flusso di muoni nell’atmosfera. In particolare, si calcola e si confronta 
con i dati sperimentali lo spettro d’energia dei muoni al livello del mare nella banda 


d’energia 30-1000 Mc?. 


(*) Traduzione a cura della Redazione. 


LETTERE ALLA REDAZIONE. 


(La responsabilità scientifica degli scritti inseriti in questa rubrica è completamente lusciata 


dalla Direzione del periodico ai singoli autori) 


Bilancio energetico della radiazione cosmica presso l’Equatore. 


V. BENZI 


Istituto di Fisica dell’Università - Cagliari 


(ricevuto il 15 Aprile 1954) 


Sulla base di precedenti bilanci per 
latitudini settentrionali (4:2) si è effet- 
tuato un bilancio energetico per il flusso 
verticale della radiazione cosmica vicino 
all'equatore (A=3° N), tenendo conto, al- 
meno in parte, della produzione di me- 
soni K oltre a quella dei mesoni 7 e w. 

L'intensità verticale della componente 
mesonica di quest’ultimo tipo è stata sti- 
mata considerando separatamente il 
gruppo dei mesoni veloci (pe = 7-10? 
MeV) da quello dei mesoni lenti. 

L'intensità del primo gruppo in fun- 
zione dell’altezza è stata ottenuta sot- 
traendo dalla curva della componente 
dura, stimata da Rossi (!) per l’inter- 
vallo 1030<x<250 g/em?, e da BALA- 
SUBRAHMANYAN et al. (*) per i rima- 
nenti 250 g/em?, una curva del tipo 
exp [— x/125] normalizzata in modo da 
dare l’intensità della componente dura 
per x = 0. 

Per la componente dei mesoni lenti 
si sono usati i risultati di CONVERSI (4) 
estrapolati nella parte più alta dell’atmo- 


(1) B. Rossi: Rev. Mod. Phys., 20, 537 (1948). 

@ G. Puppi: Suppl. Nuovo Cimento, 10, 115 
(1953). 

(8) V. K. BALASUBRAHMANYAN, G. S. GOK- 
HALE, A. W. PEREIRA e A. $. Rao: Rep. Bombay 
Conf. (1951), p. 75. 

(4) M. CONVERSI: Phys. Rev., 79, 749 (1950). 


Ch ossi 


sfera secondo il criterio indicato da 
Rossi (1). 

Applicando una correzione del 15% 
per tenere conto della divergenza ango- 
lare dei mesoni rispetto alla verticale (?) 
si ha che per la produzione dei mesoni {1 
da parte della radiazione primaria inci- 
dente verticalmente a latitudine 3° N 
vengono dissipati 


0,213 GeV cm ?s ?sterad-1. 


Per valutare la perdita di energia 
sotto forma di neutrini, sono stati con- 
siderati 1 seguenti casi: 


a) Gli unici processi di produzione 
per i neutrini sono i decadimentiz— w+v 
e u>e- dy. 

b) Oltre ai precedenti ed in compe- 
tizione con essi si deve considerare il 
decadimento K4u +Y +. 


c) Caso identico a b) ma con il 
decadimento K —>w+2v. 


I casi b) e c) si sono trattati assu- 
mendo che 1/3 dell’energia apportata dai 
nucleoni primari di energia £ > 50 GeV 
vada nella produzione dei. K, e che 
nel decadimento di questi mesoni 1/3 
dell'energia disponibile vada attribuita 
ai mesoni uw. 

I risultati sono riportati in tabella, 
insieme con i valori della potenza dissi- 


Sur 
he E 


Pe SENN ca Ur 


hive, t 


pata in ionizzazione da tutte le compo- 
nenti secondo NEHER (5), e con le per- 
dite sottoterra e nelle evaporazioni nu- 
cleari. La stima delle perdite in quest’ul- 
timo caso, fatta in base ai risultati otte- 
nuti per latitudini settentrionali (2) e al- 
l’effetto di latitudine per la componente 
neutronica (6), è piuttosto incerta. 


(a) (b) (c) 
| 
Tonizzazione | 0,400 | 0,400 | 0,400 
Sottoterra 02037 +) 705037 0,037 
| Ev. nucleari | 0,040 | 0,040 | 0,040 | 
| Neutrini 0,200 | 0,260 | 0,330 | 
Totale 0,677 072% 0,807 


GeV cm? s-1 sterad 1 


Per il confronto con la potenza d’in- 
gresso si sono considerate separatamente 
la componente protonica primaria e la 
componente composta di nuclei più pe- 
santi (in maggior parte He). Mentre per 
quest’ultima si sono uasti gli spettri 
integrali indicati da KAPLON et al. (7), per 
la prima si sono impiegati i risultati di 
VAN ALLEN e SINGER (8) per l’intervallo 


- di energia E <14 GeV, dopo aver sot- 


tratto il contributo della componente 
pesante ed il fondo pseudoprimario do- 


(6) H. V. NEHER: Prog. in C.R. Phys., 1, 
cap. V (1951). 

(6) J. A. SIMPSON e W. C. Facor: Phys. 
Rev., 90, 1068 (1953). 

(?) M. F. KAPLON, D. M. Ritson, B. PETERS 
e H. L. REYNOLDS: Phys. Rev., 85, 295 (1952) 

(8) J. A. VAN ALLEN e S. F. SINGER: Phys. 
Rev., 78, 819 (1950). 
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vuto all’albedo, secondo la stima fat- 
tane da TREIMAN (9). Si è poi usato uno 
spettro di potenza del tipo #6 per 
H>4500 GeV, normalizzato in base ai 
risultati ottenuti con le lastre (19), mentre 
per le energie comprese tra le precedenti 
si è usato uno spettro del tipo 2-14 op- 
portunamente normalizzato. In tal modo 
con un taglio di latitudine di 13,8 GeV 
si ottiene una potenza d’ingresso di 


1,00 GeV em? s71 sterad* . 


Le discrepanze con i valori riportati 
nelle colonne a), d) e c) rispettivamente 
risultano pertanto del 33%, 25% e 20%. 
Si noti che con l’impiego di un unico 
spettro di potenza del tipo Æ-11 si ot- 
terrebbe una potenza d’ingresso circa 
tre volte maggior di quella dissipata 
nell’atmosfera. 

Riguardo alle discrepanze riportate, 
quella relativa allo schema c) può attri- 
buirsi senz'altro all’errore sperimentale, 
che è stimato dell’ordine del 20%, mentre 
quella relativa allo schema a) potrebbe 
forse suggerire che tale schema è inade- 
guato ad una corretta descrizione dei 
processi di interazione della radiazione 
cosmica con l'atmosfera, quando si con- 
siderino latitudini in prossimità del- 
l’equatore. 


Ringrazio vivamente il prof. B. FER- 
RETTI per avermi suggerito l’argomento 
trattato e per i numerosi consigli, ed il 
prof. G. FRONGIA per il suo costante inte- 
ressamento. 


(®) S. B. TREIMAN: Phys. Rev., 91, 957 
(1953). 

(1°) M. F. Kapton, D. M. Rrrson e E. P. 
WooDRUFF: Phys. Rev., 85, 933 (1952). 
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Un metodo per i problemi d’urto 
con un potenziale Coulombiano medificato. 


D. BARONCINI 


Istituto di Fisica dell’ Università - Bologna | 


(ricevuto il 19 Aprile 1954) 


| 
Dovendo trattare con il metodo delle onde parziali l’urto tra due particelle | 
cariche, che interagiscano per effetto di forze nucleari e di forze Coulombiane, può 
convenire di eliminare il potenziale Coulombiano. A partire dall’equazione radiale, 
se ne sviluppa la soluzione in serie delle soluzioni Coulombiane radiali; si elimina | 
"I in tal modo il potenziale Coulombiano e si trasforma l’equazione radiale in un’equa- | 
zione integrale, che contiene come kernel il trasformato del solo potenziale nucleare | 
secondo le funzioni Coulombiane, e come incognita il coefficiente dello sviluppo in 
serie. L’equazione cosi ottenuta descrive lo scattering da parte del potenziale nu- 
cleare trasformato e, perciò, la sua soluzione ha il convenzionale comportamento 
asintotico di un’onda piana sovrapposta a un’onda diffusa. La fase di quest’ultima 
è data dalla parte singolare del termine che rappresenta l’onda diffusa. 

| Esprimendo la soluzione dell’equazione radiale di partenza mediante i coeffi- 
‘cienti così ricavati, lo sviluppo in serie delle funzioni Coulombiane risulta nella 
sovrapposizione di una funzione Coulombiana regolare e di un’onda Coulombiana 
a diffusa, sfasata rispetto alla prima. Si verifica che la fase che compare nella seconda 

é la fase stessa che caratterizza la soluzione dell’equazione integrale. 


1. — Spazio delle coordinate, 


Sia V(r) il potenziale, che modifica il potenziale Coulombiano. L’equazione 
radiale di Schrodinger é: 


ae, 24k, |, +1) nk 
(1) A k i 2 + U(r)| ki = 0, 
dove 
k? = 2uE/h?, U(r) = 2uV(rÿ#?, n=efto. 


Sia L,(kr) la soluzione radiale Coulombiana, regolare in r=0; Z,(kr) risolve 


van 
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l'equazione radiale 


Cl dr Ul + 1) nk 
I II 2 
AO Ê 2 710 
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e soddisfa la condizione asintotica: 
LD, (kr) = (kr)~? sin (ar — 2 +n,—n log 2) ; 
Se si sviluppa la soluzione della (1) in serie delle funzioni L,(kr) ponendo 
(2) R,(r) = [uz die 
la (1) si Ta nella seguente equazione integrale: 
P (k? — k'?)u,(k') + faire ir, k') = 0. 


La (3) è l’equazione integrale per scattering dovuto a un potenziale il cui kernel 
è ottenuto trasformando U(r) mediante le soluzioni Coulombiane: 


(3) W,(&', k") me dr Ly (k') U(r) L,(k"r) . 


Indicando con 
(4) u(k') = d(k —k') + xi(k, k') 


la soluzione di (3), la (2) esprime le soluzioni della (1) nella forma 


(5) Ri(r) = L(kr) + | 4%, k)Li(k'r)dk", / 


dove il secondo termine, che rappresenta la parte diffusa dell’onda Coulombiana, 
asintoticamente deve potersi esprimere nella forma 


Il l 
(6) tg Or cos (i — = + m —n log 2 - 


Nella (6) (1/kr) cos (kr — la/2+ m —n log 2kr)) è Vandamento asintotico della 
soluzione Coulombiana irregolare, K,(kr), e Ô, rappresenta lo sfasamento dell’onda 
Coulombiana diffusa rispetto all’onda Coulombiana regolare L,(kr). 

Il fattore tg 6,; che compare nella (6), si ricava direttamente dalla (4). 

Infatti solo la parte singolare della funzione 7,(k, k’) dà contributo = 0 all’in- 
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= + 


tegrale a secondo membro di (5): questa parte singolare è, d’altra parte, il termine 
che definisce lo sfasamento della soluzione (4) della (3) rispetto all’onda piana 6(k—k’). 
Se si indica con 6, tale fase, si ha: 


(7) tgò, = — a lim (k — k') x,(k, k’) 


kik 


e si può scrivere la parte singolare della funzione y,(k, k') nella forma 


a s i 
PATES 
da cui si ha: 


1 


k—k!'° 


1 ; 
EC k')L\(k'r) dk'= — — tg di | run) 
| 3 


Poichè la soluzione regolare D,(k’r) ha la forma 
exp [ik'r]f,(k’r) + exp [— ik'r]fa(k'r) 


. e la soluzione irregolare K,(k'r) 


} 


i[exp [tk'r]f,(k'r) — exp [— th’r]f,(k'r) ] 
si verifica che per grandi valori di r 


! 


Lt: = tg OK (kr) 


1 ; 
—- tg di | Lak 
It 


e quindi che le due fasi coincidono: 6, = 6, 

La soluzione del problema di scattering senza potenziale Coulombiano descritta 
dalla (3) dà direttamente la fase dell’onda Coulombiana diffusa, come richiesto 
dalla (6). 

2. — Spazio dei momenti. 

Sia 


(8) (kK? — p*)pi(p) + [aa [K,(q. p) + K.(p, gle.(q)= 0, 


l'equazione radiale nello spazio dei momenti e siano K,(p, q) e K.(p, q) rispettiva» 
mente i kernel ottenuti trasformando il potenziale V(r) e il potenziale Coulombiano 
allo spazio dei momenti. La trasformata £;(k, p) della soluzione radiale Coulombiana 


perni & 


ory. risolve l'equazione. STE 


7 . \ ~ i 7 È { 3 5 = 
sai e an pa) (ko py / q?dgK .(p, q)E.(k, q) = 0. eo. 
| Se si pone n i à 
Ore gp) = fe Hi, pal) aR’, ESTA ga 
si ha da (8) l’equazione integrale pr 
ANI Ge R')u(k') + | ul) (ke, ") AR" = 0 
che coincide con la (3). La soluzione di (8) è quindi Re: È 
* (11) i vi(p) = E. (k, p) + fa k')E,(k', p) ak’ , . 
se È, A N PR, 
5 i Le 


dove il secondo termine rappresenta la parte diffusa dell’onda Coulombiana, la « ui 
fase, rispetto alla soluzione Tesoro 0: é data direttamente dalla soluzione (4) della (10), 
ossia della (3). È \ 
È In una stesura più accurata, si mostrerà l'applicazione del procedimento al 
| problema dell’urto mesone-nucleone (7=3/2) a basse energie. 
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Some Measurements on Overheated Liquids. 


L. BERTANZA, G. MARTELLI and A. ZACUTTI 


Istituto di Fisica dell’Università - Pisa 


Istituto Nazionale di Fisica Nucleare - 


Sezione aggregata di Pisa 


(ricevuto il 26 Aprile 1954) 


It has been recently shown by GLA- 
SER (1) that overheated liquids in part- 
icular thermodynamical conditions may 
be used as detectors of ionizing particles. 
The device based on this principle has 
been called by GLASER a « bubble cham- 
ber »: with it he has been able to pho- 
tograph some tracks of particles cross- 
ing the liquid. 

The great importance of the deve- 
lopment of such a technique in relation 
with the high stopping power of the 
liquid is self-evident. 

On the other hand there are good 
theoretical reasons for not expecting a 
working bubble chamber with a sensitive 
volume larger than some tenths or some 
hundreds of cubic centimeters, depending 
upon the nature of the liquid and the 
operation conditions. 

This fact might be a very ‘strong 
limitation to the field of applicability of 
the bubble chamber, because of the diffi- 
culty of overcoming the ‘competition 
between scattering and magnetic cur- 
vature. 

Since the bubble chamber of GLASER 


(4) D. ‘A. GLASER; 
(1952); 91, 762 (1953), 


Phys. Rev., 87, 665 


was rather small (only few cubic centi- 
meters) we built a bubble chamber having 
a sensitive volume of about 60 cm3, in 
order to investigate what are the possi- 
bilities of increasing the size of a bubble 
chamber, without reducing the efficiency 
of the device itself, and to enable us to 
carry on systematic measurements. 

A detailed description of the device, 
and some photographs of tracks, will be 
given in Il Nuovo Cimento. 

This bubble chamber has been filled 
under vacuum with spectroscopically pure 
pentane, and the measurements have 
been carried out at temperature of 
~ 138 °C, as suggested by a theoretical 
calculation (2). 

A pressure of about 20 atmospheres 
was applied to the liquid. In the pre- 
sence of a radioactive source of y-rays 
the liquid always erupted as soon as 
the pressure was released. An-average 
time of 3 seconds was recorded as wait- 
ing time, in the absence of the source. 

The average time between successive 
cosmic ray particles crossing the chamber 


(*) L. BERTANZA, G. MARTELLI and A. 
ZACUTTI: On the theory of the bubble chamber 
(to be published in Nuovo Cimento). 
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has been estimated to be about 2 seconds. 
In this estimate the contribution of the 
electronic component of cosmic rays has 
been neglected, since the experiment has 
been carried out under about 150 g/cm? 
of roof. 

This first. check of the theory of the 
formation of bubbles in a radiation sen- 
sitive liquid confirms the possibility of 
using larger bubble chambers as detectors 
of ionizing particles: further measure- 
ments are now in course, in order to check 
quantitatively the temperature depen- 
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dence of the sensitivity of the device for 
bubbles having different numbers of ele- 
mentary charges. 


The authors wish to express their 
gratitude to Prof. C. BALLARIO, who 
brought their attention to the subject. 
They are deeply indebted to Sig. GIAN- 
LUIGI RapiciorTI for his help in cons- 
tructing the glass parts, and to Sig. 
ARMANDO PECCHI for his invaluable 
cooperation in building the mechanical 
parts of the apparatus. 
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Quantizzazione della massa e formalismo non locale. 


E. MINARDI 


Aosta 


(ricevuto il 29 Aprile 1954) 


Recentemente è stato considerato da 
alcuni autori il problema della quantiz- 
zazione della massa, problema imposto 
dal gran numero di particelle con masse 
differenti recentemente scoperte. 

In questa nota si presenta un esempio 
di quantizzazione delle masse, nel caso 
di particelle neutre, sia con spin 0 che 
con spin 1/2. Come suggerito da WATA- 
GHIN (!) e successivamente da altri autori, 
tale problema è considerato in connes- 
sione con la teoria dei campi non locali 
ed una lunghezza universale ! viene in- 
trodotta nel formalismo. Inoltre tale pro- 
blema è trattato in condizioni stazionarie, 
mediante una condizione supplementare 
già introdotta da YUKAWA (?) e come 
è stato suggerito recentemente da Ray- 
SKI (3). 

Introduciamo le variabili di Yukawa 


(1a) Ue == bo, + oy), 


(1b) Ny = ©, — B, 
Come è stato osservato da RAYSKI 


(*) G. WATAGHIN: Nuovo Cimento, 10, 500 
(1953); 10, 1618 (1953). 

(2?) H. YUKAWA: Phys. Rev., 77, 219 (1950). 

(3) J. RAYSKI: Nuovo Cimento, 10, 1792 
(1953). 


in un punto di vista non locale è molto 
naturale assumere una dipendenza delle 
variabili di campo dalle », oltre che 
dalle x,; poichè in tal caso ambedue le 
forme linearmente indipendenti [y, d,]_ e 
[y, d,]+ nelle y e d, (d, essendo gli ope- 
ratori di spostamento) entrano in gioco 
nella teoria, e non soltanto la prima, come 
nel caso locale. In generale quindi, è 
plausibile che il campo di particelle debba 
essere descritto da un’equazione genera- 
lizzata, somma di due operatori, funzioni 
rispettivamente delle variabili x, e n,. 

Nel caso di una particella scalare o 
pseudoscalare il moto ordinario è de- 
scritto da un’equazione ottenuta dall’in- 
variante pl°p""—— m?, con la ben nota 
sostituzione operatoriale dei p;. 

Assumiamo che le proprietà interne 
della particella pseudoscalare siano de- 
scritte da una funzione g(7,7”) che sod- 
disfa ad un’equazione ottenuta dall’inva- 
riante pp", ove i p sono autovalori 
degli operatori 
(2) a = 

> —I—. 
on, 


Come valore dell’invariante assu- 
miamo 


(3) pe? pl? == mi, 


cosicchè, se (x) è la funzione d’onda 
associata al movimento ordinario, l’equa- 
zione generalizzata è 


a o? 
(4) (a an, … Ve) sO» 
con 

p(x, n) = u(æ)p(n,n") « 


Accanto alla precedente poniamo la 

seguente condizione di Yukawa (?): 

(5) Magn") = 0 - 

Nel sistema del baricentro della parti- 
cella è po? =p! =p;”=0 e per la (5) la 
g(n,n) è ovunque nulla tranne che sul- 
l’iperpiano 7,=0; l'equazione spaziale 
descrivente le proprietà interne diviene 
allora 
(6) A,p(r) + mip(r) = 0, 

con r= (n ++). 

Cerchiamo una soluzione della pre- 
cedente che sia nulla sulla superficie di 
una sfera di raggio 1, ove | è una co- 
stante universale, e con centro nel punto 
r=0. Inoltre essendo 


(7) “a? =%?; X(r)m nPn(T) 


la densità di massa della particella di 
massa totale m,, scegliamo quella solu- 
zione che dia la quantità di massa mag- 
giore nell’elemento di volume dello 8S; 
contenente il baricentro r=0 della par- 
ticella. 

All’esterno della sfera assumiamo la 
soluzione nulla della (6). Perciò le pro- 
prietà interne della particella descritte 
dalla g(r), riguardano, nel suo sistema 
del baricentro, solo il volume dello 8S, 
contenuto entro la sfera. 

Entro tale sfera si hanno le auto- 
soluzioni - 


(8) 


€, essendo una costante di normalizza- 


© 
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zione; lo spettro di masse è: r NEI O. 

i i 

= = 

9 My =\(2n +71) =. È 

(9) n = (2n + 1) ol a 

Nel caso di particelle con spin 1/2 

l’invariante analogo a pp è y,p; vi 

poniamo a 

(10) yp” = — imp. 

L’equazione generalizzata é allora 

lo) © ; < hl : 

EDS sei wr)g(n) = per ie 

m(B + lu(x)p(n) , N 


w(x) essendo lo spinore associato al moto à 
ordinario. EE: 
Nel sistema del baricentro della par- 
ticella, ricordando la (5), la p(y) è diversa 
da zero solo se 7, = 0. di I 
La equazione spaziale descrivente le i 
proprietà interne si scrive 


: E 
(12) 1&, By (N) — m(n) = 0, pe 
(r= Dae De. 

Per risolvere la precedente, introdu- 3 
ciamo, come nella trattazione di Dirac (*) SIN 
di una particella spinoriale in campo cen- _ i 
trale, un operatore SE 
(13) k= B(o'L™ + 1), Roca 
ea 

L” essendo il momento angolare nello ie? 
spazio delle n, e fisici: 


(14) 


ove le componenti dio sono le note ma- Fe 
trici di Pauli. 

Com’é noto, se J(J +1), con J semin- : 
tero dispari positivo, sono gli autovalori 


(4) P. A. M. Drrac: Principles of Quantum 
Mechanics (Oxford, 1947). 
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di (L®+10'), k? ha gli autovalori 
(J+4)?, cosicchè % può essere +1, +2,.... 

Serivendo la (12) in coordinate po- 
lari (3) e ponendo 


374 
(15) p(r) = ale 


essa conduce al sistema di equazioni: 


1 d?Y 247. 
SET oredr 


HO) 7e — 


Si i wy (e+ + 7) ey 
Te 


La soluzione non singolare nell’ori- 
gine della (35a) è la funzione sferica di 
Bessel 


17 Me as 
(17) x(mr) = 


Jx+3(mr) essendo la funzione di Bessel 
di ordine k+4 (k positivo) e ¢ una co- 
stante di normalizzazione. 

Dobbiamo ora porre che Y e Z siano 
nulle sulla superficie della sfera di rag- 
gio l. 

Per la componente Z ciò è ottenuto 
ponendo: 


Jx+3(Mr) , 


1 


(18) c= | my - 


k—1 1+k)Y 
| (4 S40") 
MY r 


Nel caso della Y,(mr) deve invece 
essere soddisfatta la condizione: 


dr 


(19) Jx+3(ml) = 0. 


Questa equazione fornisce uno spettro 
di masse per ogni valore intero positivo 
di k. 

Sia lo spettro di masse (9) delle par- 
ticelle con spin 0, che gli spettri ottenuti 
dalla (22), devono riguardare particelle 
neutre poichè il campo elettromagnetico 
non entra nella presente formulazione. 
Perciò la (9), per n=0, deve dare la 
massa del mesone 7°; assumendo per 
questa il valore 265 m, (m,, massa del- 
l’elettrone), la lunghezza universale 
risulta allora dalla (9) in unità C.G.S.: 


ah 


(20) t= 
2c 265 m, 


12,9 Om cme. 
Con questo valore di / le prime cinque 


radici della (22), nel caso k=1, danno 


mz = 1305 m,, 
ms = 2365 m,. 


m=O, m, =758 m,, 
m,== 1839 m, , 


Nel caso k=2, le prime radici danno 


ie Ol, DOTE 


Mg = 1558 m,, : m= 2072 m,. 

Come si vede non si ottiene la massa 
dell’elettrone e del mesone yu. Il fatto 
che non si ottenga la massa dell’elettrone 
è in accordo con la veduta classica se- 
condo la quale essa è di origine elettro- 
magnetica. Il metodo ora presentato deve 
essere perfezionato in maniera da inclu- 
dere l’interazione tra il campo (7) e il 
campo elettromagnetico e studiare il con- 
tributo di quest’ultimo al valore della 
massa. È quindi in ogni caso molto pre- 
maturo, come si comprende, concludere 
che il presente tentativo possa condurre 
a rappresentare completamente i fatti 
fisici. 


L’autore ringrazia il prof. G. WA-_ 
TAGHIN per le utili discussioni critiche 
sull’argomento. 
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Interaction of 7* Mesons in Photographie Plates (II). 


A. Mineuzz1, G. Puppi and A. RANZI 


Istituto di Fisica dell’ Universita - Bologna 
Istituto Nazionale di Fisica Nucleare - Sezione di Padova 


(ricevuto il 12 Aprile 1954) 


In a sistematic survey of photographic 
plates exposed to the 62 MeV xt meson 
beam of the Nevis cyclotron, the various 
phenomena occurring have been exten- 


great scattering angles have been pointed 
out also in other experiments, with cloud 
chambers (?). The particular angular de- 
pendence, shown by all authors, has sug- 


TABLE: I, 

Elasti Elasti lasti 

Inelastic |Stars plus AR Sa DI 

Stars tt z tt scatt. scatt. scatt. 
PCA SCARE ta 1209400 lg 1990 

No. of events 291 ve 363 125 35 24 
Mean free path (*) 39.9 135.5 30.8 53.9 225:2 369.5 
in Ilford G5 (in em) | + 3.6 + 22.3 + 2.7 + 5.6 + 47.9 | + 101.2 
(*) The (x+ + P) have been excluded. 


sively investigated. We wish to discuss- 


in this letter the problem of the scattering 
by nuclei. In table I the relevant data 
have been collected. 

As is well known from the first re- 
sults of BERNARDINI and cow. (1) both 
elastic and inelastic scattering of pions 
occour in the photographic plates. The 
peculiar characteristics of the inelastic 
scattering with large energy losses and 


(1) G. BERNARDINI, E. T. BoorH and L. 
LEDERMAN: Phys. Rev., 83, 1277 (1951); G. 
BERNARDINI and F. Levy: Phys. Rev., 84, 
610 (1952). 


gested a strong correlation of this phe- 
nomenon with the elementary meson- 
nucleon scattering, but no definite con- 
clusion has been drawn. The aim of this 
paper is to present new experimental 
results on this field as well as a quanti- : 
tative interpretation of the inelastic scat- 
tering by nuclei as due to single scatter- 
ing of the pions against a nucleon inside 
the nuclear matter. 

The experimental results refer to an 


() H. BYFIELD, J. KESSLER and L. LEDER- 
MAN: Phys. Rev., 86, 17 (1952). 
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analysis of 488 events, such as stars, 
inelastic scatterings and elastic scatter- 
ings with a projected angle greater 
than 10°. From these events, 197 show 
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Fig. 1. — Distribution of all scatterings with a 
track length larger then 300u, corrected for 
geometrical losses. 


‘a pion coming out and these latter have 


been analyzed as a function of the para- 
meter AB/E (i.e. the fractional loss of 
energy of the pion). This criterion, in a 
beam of pions not strictly monochromatic, 


_as in our case, permits an analysis of 


the scattering phenomenon that is lar- 


gely independent from the energy of the 
incoming pion. The results are shown 
in Fig. 1, as a block diagram; the heavy 
shaded area refers to the scatterings in 
the backward direction. As it is easy to 
see, our method allows a statistical separ- 
ation of the scattering in two groups: 
an elastic group, symmetrically distri- 
buted with respect to 4H/E=0, and a 
broad inelastic group with practically 
any fractional loss. The differential cross- 
section for the elastic group is shown in 


Fig. 2, where the continuous curve re- 
presents the mean Rutherford scattering, 
i.e. the Coulomb scattering cross-section 
by point nucleus averaged over the 
nuclei of the emulsion for the mean 
energy of the beam. From Fig. 1 statisti- 
cally interpreted, the existence of a group 
of backward directed elastie scatterings 
becomes clearly evident, as shown also by 
other experiments in which the discrim- 
ination between elastic and _ inelastic 
scattering was made with an arbitrary 


cut-off AB/E. 


As emphasized also by others it is 
impossible to describe this group of 
elastic scatterings with a continuous mo- 
del of the nucleus i.e. by a complex po- 


tential of central forces experienced by 
the meson, of the type 


(V+ i8)f(r), 


with V and S real scalar quantities. 
It results from the experiment that if 
the pion-nucleus forces must be built 
up from the pion-nucleon forces, then 
the process of averaging over all the 
nucleons of the nuclei does not lead to 
central forces (3). We have not been able 
to find clearly any effect of interference 
between Coulomb and nuclear potential, 


3 


as shown by LEDERMAN and coworkers 
at the same energy, but for heavy ele- 
ments this is more difficult to see and 
our statistic is not sufficiently high for 
this purpose. The existence of the elastic 
group in the backward direction on the 
contrary is confirmed also by other un- 
published results of our laboratory. 

Let us now discuss the problem of the 
inelastic scatterings, whose angular dis- 
tribution is plotted in Fig. 3. We have 
tried to interpret the inelastic scatterings 
as due to single or plural scatterings of 
the pion against a nucleon inside the 
nuclear matter, supposing that the cha- 
racteristics of these scatterings are the 
same, if allowed by the Pauli principle, 


(3) D. C. PEASLEE: Phys. Rev., 86, 269 (1952); 
K. M. WATSON: Phys. Rev., 89, 575 (1953). 
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as those against the free nucleons. From 
the known characteristics of the meson- 
nucleon scattering (*) we know that the 
most important process, at this energy, 
is the scattering of positive pions on pro- 
tons. A Montecarlo calculation has been 
done for this process inside the nuclear 


matter with a Fermi gas model of the . 


nucleus taking into account the motion 
of nucleons, for his implication on the 


total cross-section (which is strongly © 


energy dependent) and on the kinematics 
of the scattering process. The calculat- 
ions have been performed for a pion’s 
energy of 65 MeV inside the nuclear 
matter; and for this 
57% of the scatterings are allowed by 
the Pauli principle giving a mean cross- 
section of 21 mb per proton, while with- 
out the Pauli principle the mean cross- 
section per proton would be of 27 mb, 
the cross-section assumed for the scat- 


tering against protons at rest being of 
20 mb. 

We have also investigated the in- 
fluence of a real attractive potential 


(4) D. BODANSKY, A. SacHs and J. STEIN- 
BERG: Phys. Rev., 90, 996-997 (1953). 
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V= 15 MeV experienced by the pion 
by crossing the nuclear boundary: this 
potential is required for the explanat- 
ion of the interference effect on lght 
elements. The effect of this potential 
on inelastic scatterings is to increase the 
fractional energy losses of the pion, by 
a factor (14+V/E). Fig. 4 shows the 
calculated angular distribution in the 
lab. system of pions scattered by protons 
ins'de the nucleus together with a con- 
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tinuous curve which represents the scat- 
tering by a proton at rest; the two di- 
stributions are approximately the same. 

The broad distribution of the scat- 


terings as a function of AE/E predicted 
by the calculation (Fig. 5) is also in 
agreement w:th the experimental distri- 
bution if we single out a group around 
the fractional energy loss of AE/E ~0.30 
which is, in part at least, due to the 
scatterings on neutrons. 

Of the total of inelastic scatterings 
experimentally observed 15% roughly 
show large energy losses and forward 
direction; they can be classified as double 
scatterings and have been not included 
in the previous comparison in order to 
isolate the results of the first scattering 
which takes place inside the nuclear 
matter. 

At this stage of the research the most 
reasonable interpretation of the inelastic 
scatterings of pions against nuclei seems 
to be in terms of individual meson nu- 
cleon scattering in which the nucleon 
preserves the properties of interaction of 
the free state, so that non linear effects 
seem to play no important role. 

A much more profound analysis is 
necessary in order to know to what 
degree of approximation this extreme 
model is true. 
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